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Chapter 1. The Orbits of One-Dimensional Maps. 

1.1 Iteration of Functions and Examples of Dynamical Systems. 

Dynamical systems is the study of how things change over time. Examples include 
the growth of populations, the change in the weather, radioactive decay, mixing of 
liquids and gases such as the ocean currents, motion of the planets, the interest in a 
bank account. Ideally we would like to study these with continuously varying time 
(what are called flows), but in this text we will simplify matters by only considering 
discrete changes in time. For example, we might model a population by measuring 
it daily. Suppose that x n is the number of members of a population on day n, where 
x 0 is the initial population. We look for a function / : R -» R, (where R - set of all 
real numbers), for which 

x x = f(x 0 ),x 2 = f(xi) and generally x n = f(x n . 1 ), n = l,2,.... 
This leads to the iteration of functions in the following way: 
Definition 1.1.1 Given x 0 e R, the orbit of x 0 under / is the set 

O(x 0 ) = {x 0 J(x 0 ),f(x 0 ),...}, 
where f 2 (x 0 ) = f(f(x 0 )), f(x 0 ) = f(f 2 (x 0 )), and continuing indefinitely, so that 
f n ( x ) = f ° / ° / ° • • • o f(x); (n-times composition). 
Set x„ = f n {x 0 ), Xi = f(x 0 ), x 2 = f 2 (x 0 ), so that in general 
x n+ i = f n+1 (x 0 ) = f(P( Xo )) = f( Xn ). 

More generally, / may be defined on some subinterval / of R, but in order for the 
iterates of x € I under / to be defined, we need the range of / to be contained in / 
so / : / -4 /. 

This is what we call the iteration of one-dimensional maps (as opposed to higher 
dimensional maps / : R" -> R» n > 1, which will be studied briefly in a later 
chapter) . 

Definition 1.1.2 A ( one dimensional) dynamical system is a function f : I -* I 
where / is some subinterval of R. 

Given such a function /, equations of the form x n+l = f(x n ) are examples of 
difference equations. These arise in the types of examples we mentioned above. For 
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example in biology x n may represent the number of bacteria in a culture after n hours. 
There is an obvious correspondence between one-dimensional maps and these differ- 
ence equations. For example, a difference equation commonly used for calculating 
square roots : ^ iv^ew €Y *V 

*f 2- x n+1 = -(x n + — ), 

corresponds to the function f(x) = \(x + | ). If we start by setting x 0 = 2 (or in 
fact any real number), and then find Xi, x 2 etc., we get a sequence which rapidly 
approaches y/2 (see page 9 of Sternberg [49] ) . One of the issues we examine is what 
exactly is happening here. 

Examples of Dynamical Systems 1.1.3 

1. The Trigonometric Functions Consider the iterations of the trigonometric 
functions starting with / : R — > R, f(x) = sin(x). Select x 0 G R. at random, e.g., 
Xq = 2 and set x n +\ — sin(x* n ), n = 0, 1, 2, — Can you guess what happens to x n as 
n increases? One way to investigate this type of dynamical system is to enter 2 into 
our calculator, then repeatedly press the 2nd, Answer and Sin keys (you will need 
to do this many times to get a good idea - it may be easier to use Mathematica, or 
some similar computer algebra system). Do the same, but replace the sine function 
with the cosine function. How do we explain what appears to be happening in each 
case? These are questions that we aim to answer quite soon. 

2. Linear Maps Probably the simplest dynamical system (and least interesting from 
a chaotic dynamical point of view), for population growth arises from the iteration of 
linear maps: maps of the form f(x) = a • x. Suppose that x n = size of a population 
at time n, with the property 

x n _|_i — a - x n 

for some constant a > 0. This is an example of a linear model for the growth of the 
population. 

If the initial population is Xo > 0, then xi = ax 0 , x 2 = axi = a 2 x 0 , and in general 
x n = a n xo for n = 0, 1,2, . . .. This is the exact solution (or closed form solution) 
to the difference equation x n+ i = a • x n . Clearly f(x) = ax is the corresponding 
dynamical system. We can use the solution to determine the long term behavior of 
the population: 

x n is very well behaved since: 

(i) if a > 1, then x n — > oo as n — » oo, 



(ii) if 0 < a < 1 then x„^Oasn^oo (i.e., the population becomes extinct), 
(hi) if a = 1, then the population remains unchanged. 

3. Affine maps These are functions / : M — ► K of the form)7(x) =~aV+b\(a ^ 0), 
for constants a and b. Consider the iterates of such maps: 



f(x) = f(ax + b) = a(ax + b) + b = a 2 x + ab + b, 
f{x) = a 3 x + a 2 b + ab + b, 
f 4 (x) = a 4 x + a 3 b + a 2 b + ab + b, 

and generally 

r(x) = a n x + a n - l b + a n -% + ■ • • + ab + b. 
Let x Q 6 K and set x n = f n (x 0 ), then we have 

x n = a n x 0 + (a"" 1 + a n ~ 2 + ■ . . + a + 1)6 



= a"x 0 + 6 

or 



'a n -l\ 



1 



Xn=(xo + ^)a« + r A-l i fa^l ) 



is the closed form solution in this case (here we have used the formula for the sum of 
a finite geometric series: 

n-l 

when r ^ 1). If a = 1, the solution is x n = x 0 + nb . 

We can use these equations to determine the long term behavior of x n . We see 
that: 

(i) if \a\ < 1 then a n -> 0 as n -> oo, so that 

b 



lim x n = , 

n-+oo 1 — a 

(ii) if a > 1, then lim^^ x n = oo for 6, x 0 > 0, 

(iii) if a = 1, then li m7WOO x n = oo if b > 0. 
The limit won't exist if a < -1. 
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4. The Logistic Map : R — > R, ^(^) = - x) was introduced to model 
a certain type of population growth (see [29]). Here \i is a real parameter which 
is fixed. Note that if 0 < // < 4 then is a dynamical system of the interval 
[0, 1], i.e. L M : [0, 1] -+ [0, 1]). For example, when /i = 4, L 4 (x) = 4x(l - x), with 
^([0, 1]) = [0, 1] with graph given below. If fi > 4, then L M is no longer a dynamical 
system of [0, 1] as Z/ M ([0, 1]) is not a subset of [0, 1]. 




0.0 0.2 0.4 0.6 0.8 1.0 



Recurrence Relations 1.1.4 

Many sequences can be defined recursively by specifying the first term (or two), 
and then stating a general rule which specifies how to obtain the nth term from the 
(n — l)th term (or other additional terms), and using mathematical induction to see 
that the sequence is ik well defined" for every n G Z + = {1,2,3, . . .}. For example, 
n\ — n-factorial can be defined in this way by specifying 0! = 1, and n! = n • (n — 1)!, 
for ii e Z+. The Fibonacci sequence F n , can be defined by setting 

F 0 = l, Fi = l and F n+ i = F n + F n - U for n G Z + , 

so that F 2 = 2, F 3 = 5 etc. 

The orbit of a point x 0 G K under a function / is then defined recursively as follows: 
we are given the starting value x 0 , and we set 

Zn = /(Zn-i), for nGZ + . 

The principle of mathematical induction then tells that x n is defined for every n => 0, 
since it is defined for n = 0, and assuming it has been defined for k = n — 1 then 
x n = f(x n -i) defines it for k = n. 
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Ideally, given a recursively defined sequence x n , we would like to have a specific 
formula for x n in terms of elementary functions (so called closed form solution), but 
this is often very difficult, or impossible to achieve. In the case of affine maps' and 
certain logistic maps, there is a closed form solution. One can use the former to study 
problems of the following type: 

Example 1.1.5 An amount $ T is deposited in your bank account at the end of each 
month. The interest is r% per period. Find the amount A(n) accumulated at the 
end of n months (Assume A(0) = T). 

Answer. A(n) satisfies the difference equation 

A(n + 1) = A(n) + A(n)r + T, where A(0) = T, 



or 



A(n + l) = A(n)(l + r) + T, 

so we set x 0 = T, a = 1 + r and b = T in the formula of Example 1.1.3, then the 
solution is 

A(n) = (1 + r) n T + T ( (1+r) "~M 

V 1 + r - 1 ) 

= (l + r) n T+j((l + r ) n -l). 

Remark 1.1.6 It is conjectured that closed form solutions for the difference equation 
arising from the logistic map are only possible when n = ±2,4 (see Exercises 1.1 # 
3 for the cases where n = 2, 4, and # 7 for the case where fi = -2 and also [53] for a 
discussion of this conjecture). 



Exercises 1.1 

1. If L^x) = fi X {l - x) is the logistic map, calculate Lj(z) and L 3 M (x). 

2. Use Example 1.1.3 for affine maps to find the solutions to the difference equations: 

X 

(i) x n+1 - y = 2, x 0 = 2, 



11 



(ii) x n +i + 3x n = 4, x 0 = -1. 



3. A logistic difference equation is one of the form x n +i = fix n (l - x n ) for some fixed 
/i6R. Find exact solutions to the logistic equations: 

(i) x n +i = 2x n (l - x n ). Hint: use the substitution x n = (1 - y n )/2 to transform the 
equation into a simpler equation that is easily solved. 

(ii) x n+ i = 4x n (l - x n ). Hint: set x n = sin 2 (# n ) and simplify to get an equation that 
is easily solved. 



4. You borrow $ P at r % per annum, and pay off $ M at the end of each subsequent 
month. Write down a difference equation for the amount owing A(n) at the end of 
each month (so A(0) = P). Solve the equation to find a closed form for A(n). If 
P = 100, 000, M = 1000 and r = 4, after how long will the loan be paid off? 



5. If Tpfe) = | x j If 2/2 < x </l 3 s ^ low ^ iat * s a dynamical system of 

[0, 1] for fi e (0, 2] 



6. Let / : R — > R be the function defined below. For each of the intervals given, 
determine whether / can be considered as a dynamical system /:!—►/: 

(a) f(x) = x 3 - 3x, 

(i) / = [-ML (n) /=["2,2]. 

(b) f(x) = 2x 3 - 6x, 

(i) /=[-!,!], (") /= [-Vf'\/| ] ' /= [- 4 ' 4 ]- 



7. For the following functions, find / 2 (x), / 3 (x) and a general formula for f n (x): 

mm-*, («)/(*) = i* (in) /w - { , » <Jl 
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8. Use mathematical induction to show that if f(x) = 2 then 

x + 1 ' 

r , ) = 2"(x + 2) + (-l)"(2x-2) 
2 n (ar + 2)-(-l)»(ar-l) ' 



9*. Show that a closed form solution to the logistic difference equation when 
is given by 

x n = \[^-f[r n r 1 (l-2x 0 )]], where r = -2 and f(0)=2co B ' 



(Hint: Set x n = — 2- an d use steps similar to 3(H)). 



2 
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1.2 Newton's Method and Fixed Points 

Given a differentiate function / : R — > R, Newton's method often allows us to 
find good approximations to zeros of /(x), i.e., approximate solutions to the equation 
f(x) — 0. The idea is to start with a first approximation xq and look at the tangent 
line to f(x) at the point (#o, /(#o))- Suppose this line intersects the x-axis at Xi 3 
then we can check that 



For reasons that we will make clear, it is often the case that X\ is a better approxima- 
tion to the zero x = c than xq was. In other words, Newton's method is an algorithm 
for finding approximations to zeros of a function f(x). The algorithm gives rise to a 
difference equation 

_ „. /O^) 

where Xo is a first approximation to a zero of f(x). The corresponding real function 
is 



X i = x 0 — 



/(so) 



if /'(*o)^0. 



iV/(x) = x - 



/(*) 



(the Newton function). 



/'(*)' 



Note that if f(x) = x 1 — a, then f(x) = 2x and 




so that 




the difference equation we mentioned in Section 1.1 that is used for approximating 
y/2 when a = 2. 



L4 




Note that we are looking for where /(*) = 0 , and this happens if and only if 

N f (x) = x. J 

Definition 1.2.1 For a function / : R - R, a point c E R for which /(c) = c is called 
a fixed point of /. It is a point where the graph of f(x) intersects the line y = x. We 
denote the set of fixed points of / by Fix(/). 



/ 


/\ 




/ \ 


/ 













Fixed points occur where the graph of f(x) intersects the line y = x. 

Example 1.2.2 Suppose that f(x) = x\ then x* = x gives x(x - 1) = 0, so has fixed 
points c = 0 and c = 1, so Fix(/) = {0, 1}. If f(x) = x* - x, then x* - x = x gives 
x(x - 2) = 0, so the fixed points are c = 0 and c = ±y/2, Fix(/) = {0, ±V2}. 
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The logistic map f(x) = 4x(l - x) = 4x - 4x 2 , 0 < x < 1 has the properties: 
/(0) = 0, (a fixed point), f(l) = 0, a maximum at x = 1/2 (with /(1/2) = 1). 
Solving /(x) = x gives 4x - 4x 2 = x, so 4x 2 = 3x, so the fixed points are c = 0 and 
c - 3/4. 

This map has what we call eventual fixed points: /(l) = 0 and /(0) = 0, so we say 
that c = 1 is eventually fixed. Also /(1/4) = 3/4, so c = 1/4 is eventually fixed, as is 
c = (2 + >/3)/4. 

Definition 1.2.3 x* G R is an eventual fixed point of f(x) if there exists a fixed point 
c of f(x) and r G Z + satisfying / r (x*) = c, but / 5 (x*) 7^ c for 0 < s < r. 

Example 1.2.4 The Tent Map 

Define a function T : [0, 1] — > [0, 1] by 



T(x) is called the ieni map, it has fixed points: T(0) = 0 and T(2/3) = 2/3. Since 
T(l/4) = 1/2, T(l/2) = 1 and T(l) = 0, 1/4, 1/2 and 1 are eventually fixed. It is not 
difficult to see that there are many other eventually fixed points (infinitely many). 



T(x) 




2x :0<x<l/2 
2(1 -x) :1/2<x<1 



= l-2|x-l/2|. 



1.0 




0.8 



0.6 



0.4 



0.2 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



Note that some maps do not have fixed points: 
Example 1.2.5 f(x) = x 2 + 1 has no fixed points. 
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Question 1.2.6 For what values of a £ R does Q a {x) = x 2 + a have fixed points? 
We see below that certain functions always have fixed points: 

Theorem 1.2.7 Let f : I -> / be a continuous function, where I = [a,b], a < b is a 
bounded interval. Then f(x) has a fixed point c € /. 

Proof. Set g(x) = f(x) - x. We may assume that f(a) ? a and f(b) f b, so that 
/(a) > a and f(b) < b. 




The graph of f(x) always intersects the line y = x. 

It follows that 

9(a) = f(a) - a > 0, and g(b) = f(b) - b < 0, 
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so that g(x) is a continuous function which is positive at a and negative at 6, so by 
the Intermediate Value Theorem (IVT), there exists c G (a, b) (open interval) with 
g(c) = 0, i.e., /(c) = c, so c is a fixed point of f(x). □ 

Remark 1.2.8 The above is an example of an existence theorem, it says nothing 
about how to find the fixed point, where it is or how many there are. It tells us that 
if f(x) is a continuous function on an interval / with /(/) C /, then f(x) has a fixed 
point in /. It is also true that if instead f(I) D I then f(x) has a fixed point in I as 
we see now: 

Theorem 1.2.9 Let f : / — > M (J = [a, 6], a < b) be a continuous function with 
f(I) 5 I, then f(x) has a fixed point in I. 

Proof. As before, set g(x) = f(x) - x, then there exists C\ G (a, b) with /(ci) < C\ 
(in fact f(ci) < a < Ci). Also there is c 2 G (a, 6) with /(c 2 ) > c 2 . 



Then ^(ci) < 0 and #(c 2 ) > 0 and since g(x) is a continuous function, it follows by 
the Intermediate Value Theorem that there exists c G /, (c\ < c < c 2 or c 2 < c < Ci), 
with 0(c) = 0; /(c) = c. □ 

Remark 1.2.10 It is often difficult to find fixed points explicitly: 

Example 1.2.11 Set f(x) = cosx. If /(c) = c, then cose = c. It is possible to find 
an approximation to the fixed point c = .739085... using (for example) Newton's 
method. 
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f(x) = cosx has a fixed point in [0, tt/2]. 



Exercises 1.2 



1. Find all fixed points and eventual fixed points of the map f(x) = \x - 1|. 

2. Let / : R -» R be such that for some n e Z+, the nth iterate of / has a unique 
fixed point c (i.e., /"(c) = c and c is unique). Then show that c is a fixed point of /. 

3. Use the Mathematica command NestList (see below) to find how the following 
functions behave when the given point is iterated (comment on what appears to be 
happening in each case): 

(i) Li(x) = x(l - x), with starting point x 0 = .75. 

(ii) L 2 (x) = 2x(l-x), with starting point x 0 = .1. 
(in) L 3 (x) = 3x(l - x), with starting point x 0 = .2. 

(iv) L 3 . 2 (x) = 3.2x(l-x), with starting point x 0 = .95. 

(v) f(x) = sin(x), with starting point x Q = 9.5. 

(vi) g(x) = cos(x), with starting point x 0 = -15.3. 

4. Show that the logistic map L^x) = (jlx{\ - x), x G [0, 1], for 0 < // < 4, the fixed 
point x = 0 has no eventual fixed points other than x = 1. 

Show that for 1 < (x < 2, the fixed point x = 1 - 1/ M only has x = 1/// as an 
eventual fixed point, but this is not the case for 2 < fj, < 3. 

To use Mathematica,, first define your function, say f(x) = sin(x) in an input 
cell using the syntax 

f [x_] := Sin[x] . 

Execute the cell by doing Shift and Enter simultaneously. The command 

NestList [f, 9.5, 100] 
when executed in a new input cell will give 100 iterations of sin(rc) with starting 
value x 0 = 9.5. If you want to graph f{x) on the interval [a, b], use the command: 

Plot[f[x], { x, a, b } ]. 
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5. (a) Let f(x) = (1 + x)" 1 . Find the fixed points of / and show that there are 
no points c with f 2 (c) = c and f(c) ^ c (period 2-points). Note that f(-l) is not 
defined, but points that get mapped to -1 belong to the interval [-2, -1) and are of 
the form 



"n = =— , n>l, 



where {F n }, n > 0, is the Fibonacci sequence (see 1.1.4). Note that v n -> -r as 
n -> oo, where -r is the negative fixed point of / (see [7] for more details). 

F x ~\~ F 

(b) If x 0 e (0, 1], set x n = — Use mathematical induction to show that 

r n x 0 + t n+ i 



In+1 = fM = f " x i +f ™ 



Deduce that asn-^oc, x n — >l/r, the positive fixed point of f. 



1.3 Graphical Iteration 

It is possible to follow the iterates of a function f(x) at a point x 0 using graphical 
iteration (sometimes called web diagrams). * ^b^Jeb r <?( £-\a\.Y g 

We start at x 0 on the x-axis and draw a line vertically to the function. We then 
move horizontally to the line y = x, then vertically to the function and continue in this 
way. Notice that in some examples the iterations converges to a fixed point. In others 
it goes off to oo, whilst in others still, it oscillates between two points indefinitely. 

Example 1.3.1 f(x) = x(l - x). In this example, an examination of graphical 
iteration seems to suggest that the orbits of any point in [0, 1] approach the fixed 
point x = 0. This is an example of an attracting fixed point. 
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Example 1.3.2 Let f(x) = x/2 + 1. This is an affine transformation with a = 1/2, 
b=l. According to what we saw in Example 1.1.3, the iterates should converge to 
jz^ = 2 since \a\ < 1. What is actually happening is that c = 2 is an attracting fixed 
point of f(x) with the property that it attracts all members of R (said to be globally 
attracting. 




Example 1.3.3 We see from the graph of f(x) = sinx that c = 0 seems to be an 
attracting fixed point. We shall show later that it is actually globally attracting. 
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We see two basic situations arising (together with some variations of these). In 
particular, we notice that if the sequence x n = f n (x 0 ) converges to some point c as 
n —» oo, then c is a fixed point. 

(i) Stable orbit, where the web diagram approaches a fixed point. 

(ii) Unstable orbit, where the web diagram moves away from a fixed point. 




Proposition 1.3.4 If y = f(x) is a continuous function and lining f n (x Q ) = c, 
then f(c) = c (i.e., if the orbit converges to a point c, then c is a fixed point of f) . 

Proof. We see that 

Km f n (x 0 ) = c /(Um r(x 0 )) = /(c), 
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and since / is continuous, 

to /•«(>>)« /(c). 

But clearly c = lining f n+1 (x 0 ), so that c = /(c) by the uniqueness of limit. □ 

1.4 Attractors and Repellers 

Suppose that / : X -> X is a function (where we assume that X is some subset of 
R such as an interval or possibly R itself). Let c e X be a fixed point of /. We define 
the notions on attracting fixed point and repelling fixed that we looked at graphically 
in the last section. 

Definition 1.4.1 Let / : X -* X with /(c) = c. . / ,( V^,^, ^ ?A 

(i) c is a ^ble fixed jpoint if for all e > 0 there exists (5 > 0 such that if x € X and 
Is - c| < 5, then \ f n {x) - c\ < e for all n £ Z+. 

If this does not hold, c will be called unstable. 

(ii) c is said to be attracting if there is a real number n > 0 such that 

|s - c| < 77 =>■ lim / n (rr) = c. 
n— »oo 

(iii) c is asymptotically stable if it is both stable and attracting. 

Remark 1.4.2 We will show that a fixed point c of / : X -> X (suitably differ- 
entiate) with the property \f'(c)\ < 1 is an asymptotically stable fixed point. If c 
is an unstable fixed point, we can find e > 0 and x arbitrarily close to c such that 
some iterate of x, say f n {x), will be greater than distance e from c. This is the case 
when | /'(c) | > 1 as the next theorem shows. We call c a repeller (c is a repelling 
fixed point), since the iterates move away from the fixed point (c is unstable). We will 
also see that a fixed point can be stable without being attracting, and that it can be 
attracting without being stable. 

Definition 1.4.3 A fixed point c of / is hyperbolic if |/'(c)| ^ 1. If |/'(c)| = 1 it is 
non-hyperbolic. 

Thus if c is a non-hyperbolic fixed point, then /'(c) = 1 or /'(c) = -1, so the graph 
of f(x) either meets the line y = x tangentially, or at 90°: 



is aw open \Mtr**\ 
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Up 1 



U4\ 



/(#) = — 2x 3 + 2x 2 + x has both types of non-hyperbolic fixed points.' 



Hyperbolic fixed points have the following properties: 



Theorem 1.4.4 Let f : X — > X be a differentiable function with continuous first 
derivative (we say that f is of class C 1 ). 

(i) If a is a fixed point for f(x) with < 1, then a is asymptotically stable. 
The iterates of points close to a, converge to a geometrically (i.e., there is a constant 

0 < A < 1 for which \f n (x) — a\ < X n \x — a\ for all n G Z + and for all x G X 
sufficiently close to a). 

(ii) If a is a fixed point for f(x) for which |/'(a)| > 1, then a is a repelling fixed point 
for f. 

Proof, (i) We may assume that X is an open interval with a G X. Suppose \f(a)\ < 
A < 1 for some A > 0, then using the continuity of f'(x), there exists an open interval 

1 with < A < 1 for all xel. 

By the Mean Value Theorem, if x G / there exists c G /, lying between x and a, 
satisfying 



/'(c) = 



/(*) - /(«) 



x — a 



so that 

\f(x)-a\ = \f'(c)\\x-a\<\\x-al 
(i.e., f(x) is closer to a than x was). 

Repeating this argument with f(x) replacing x gives 

A \f(x) -a\< X 2 \x - a|, . . . , \f n {x) - a| < X n \x - a\. 



> - 
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Since X n -> 0 as n -» oo, it follows that / n (x) -> a as n -> oo. 
The proof of (ii) is similar. r-. 

Example 1.4.5 Denote the Logistic map by L^x) = fix(l - x) = fix - fix 2 We 
are interested in this map for 0 < x < 1 and 0 < fi < 4, since for these parameter 
values fx, is a function which maps the interval [0, 1] into the interval [0, 1], so it 
necessarily has a fixed point and iterates stay inside [0, 1] (if fj, > 4, L^x) > 1 for 
some values of x in [0, 1] and further iterates will go to -oo). 

Solving L,(x) = x gives x = 0 or x = 1 - 1/fi. If 0 < n < 1, then 1 - l//x < 0, so 
c = 0 is the only fixed point in [0, 1] in this case. 

For 0 < fi < 1, L'^x) =fi- 2fix, so L;(0) = fx and 0 is an attracting fixed point 
for 0 < fL < 1. 




/i > 1, then 0 and 1 - 1/// are both fixed points in [0, 1], but now 0 is repelling. 
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Also, 




- l//i) = /i- 2/i(l - = 2 - /i, 

so that 

|l£(l-l/M)IH2-Ji|<l iff K/i<3. 

In this case c = 1 - l//i is an attracting fixed point, and repelling for \i > 3. Note 
that L[(0) = 1 and £3(2/3) = 1 so we have non-hyperbolic fixed points. 

We will examine in Chapter 2 what happens in the cases where fi = 1 and /i = 3. 

Example 1.4.6 Let f(x) = l/x. The two fixed points x = ±1 are°hyperbolic since 
f'(x) = -l/x 2 , so /'(±1) = -1. However, they ar e stable but not attracting since 
P( x ) = /(l/^) = and we see points close to ±1 neither move closer nor further 
apart. 



4 ^ 

: ! 

1 

3 - 1 

\ 
\ 

I - 






- 











0 12 3 4 
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Example 1.4.7 Consider instead the function f a (x) = 



—2x if x < a 



0 



if x > a 



, where 



a is an y positive rea l number. Then c = 0 is an unstable (repelling) fixed point of 
f a which is attracting i.e., lim^oo f n (x) = 0 for all x e R. A fixed point having 
this latter property is said to be globally attracting . Note that the function f a is not 
continuous at x = a. It has been shown by Sedaghat [44] that a continuous mapping 
of the real line cannot have an unstable fixed point that is globally attracting. 





r 

6) <-?'J>~<^ 



The map /i (in bold) where x = 0 is an attracting but unstable fixed point. 

£(*) tofc£c*4rtVV>*S 



Example 1.4.8 Newton's Method Revisited. 

Suppose that /:/—>/ is a function whose zero c is to be approximated using 
Newton's method. The Newton function is N f (x) = x - f(x)/f(x), where we are 
assuming /'(c) ^ 0. Notice that since /(c) = 0, N f (c) = c, i.e., c is a fixed point of 
N f . Consider N' f (c): 

M'( x ) = , _ /WW ~ /(*)/"(*) _ /(*)/"(*) 



so that 



since f(c) = 0. 



Nl{c) = mnc ) = 
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It follows that \N' f (c)\ = 0 < 1, so that c is an attracting fixed point for N f , and 
in particular 

lim A7(x 0 ) = c, 

provided x 0 , the first approximation to c, is sufficiently close to c. 

Definition 1.4.9 A fixed point c for f(x) is called a super- attracting fixed-point if 
/'(c) = 0. This gives a very fast convergence to the fixed-point for points nearby. 

Remark 1.4.10 Suppose that /'(c) = 0, then N f (x) is not defined at x = c, since the 
quotient f(x)/f'(x) is not defined there. Suppose we can write f(x) = (x - c) k h(x) 
where h(c) ^ 0 and k G Z+ (for example if /(x) is a polynomial with a multiple root), 
then we have 

/(*) = ~ c)*ft(s) (x-c)/i(x) 

/'Or) fc(x - c)*" 1 ^) + (x - c)^'(x) + (x - c)/i'(x) " 3 

when x = c, so that iV/(x) = c has a removable discontinuity at x = c (removed by 
setting N f (c) = c). Then we can find N' f (x) and show that N f f (x) has a removable 
discontinuity at x = c, which can again be removed by setting N' f (c) = (k - l)/k, so 
that \N' f (c)\ < 1 (see the exercises). 

We summarize the above with a theorem: 

Theorem 1.4.11 For a differ entiable function f : I —> I, the zero c of f(x) is a 
super- attracting fixed point of the Newton function N f , if and only if f(c) ^ 0. 

Example 1.4.12 Suppose f(x) = x 3 - 1, then /(l) = 0 and 

so N f (l) = 1 and N' f (x) = § - so that N' f (l) = 0. We see from graphical iteration 
(next section), very fast convergence to the fixed point of Nr. 
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Very fast convergence to the fixed point. 
Exercises 1.4 

1. Find the fixed points and determine their stability for the function f(x) = 4 - ^. 

2. For the family of quadratic maps Q c (x) = x 2 + c, x G [0, 1], use the Mathemat- 
ica program webPlot below to give graphical iteration for the values shown (use 20 
iterations): 

(i) c = 1/2, starting point x 0 = 1, (ii) c = 1/4, starting point x 0 = .1, (iii) c = 1/8, 
starting point x 0 = .7. 

Mathematica Program: Use as in: 
webPlot[Cos[x] , {x, 0, Pi}, {.1, 20}, AspectRatio->Automatic] 

webPlot[g_, {x_, xmin_, xmax_}, {al_, n_}, opts_] := 
Module [{seq, r, pts, web, graph}, r[t_] := N[g /. x -> t] ; 
seq = NestList[r, al, n] ; pts = Flatten[ 
Table [{{seq [[i]], seq [[i + l]]},{seq[[i + l]],seq[[i + l]]}},{i, 1, n}] , 1] 
web = Graphics [{Hue [0] , Line [PrependTo [pts, {seq[[l]], 0}]]}]; 
graph = Plot[{x, r[x]}, {x, xmin, xmax}, 

DisplayFunction -> Identity]; Print[ M last iterate =", Last[seq]]; 
Show [web, graph, opts, Frame -> True, 



PlotRange -> {{xmin, xmax}, {xmin, xmax}}]] 
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3. Let S^x) = /isin(x), 0 < \i < 2tt, 0 < /i < 7r and C^(x) = /icos(x), -n < x < it 
and 0 < \± < 7i. 

(a) Show that has a super-attracting fixed point at x = 7r/2, when \i = tt/2. 

(b) Find the corresponding values for having a super-attracting fixed point. 

2 

4. Show that the map f(x) = has no periodic points of period n > 1 (Hint: 

x + 1 

Use the closed formula in Exercise 1.1.6). 

5. Let N f be the Newton function of the map f(x) = x 2 + 1. Clearly there are no 
fixed points of the the Newton function as there are no zeros of /. Show that there 
are points c where Nj(c) = c (called period 2-points of Nf). 

6. (a) Suppose that f(c) = f'(c) = 0 and /"(c) ^ 0. If f"(x) is continuous at x = c, 
show that the Newton function iV/(x) has a removable discontinuity at x = c (hint: 
apply L'Hopital's rule to Nf at x = c). 

(b) If in addition, / 7// (^) is continuous at x = c with /"(c) ^ 0, show that iV}(c) = 1/2, 
so that x = c is not a super- attracting fixed point in this case. 

(c) Check the above for the function f(x) = x 3 - x 2 with c = 0. 

7. Continue the argument in Remark 1.4.11 (generalizing the last exercise), to show 
that the derivative of the Newton function N f has a removable discontinuity at x = c, 
which can be removed by setting Nj(c) = (k - l)/k. 

8. Let / be a twice differentiate function with f(c) = 0. Show that if we find the 
Newton function of g{x) = f(x)/f'(x), then x = c will be a super-attracting fixed 
point for N g , even if f(c) = 0 (this is called Halley's method). 
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1.5 Non- hyperbolic Fixed Points. 

Example 1.5.1 We have seen that it / : X X (where usually X is some subinterval 
of R) and a € X with /(a) = a, \f'(a)\ < 1, with suitable differentiability conditions, 
then a is a stable fixed point for /. This term is used because points close to a will 
approach a under iteration. If we look at a function like f(x) = sinx, we see that 
c = 0 is a fixed point but in this case |/'(0)| = 1, so it is a non-hyperbolic fixed point. 
However, graphical iteration suggests that the basin of attraction of / is all of R, so 
c = 0 is a stable fixed point. Before considering non-hyperbolic fixed points in more 
detail, let us prove the last statement analytically: 

Proposition 1.5.2 The fixed point c = 0 of f(x) = sinx is globally attracting . 

Proof. First notice that c = 0 is the only fixed point of sinx. This is clear for if 
sinx = x then we cannot have \x\ > 1 since | sinx| < 1 for all x. If 0 < x < 1, the 
Mean Value Theorem implies there exists c e (0, x) with 

/'(c) = f ^ ~ = sinx 
x-0 x ' 

so that 

sinx = xcosc < x, 
since | cosc| < 1 for c € (0, 1). Similarly if -1 < x < 0. 

To show that /(x) is globally attracting, let x e R. We may assume that -1 < 
x < 1, since this will be the case after the first iteration. 

Suppose that 0 < x < 1, then we note that 0 < /'(*) < 1 on this interval. By the 
Mean Value Theorem, there exists c e (0, x) with 

tu \ f(x) - /(0) 

f W = or 0<f(x) = f'(c)x<x. 

Continuing in this way, we see that 

0<r(x)</"- 1 (x)<...</( 2; )< a ;, 

so we have a decreasing sequence x n = f n (x) bounded below by 0. It follows that this 
sequence converges. But Proposition 1.4.5 implies that if this sequence converges, it 
must converge to a fixed point, c = 0 being the only fixed point gives f n (x) -> O as 
n — > oo. A similar argument can be used for -1 < x < 0. □ 

Example 1.5.3 It is possible for the fixed point to be unstable, but to have a one- 
sided stability (to be semi-stable). For example, consider /(x) = x 2 + 1/4 which has 
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the single (non- hyperbolic) fixed point c = 1/2. This fixed point is stable from the 
left, but unstable on the right. 




In the following we give some criteria for non-hyperbolic fixed points to be asymp- 
totically stable/unstable etc. It also gives a criteria for semi-stability. 

Theorem 1.5.4 Let c be a non-hyperbolic fixed point of f(x) with /'(c) = 1. If f'(x), 
f"(x) and f'"(x) are continuous at x = c, then: 

(i) if f "(c) ^ 0, the c is semi-stable, 

(ii) if /"(c) = 0 and f'"(c) > 0, then c is unstable, 

(hi) If f"{c) = 0 and f"'(c) < 0, then c is asymptotically stable. 

Proof, (i) If /'(c) = 1 then f(x) is tangential to y = x at x = c. Suppose that 
f"(c) > 0, then f(x) is concave up at x = c and the picture must look like the 
following: 
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We see this gives stability on the left and instability on the right. 

More formally, since the derivatives are continuous, and /"(c) > 0 this will be 
true in some small interval (c-S,c + S) surrounding c. In particular, the derivative 
function f'(x) must be increasing on that interval, so that since f'(c) = 1 then 

f'(x)<l for all x€(c-S,c), and f'(x) > 1 for all xe(c,c + S), 

for some 6 > 0. Also, from the continuity of /'(*), we can assume that fix) > 0 in 
this interval. 

Now by the Mean Value Theorem applied to the interval [x, c] C (c - S c) there 
exists q € (x,c) with ' J ' 



x — c 

Now since 0 < f(q) < 1 and c > x, we have 

o<MzM <1 

X — c 

or 

X < f( x ) < c. 



Repeating this argument, we see that the sequence /»(*) is increasing and bounded 
above by c, so must converge to a fixed point. There can be no other fixed point (say 
d ? c), m this interval as the Mean Value Theorem would give f(q x ) = l for some 
</i G (x, c), a contradiction. Consequently we see that c is stable on the left 

On the other hand, if [c, *] c (c, c + 6), then applying the Mean Value Theorem as 
above gives 

f'(o) - ~ f ^ ^ i ti ^ 

J {q) X ~Z~ C — > X ' so f(x) >x>c, 

since x - c> 0. This tells us that the point moves away from c under iteration so 
the fixed point is unstable on the right, Similar considerations can be used when 
/ (c) < 0 and the graph is concave down at x = c. 

(ii) is similar to (iii), so is omitted. 

(iii) In this case /'"(c) < 0, /"(c) = 0 and f'(c) = 1. We will show that we have a 
point of inflection at x = c as in the following picture: 
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0.5 - 



By the second derivative test, f'(x) has a local maximum at x = c (the continuous 
function f'(x) is concave down). It follows that 

f'(x) < 1 for all x e (c - 5, c + 6), x ^ c, 

for some 6 > 0 (f"(x) > 0 for x e (c-8, c), so f'(x) is increasing there, and f"{x) < 0 
for x G (c, c + 6), so /'(x) is decreasing there). In particular f(x) ^Wix + c. 
Now use an argument similar to that of (i) above to deduce the result. □ 

Example 1.5.5 Returning to the function f(x) = sin re, we see that /'(O) = 1, 
/"(0) = 0 and f"'(0) = -1, so the conditions of Theorem 1.5.4 (iii) are satisfied and 
we conclude that x = 0 is an asymptotically stable fixed point. 

If f(x) = t&nx, then /'(0) = 1 and we can check that Theorem 1.5.4 (ii) holds 
so that the fixed point x = 0 is unstable. For f(x) = x 2 + 1/4, with /(1/2) = 1/2, 
/'(1/2) = 1, we can apply Theorem 1.5.4 (i). 

How do we treat the case where /'(c) = -1 at the fixed point? We use: 

Definition 1.5.6 The Schwarzian derivative Sf(x) of f(x) is the function 



Sf(x) = 



/"'(*) 
/'(*) 



3 
2 



7"(aQ 
./'(*) J 



so that 



Sf(x) = -f"'{x) - |[/"0r)] 2 , when /'(*) = -1. 
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Theorem 1.5.7 Suppose that c is a fixed point for f(x) and /'(c) = -1. // f'(x), 
f"(x) and f"'(x) are continuous at x = c then: 

(i) tf Sf{c) < 0, then x = c is an asymptotically stable fixed point, 

(ii) if Sf(c) > 0, then x = c is an unstable fixed point. 

Proof, (i) Set g(x) = f 2 (x), then g(c) = c and we see that if c is asymptotically 
stable with respect to g, then it is asymptotically stable with respect to /. Now 

so that g'{c) = /'(c) • /'(c) = (-1)(-1) = 1. 

The idea is to apply Theorem 1.5.4 (iii) to the function g(x). Now 

Ax) = fV(x))-f"(x) + f»(f(x)).[f'(x)]>, 

thus 

9"{c) = /'(c)/"(c) + /"(c)[/'(c)] 2 = 0, since f(c) = -1. 

Also 

9''\^-f\myf\x).f%x) + f\f(x)).f\x) + r%f(xw 

Therefore 

9"'{c) = [/"(c)] 2 (-l) - /'"(c) - /'"(c) + 2/"(c)(-l)/"(c) 

= -2f" (c) - 3[/"(c)] 2 
= 25/(c) < 0, 
and the result follows from Theorem 1.5.4 (iii). 

(ii) This now follows from Theorem 1.5.4 (ii) □ 

Remark 1.5.8 The above proof shows how the Schwarzian derivative arises from 
differentiating g = f o f = f 2 . In the case where /'(c) = -1, it follows that 

g"(c)=0 and Sf(c) = i</"(c). 

Example 1.5.9 For the logistic map L^x) = fix{l - x) we have 

L'^x) = n - 2fix, Ll(x) = -2fji and L"(x) = 0. 
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When /a = 1, x = 0 is the only fixed point and Theorem 1.5.4 (i) shows that x = 0 is 
semi-stable (attracting on the right). However, we regard this as a stable fixed point 
for L M defined on the interval [0, 1], since points to the left of 0 are not in the domain 
of L M . 

When fi = 3, c = 2/3 is fixed and 1^(2/3) = -1 giving a non-hyperbolic fixed 
point. However, 5/(2/3) = 0 - f [6] 2 < 0 (negative Schwarzian derivative), so by 
Theorem 1.5.7 (i), x = 2/3 is asymptotically stable. 

1.0 

0.8 
0.6 
0.4 
0.2 
0.0 

0.0 0.2 0.4 0.6 0.8 1.0 



Exercises 1.5 



1. Show that f(x) = -2x 3 +2x 2 +x has two non-hyperbolic fixed points and determine 
their stability. 

2. For the family of quadratic maps Q c (x) = x 2 + c, x G R, use the Theorems of 
Section 1.5 to determine the stability of the fixed points for all possible values of c. 
Find any values of c so that Q c has a non-hyperbolic fixed point, and determine their 
stability. 
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3. Find the fixed points of the following maps and use the appropriate theorems to 
determine whether they are asymptotically stable, semi-stable or unstable: 



x z x 



(i) fix) = — + - , (ii) f(x) = arctan x, (iii) f(x) = x 3 + x 2 + x, 



-«*-«*+«. w/w = { 0 .8^,) lift 



4. Let fix) =ax 2 + bx + c,a^ 0, and p a fixed point of /. Prove the following: 

(i) If f'ip) = 1, then p is semistable. 

(ii) If f'(p) = -1, then p is asymptotically stable. 



5. If f(x) = cx + d , a,b,c,d£R is the linear fractional transformation, show that 
its Schwarzian derivative is Sf(x) = 0 for all x in its domain. 



6. If Sfix) is the Schwarzian derivative of f(x), a C 3 function and Fix) = ^"^ X \ 
show that Sfix) = F\x) - iFix)) 2 /2. 



7. If / : R - R is defined by /(,) = ( « * * 0 ^ find ^ ^.^ ^ 

of / and show that for x 0 they are non-hyperbolic. Show that x = 0 is not an 
isolated fixed point (i.e., there are other fixed points arbitrarily close to 0). Is x = 0 
a stable, attracting or repelling fixed point? 



8. ([37]) Let N f be the Newton function of a four times continuously differentiable 
function /. If /(a) = 0, show that N'j'{a) = 25/(q), where Sf is the Schwarzian 
derivative of /. 
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9. (a) Use the Intermediate Value Theorem to show that f(x) = cos(x) has a fixed 
point c in the interval [0, it/ 2}. We can show experimentally that this fixed point is 
approximately c = .739085 . . for example by iterating any x 0 e R. 

(b) Show that the basin of attraction of c is all of R (Hint: You may assume that 
x e [-1, 1] - why? Now use the Mean Value Theorem to show that \ f(x)-c\ < X\x-c\ 
for some 0 < A < 1). 

(c) Does f(x) have any eventual fixed points? 

(d) Can f(x) have any points p with f 2 (p) = p other than c? 




(Sj^/fc P^^waCcs 21^ 

Mr 
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Chapter 2, Bifurcations and the Logistic Family 

2.1 The Basin of Attraction 

In this chapter we examine the basin of attraction of the logistic maps, i.e., for a 
given p and fixed point x = c we look for the set of those x 6 [0, 1] which converge 
to c under iteration by L^x) = px(l - x). 

Definition 2.1.1 The basin of attraction B f (c) of a fixed point c of f(x) is the set 
of all x for which the sequence x n = f n (x) converges to c: 

B f (c) = {xeX: f n (x) -^c, as n^oo}. = 

The immediate basin of attraction of / is the largest interval containing c, contained 
in the basin of attraction of c. We first show that this is always an open interval when 
c is an attracting fixed point. 

Proposition 2.1.2 Let f : R -» R be a continuous function having an attracting 
fixed point p. The immediate basin of attraction of p is an open interval. 

vs/p 

Proof. Since p is an attracting fixed point there is e > 0 such that for all x <E I t = 
(p-e,p + e), f n (x)-+pasn^oo. Denote by J the largest interval containing p for 
which f n (x) -> p for x <E J, as n -> oo. 

Suppose that J = [a, b], a closed interval, then there exists r e Z + with f(a) 6 J £ . 
Now / r is also a continuous function, so points close to a will also get mapped into 
7 £ , leading to a contradiction. 

More specifically, there exists 5 > 0 such that if \x-a\ < 5, then \f r {x)- f(a)\ < tj, 
where V = min{|/ r (a) - (p- e)|, \(p+e) - f r (a)\}. Thus there are points x < a close 
to a for which f r (x) 6 I t , so / rn (x) -> p as n -> oo, a contradiction. We conclude 
that a <£ J (and similarly for 6), so J is an open interval. □ 

Example 2.1.3 1. If f(x) = rr 2 , the fixed points are c = 0 and c = 1, both hyperbolic, 
the first being attracting and the second repelling. Clearly 5/(0) = (-1,1) and 
B f( l ) = {-!> !}• We sometimes regard c = oo as an attractive fixed point of /, so 
that B f (oo) = [-oo, -1) U (1, oo]. 

2. If /:/—>/ is a continuous function on a closed interval / = [a, b] having an 
attracting fixed point p e (a, b), then we cannot exclude the possibility that the basin 
of attraction includes either a or 6, or both. On the other hand, if x = a is the 



i 
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attracting fixed point, the bsain of attraction may be a set of the form [a, c) for some 
c G (a, b}. Such a set can be regarded as being open as a subset of [a, 6]. 

2.2 The Logistic Family 

The logistic maps L^(x) = jix(\ — x) are functions of two real variables fi and x. 
We usually restrict x to the interval [0, 1] and consider /i G (0,4]. fi is a parameter 
which we allow to vary, but then we study the function for specific fixed values of 
/i. As the parameter \± is varied, we see a corresponding change in the nature of the 
function L^. This is what is called bifurcation. For example, for 0 < /x < 1, Ln has 
exactly one fixed point in [0, 1], c = 0, which is attracting. As \± increases beyond 1, a 
new fixed point c = 1 — is created in [0, 1], so now L M has two fixed points, c = 0 
is now repelling and c = 1 — is attracting (for 1 < /x < 3). At /i = 3 the nature 
of these fixed points changes again as we shall see. In this section we determine the 
basin of attraction of these fixed points as /i increases from 0 to 3. We see that the 
"dynamics" (long term behavior) of L M is quite uncomplicated for this range of values 
of //. 



The function L^(x) = fix(l - x), 0 < x < 1 has a maximum value of /x/4 when 
x = 1/2. Consequently, for 0 < \± < 4, maps the unit interval [0, 1] into itself. We 
shall consider later what happens when // > 4. We start by showing that the basin 
of attraction of for 0 < // < 1 is all of the domain of namely [0, 1]. We say 
that 0 is a global attractor in this case. 

Theorem 2.2.1 LetL^x) = iix{l-x), 0 < x < 1 be the logistic map. ForO < fi < 1, 
B Lfi (0) = [0, 1] and for 1< /x < 3, B L ^\ - = (0, 1). 

We split the proof into a number of different cases: 
Case 2.2.2 0 < /x < 1. 

e&wfk; &nsldar -tU ™f 3*L-lA~l *^ Jf \ 

A I '\&-3 >? /<K<4 
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0.0 0.2 0.4 0.6 0.8 1.0 

For 0 < n < 1, the only fixed point is 0. 

We have seen that for n £ (0, 1), has only the one fixed point x = 0 in [0, 1] 
(the other fixed point is 1 - < 0). For fi < 1 this fixed point is asymptotically 
stable, but in any case 

0 < fi < 1, 0 < 1 - x < 1 0 < fi(l - x ) < 1 
0 < L^x) = ^(1 - z) < x, xe (0, 1], 
and in a similar way L^x) < L^x) etc., so the sequence L^x) is decreasing, bounded 
below by 0, and hence must converge to the only fixed point, namely 0. It follows 
that the basin of attraction is #l m (0) = [0, 1] (L„(l) = 0). 

Case 2.2.3 1 < /x < 3. 

1.0 r 
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1 - < 1/2 for 1 < fi < 2. 1 - 1/// > 1/2 for 2 < /i < 3. 

We have seen that for /x > 1 the fixed point 0 is repelling, but a new fixed point 
c = 1 — has been "born" which is attracting (for 1 < /i < 3), so by Proposition 
2.1.2 there is a largest open interval / = (a, 6) containing the fixed point for which 
x e I implies L^x) — > c as n — > oo. If the basin of attraction of c is £? M (c), then 
0, 1 £ J5^(c) because L M (0) = 0 and 1^(1) = 0, so J5 M (c) ^ [0, 1]. Furthermore, clearly 
a, 6 £ S M (c). 

From the Intermediate Value Theorem, L^(a, b) is an interval which must be con- 
tained in (a,/;), for if x G (a, 6), L^(L M (x')) — > c as n -> oo. Let x n be a sequence 
in (a, 6) with lim n _>oo z n = a, then by the continuity of L M , linin-oo L^{x n ) = L M (a). 
Since x n G (a, 6) for every n G Z+, we have L^{x n ) G (a, 6) for all n G Z + . Since 
L^(a) ^ (a, 6), the only way this is possible is if L M (a) = a or L M (a) = 6, and similarly 
for 6. This is only possible if a and b are fixed points, are eventual fixed points or 
L M (a) = 6,1^(6) = a (we call {a, 6} a 2-cycle - see Section 2.3). We show that the 
latter case cannot happen, and this leads to the conclusion that a = 0 and 6=1, 
since there are no other fixed or eventual fixed points in [0, 1] that can satisfy these 
condition. Consequently, we must have £? M (1 - l//x) = (0, 1). 

Now we can check that 

Ll(x)-x = fi 2 x(l-(/j+l)x + 2fix 2 -fix 3 )-x = x(/i 2 x 2 -/i(/i+l)x+/i+l)(//x-/i+l), 

and a and b must satisfy this equation. We can disregard the linear factors as they 
give the two fixed points. The discriminant of the quadratic factor is // 2 (/i + l) 2 - 
4// 2 (//+l) = // 2 (/i+l)(//-3) < 0 for 1 < fi < 3, so there is no 2-cycle when 1 < /z < 3. 
When fi = 3, the discriminant is zero, the fixed point is c = 2/3 and the quadratic 
factor gives rise to no additional roots, so again there is no 2-cycle. □ 

Remark 2.2.4 Note that when /j, = 2, x = 1/2 is a super attracting fixed point (since 
14(1/2) = 0). As we saw above the basin of attraction will be (0, 1). 



Exercises 2.2 

1. If L fl (x) = fix(l - x) is the logistic map, show that x = 1/2 is the only turning 
point of L 2 (x) for 0 < // < 2, but when (i > 2, two new turning points are created. 



42 

Use this to show that for 2 < n < 3, the interval [1/fi, 1-1///] is mapped by I? onto 
the interval [1/2, 1 - l//i]. 



2.3 Periodic Points 

Points with finite orbits are of importance in the study of dynamical systems and 
their long-term behavior: 

Definition 2.3.1 Let / : X -> X be a function with c e X. 

(i) c is a periodic point of f(x) with period r 6 Z+ if / r (c) = c and / fc (c) ^ c for 
0 < fc < r (in particular, c is a fixed point of f r ). We call r the period of c and the 
set 0(c) = {c, /(c), / 2 (c), . . . , .f - x (c)} is an r-q/c/e. We write 

Pev r (f) = {xeX:f(x) = x}, 

so that Fix(/) C Per„(/), n = 1,2,..., since the points in Per„(/) may not be of 
period n, but of some lesser period. 

(ii) c is eventually periodic for / if there exists m e Z + such that f m (c) is a periodic 
point of / (we assume that c is not a periodic point). 

(iii) c is stable (respectively asymptotically stable, unstable etc.) if it is a stable fixed 
point of f r . 

The following criteria for stability now immediately follows from Theorem 1.3.3: 

Theorem 2.3.2 Suppose that c is a point of period r for f and that f'(x) is continuous 
at x = c. If a = f(c), i = 0, 1, . . . ,r - 1 then: 

(i) c is asymptotically stable if 

l/'W'/'(c)'/'(c 2 )-/'(c r _i)l<l. 



(ii) c is unstable if 

\f(c 0 )-f(c 1 )-f(c 2 )...f(c r . 1 )\>l. 



Proof. Let us look at the case where r = 3 as this is typical; 
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O(c) = {cj(c),f(c)} = {co,c u c 2 }. 

Then 

i if{x)) = ^ {f{f2{x))) = nf\x))(f\x))> = f(f(x))f'(f(x))f'(x) 



= /"(c 2 )/ / (c 1 )/ , (co), when x = c. 

The result now follows from 1.3.3. □ 

Example 2.3.3 If we look at the graph of the tent map T(x), we see it has two fixed 
points, c = 0 and c = 2/3. If we graph T 3 it has eight fixed points: These arise from 
two 3-cycles: {2/7,4/7,6/7} and {2/9,4/9,8/9}, together with the two fixed points, 
so that 

Per 3 (T) = {0, 2/3, 2/7, 4/7, 6/7, 2/9, 4/9, 8/9}. 

1.0- 
0.8 
0.6 
0.4 
0.2 



0.0 0.2 0.4 0.6 0.8 1.0 

The graph of T 2 shows that T 2 has four fixed points coming from a 2-cycle {2/5, 4/5} 
and the two fixed points. Since |T 7 (x)| = 2 and \{T 2 )\x)\ = \T'(x)\\T'{Tx)\ = 4, etc. 
(except at points of non-differentiability), all these periodic points will be unstable. 
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0.0 0.2 0.4 0.6 0.8 1.0 



Example 2.3.4 Let f(x) = 1/x, x^0,x^±l. Note that f 2 {x) = x and f(x) + x 
for all such x, giving rise to the 2-cycle {x, 1/x}. In this case 

\f'(x)f(l/x)\ = \-l/x 2 (-x 2 )\ = l, 

so the theorem is inconclusive. However we see that the periodic points are stable 
but neither attracting nor repelling. 




Example 2.3.5 Consider the quadratic function f(x) = x 2 -2. We have seen that to 
find the fixed points we solve f(x) = x, or x 2 - 2 = x, x 2 - x - 2 = (x - 2)(x + 1) = 0, 
so x = 2 or x = — 1. 

To find the period 2-points we solve f 2 (x) = x or f 2 (x) - x = 0. This is simplified 
when we realize that the fixed points must be solutions of this equation, so that 
(x -2)(x + 1) is a factor. We can then check that 



f(x) - x = x 4 - 4x 2 - x + 2 = (x - 2){x + l)(x 2 + x-l). 
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Solving the quadratic gives 

-1±a/5 



x = 



so tnat { , } is a 2-cycle. In general, finding periodic points of 

quadratics can be complicated. If f(x) is a quadratic, f n (x) - x is a polynomial of 
degree 2 n . 

To check the stability, we calculate 

- v/5)/2)/'((-l + VE)/2)\ = |(-1 - VE)(-1 + VE)\ = |1 - 5| = 4 > 1, 
giving an unstable 2-cycle. 

Remark 2.3.6 1. As before, periodic points can be stable but not attracting (as 
above with f(x) = 1/x at x ^ 1). They can also be attracting but not stable as in 
Example 1.4.7). 

2. Functions such as f(x) = sin a; can have no period 2 points or points of a higher 
period since this would contradict the basin of attraction of x = 0 being all of E. 
Similarly, the logistic map 0 < \x < 3 cannot have period n-points for n > 1. 



Exercises 2.3 

1. For each of the following functions, c = 0 lies on a- periodic cycle. Classify this 
cycle as attracting, repelling or neutral (non-hyperbolic). Say if it is super-attracting: 

W f( x ) = ^cosx, (ii) g( x ) = ~lx 3 - ^x 2 + 1. 



2. (a) Show that C^x) = //cos(x) has a super-attracting 3-cycle {0,A,tt/2} where 
ix . = A and A satisfies the equation Acos(A) = 7r/2. 

(b) Give similar conditions for S^x) = (Msm(x) to have (i) a super-attracting 2-cycle, 
(ii) a super-attracting 3-cycle. 

(c) Explain why for families of maps, say one member of a super-attracting n-cycle 
is a super- attracting fixed point (for a different value of//). 
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3. Find the fixed points and the period two points of following maps (if any) and 
determine the stability of the 2-cycle: 

(i) /(*) = ( fi )/W = fl --(M0), (iii) f(x) = ^py, (iv) f(x) = \x 



4. Let Q c (x) = x 2 + c. Show that for c < -3/4, Q c has a 2-cycle and find it explicitly. 
For what values of c is the 2-cycle attracting? 

5. Let f(x) = ax 2 + bx + c, a ^ 0. Let {x 0 , Xi} be a 2-cycle for f(x), 

(a) If f'(xo)f(x 1 ) = -1, prove that the 2-cycle is asymptotically stable. 

(b) If f(x 0 )f(x 1 ) = 1, prove that the 2-cycle is asymptotically stable. 

6. Let f(x) = ax 3 - bx + 1, a ^ 0. If {0, 1} is a 2-cycle for /(x), give conditions on a 
and 6 under which it is asymptotically stable. 

7. Let f(x) be a polynomial with f(c) = c. 

(i) If f'(c) = 1, show that (x - c) 2 is a factor of g(x) = f(x) - x. 

(ii) If /'(c) = 1, or /'(c) - -1, show that {x-cf is a factor of h(x) = f 2 (x)-x (i.e., 
if f(x) has a non-hyperbolic fixed point c, then c is a repeated root of f 2 (x) - x). 
Hint: Recall that a polynomial p(x) has (x - c) 2 as a factor if and only if both p(c) = 0 
and p'(c) = 0. 

(iii) Show in the case that f'(c) = -1 we actually have that (x - c) 3 is a factor of 
h(x) = f 2 (x)-x. 

(iv) Check that (iii) holds for the non-hyperbolic fixed point x = 2/3, of the logistic 
map L 3 (x) = 3x(l — x). 

(v) Check that (i), (ii) and (iii) hold for the (non-hyperbolic) fixed points of the 
polynomial f(x) = —2x 3 + 2x 2 + x. 
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8. Let / : R — > R be a function which is differentiable everywhere. If {x 0 , X\, . . . , } 
is an n-cycle for /, show that the derivative of f n is the same at each x u i = 
0,1,. ..,n - 1. 

9. Show that if x 0 is a peridic point of / with period n and f m (x 0 ) = x 0 m G Z+ , 
then m = kn for some A; G Z+. (Hint: Write m = qn + r for some r, 0 < r < n and 
show that r = 0) . 

10. Show that if = x and = x, and n is the highest common factor of 
p and then / n (x) = x. (Hint: Use the previous question and the fact that every 
common factor of p and q is a factor of n). 

11. (a) Let / : R -> R be an odd function (f{-x) = -f(x) for all x G R). Show that 
x = 0 is a fixed point and that the intersection of the graph of / with the line y = — x 
gives rise to points c of period 2 (when c ^ 0). 

(b) Use part (a) to find the 2-cycles of / : R -* R, f(x) = - 3x/2 and determine 
their stability. 

(c) ) Let / : R -> R be an e?;en function (f{-x) = /(x) for all x G R). Show that 
the intersection of the graph of / with the line y = -x gives rise to eventually fixed 
points c (when c^O). What are the eventually fixed points of f(x) = cos(x)? 

{xsin(l/x) if x ^ 0 
0 if x = 0 " Find the fixed points of f and determine 
their nature/stability. Show that x = 0 is a non-isolated fixed point of / (there are 
fixed points of / that are arbitrarily close to 0). Find the eventual fixed points of / 
(note that f(x) is an even function). 

(b) (Open ended question) Discuss the stability of the fixed points of f(x), their 
basins of attractions and the existence of period 2-points. Note that f(x) < x for all 

xeR. 
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/ \ tvt i. i \ f xcos(l/x) if x 7^ 0 _ — , _ „ 

(c) Now set = <^ 0 if x = 0 • nd the fixed points of 9> and show 
that in this case there are period 2-points (Hint: note that g{x) is an odd function - 
see the previous question). 



2.4 Periodic Points of the Logistic Map 

We try to find the 2-cycles of the logistic map L^x) = ftx(l - x), 0 < x < 1. To 
do this we solve the equation 

Ll(x)=x, 

or 

fix(l - x)[l - fix(l - x)] - x = 0, 

or 

-//V + 2/?x 3 - (// + f)x 2 + fx - x = 0, 

Clearly a; is a factor (since c = 0 is a fixed point of L^(x)) and similarly x - (1 - 1///) 
must be a factor so we get 

Lfa) -x = -x(/j,x l)0uV - n{n + l)x + y, + 1), 

giving a quadratic equation which has no roots if // < 3: 

fx 1 - + l)x + ix + 1 = 0. 

Solving using the quadratic formula gives 



_ (l + M)± x /(/i-3)Qu + l) 
2// 

This is real only for > 3 (called the "birth of period two"). Let us call these two 
roots ci and c 2 (dependent on //). 

This 2-cycle is asymptotically stable if 

H)'(ci)\ = |^( Ci )L;(c 2 )| < 1, 

or 

-K// 2 (l-2 Cl )(l-2c 2 ) < 1, 
-1 < (-1 - vV-2//-3))(-l + vV-2/*-3)) < 1, 
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-l<l-(/i 2 -2/i-3)<l, 
and this gives rise to the two inequalities 

/j 2 - 2/i - 3 > 0 and /i 2 - 2/i - 5 < 0, 

and solving gives 

3 < fi < 1 + y/6. 

This is the condition for asymptotic stability of the 2-cycle {c u c 2 }. 
For \± = 1 + x/6, it can be seen that 

^(ci)L;(c 2 ) = -1, and SL 2 ^) < 0, 

so Theorem 1.5.7 (i) shows that the 2-cycle is asymptotically stable. Also, the 2-cycle 
is unstable for \± > 1 + \/6. In summary: 

Theorem 2.4.1 For 3 < fi < 1 + \/6, £/ie logistic map L^x) = /ix(l - x) /ms an 
asymptotically stable 2-cycle. The 2-cycle is unstable for /i > 1 + \/6. 

The above shows that we have a bifurcation when /i = 3 where a 2-cycle is created 
and which was not previously present. There is another bifurcation at fi = 1 + \/6. 
This means that for 3 < \i < 1 + y/6, when we use graphical iteration of points close 
to ci and c 2 , they will approach the period 2 orbit, and not the fixed point (which is 
now unstable). In fact it can be shown that for this range of values of /i, the basin 
of attraction of the 2-cycle consists of all of (0, 1), (except for the fixed point 1 - l//i 
and eventual fixed points such as l/fi). When /i exceeds 1 + = 3.449499. . the 
period 2-points become unstable and this no longer happens, but we shall see that 
something different happens. We have another bifurcation when fi = 1 + \/6, with 
the birth of an attracting period 4-cycle. 

2.5 The Period Doubling Route to Chaos 

We summarize what we have determined so far and talk about what happens as fi 
increases from zero to around 3.57. 

For /x < bi := 1, c = 0 is the only fixed point and it is attracting for these values of 
fi. There is a bifurcation at bi = 1, where a non-zero fixed point c = 1- l//z is created. 
This fixed point is attracting for \i < 3 (and c = 0 is no longer attracting) and super 
attracting when = si = 2. The second bifurcation occurs when fi = b 2 = 3. The 
fixed point c = 1 - becomes unstable, and an attracting 2-cycle is created for 
3 < ii < 1 + V6 = b 3 . 
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2.5.1 A Super-Attracting Period 2-Cycle 

We saw that when ll = 2, c = 1/2 is a super-attracting fixed point for L„. We 
now look for a super-attracting period 2-cycle for when 3 < fi < 1 + ^6, as it 
illustrates an important general method that can be used for finding where period 
three is born: 

Suppose that {x h x 2 } is a 2-cycle for the logistic map L M , which is super-attracting, 
then 

Xi = /ix 2 (l - x 2 ), and x 2 = LiXi(l - Xi), 
so multiplying these equations together gives the equation 

(J?(l - xi)(l - x 2 ) = 1. 

In addition we must have 

{Ll)'{x x ) = L'rMLfa) = 0, 

so that 

At 2 (l - 2a?i)(l - 2x 2 ) = 0. 

Thus either x t = 1/2, or x 2 = 1/2, so suppose the former holds, then x 2 = fi/4. 
Substituting into the first equation gives 

/x 2 (l-/V4)(l-l/2) = l, 

or 

fi 3 - 4/i 2 + 8 = 0. 
(J. — 2 must is a factor of this cubic, so we have 

(//-2)(// 2 -2/i-4) =0, 
and this gives /x = 1 + y/§. We have shown: 

Proposition 2.5.2 When fx = s 2 = 1 + y/E, {1/2, ^-~^-} is a super- attracting 
2-cycle for . 

When fx exceeds b 3 = 1 + ^6, the 2-cycle ceases to be attracting and becomes 
repelling. In addition, a 4-cycle is created which is attracting until fi exceeds a value 
b 4 when it becomes repelling and an attracting 8-cycle is created. This type of period 
doubling continues so that when n exceeds b n , an attracting 2"- 1 -cycle is created until 
fi reaches b n+l . These cycles become super attracting at some s n (b n < s n < 6 n+1 ). 
This behavior continues with 2 n -cycles for all n € Z+ being created, until fx reaches 
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approximately 3.57. In other words, for b n < n < b n+1 , has a stable 2 n -cycle. It 
can be shown that = lim, woo b n = 3.570 approximately. 
Although b n+ i - b n — ► 0 as n — ► oo, it can be shown that 

hm -2 V = 5 = 4.6692016 

n-oc b n= i - 6,, 

(called Feigenbaums's number). Feigenbaum showed that you get the same constant 6 
in this way for any family of unimodal maps (f^x) is unimodalif f^Q) = 0, / M (l) = 0, 
is continuous on [0, 1] with a single critical point between 0 and 1). 

2.6 The Bifurcation Diagram. 

The behavior described above can be illustrated graphically using a bifurcation 
diagram. To create a bifurcation diagram we plot /i, 0 < // < 4 along the x-axis, 
and values of L^(x) along the j/-axis. The idea is to calculate (say) for each value 
of /i the first 500 iterates of some (arbitrarily chosen) point x 0 . We ignore the first 
450 iterates and plot the next 50. So for example, if 1 < // < 3, because the fixed 
point is attracting, the iterates will approach the fixed point 1 - 1///, so for n large, 
what we see plotted will be (very close to) the value 1 - For 3 < fi < 1 + 
the fixed point has become repelling, so this no longer shows up, but the 2-cyck has 
become attra cting , so we see plotted the 2 points of the 2-cycle. This continues with 
the 4-cycle, 8-cycle etc. This is called the period doubling route to chaos. 

We can create a bifurcation diagram for using Mathematica in the following 
way: 

First we define the logistic map for // = 1, and then a function of two variables: 
f [xj := x(l-x) 

h[x_,a_] := a*f [x] 



For a given a G [0,4] we pick x e (0, 1) randomly and iterate h{x,a) 100 times. 
Since this gives rise to a single number, we can define a function as the output: 

g[aj: = ( 

k[x_] :=h[x,a] ; 

A=NestList [k, Random [Real, {0,1}], 100]; 
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Return [A [[100]]] 

) 

Now we generate a list using the Table command for a € [1,4]: 
B=Table[g[j/1000] ,{j, 1000,4000}] ; 

and now we plot this list: 

ListPlot[B, PlotStyle->PointSize [0.001] ,Ticks->False] 




The bifurcation diagram for the logistic map. 



2.6.1 Where Does Period Three Occur for the Logistic Map 

If we look at the graph of for values of fj, close to 3.8, we see that is where a 
3-cycle is created: 




1 1 ' 1 1 ' ' ' 1 1 1 1 1 1 1 X 1 I . _ 

0.0 0.2 0.4 0.6 0.8 1.0 

The Graph of Lj(rr) for fi < 1 + y/S . 




0.0 0.2 0.4 0.6 0.8 1.0 



The Graph of L*(x) for fx = 1 + 
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00 0.2 0.4 0.6 0.8 1.0 

The Graph of Lj(x) for // > 1 + y/8. 

We aim to show that this happens when fi = 1 + y/8. For values of ji slightly 
smaller than this we only see the two fixed points. For /j = 1 + ^8 we see that a 
3-cycle is born, and for larger values of //, we see two 3-cycles. For a small range of 
values of fi, the 3-cycle is attracting. 

To show that the bifurcation occurs when // = 1 + y/8 we follow the argument of 
Feng [16] (see also [5], [17] and [42]): 

Period 3-points of occur where Lj(rc) = x, so we look at where 

Ll(x) -x = 0. 
To disregard the fixed points of L M we set 

y " { ' L,{x)-x 
Using Mathematica (for example) we can check that 

g,{x) = //x 6 - (jf + 3 fi 6 )x 5 + (// 4 + 4// + 3/*V - (/? + 3// + 5/i 5 + ^ 6 )x 3 
+(/x 2 + 3/? + 3/x 4 + 2/iV - (/x + 2/i 2 + 2// + //)* + 1+ ^ + tf. 
Set A = 7 + 2yu - yu 2 and let 

hn(z)=g (t (-z/fi), 

then 

M*) = * 6 + (3/i + l)^ 5 + (3/i 2 + 4/i + 1)^ 4 + (// + 5/z 2 + 3/i + 1)^ 3 
+(2// 3 + 3// 2 + 3// + l)z 2 + (/i 3 + 2/j, 2 + 2// + l)z + // 2 + /i + 1. 



55 

Then if 

M2) = { 2 3 + z i(lE±ll + 2(2 , + 3-|) + ^-|} 2 + ^ + ,)* (2 + „)., 

then we can check that k^z) = h^z) for all 2 (e.g., using Mathematica). 
Note that 

A > 0 for fi < 1 + V8, A = 0 for /x = 1 + v^, and A < 0 for n > 1 + Vs. 

This means that for fj, < 1 + V§, ^(2) > 0 for all 2 (we say ^(2) is positive definite), 
so cannot have any roots, i.e., g^x) = 0 has no solution, so L M cannot have any 
3-cycles. We summarize this as follows: 

Theorem 2.6.2 [16] (i) IfO < fx < 1 + y/8, then h M (z) is positive definite and the 
equation h^(z) = 0 does not have any real roots. Consequently, the logistic map L^x) 
does not have a 3-cycle. 

(ii) If H = 1 + VS, then h^(z) = 0 has three distinct roots, each of multiplicity two. 
These three roots constitute a 3-cycle for L^(x). 

(hi) Ifn > 1 + -y/8 (with // - (1 + V8) sufficiently small), the equation h^z) = 0 has 
six simple roots which give rise to two 3-cycles for L^{x). 

Proof, (i) If (j, < 1 + y/8, then since h^(z) is positive definite the result follows. 

(ii) If /j, = 1+ V8, then A = 0 and the equation becomes 

h^z) = (z 3 + (2 + 3y/2)z 2 + (5 + 4>/2> + 3 + v^)' . 

The resulting cubic can be solved using the cubic formula to give three real solutions, 
zi, z 2 , z 3 and these can be used to give the three solutions to g^x) = 0, corresponding 
to the 3-cycle of L^(x): 

2V7 (\ , 1 2torA 2 + 3V2 

Zk = — cos U arccos( - 27f + ~T ] ) 3—' k = °' 2 > 

(see the graph of g^(x) below). 

(iii) For A < 0 we can factor h^z) = h 1 (z)h 2 (z) using the difference of two squares, 
and then use the Intermediate Value Theorem on each of h x (z) and h 2 {z), to see that 
they each have three different roots corresponding to two 3-cycles, which can then be 
shown to be distinct (see [16] for details). □ 



g(z) 




g^(x) for fi < 1 + a/8. 

g(z) 



5 




g^x) for /i > 1 + y/8. 
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Example 2.6.3: A Super- Attracting 3-Cycle for the Logistic Map 

Recall that super-attracting periodic points occur where the derivative is zero, so 
for a super- attracting 3-cycle {ci, c 2 , c 3 }, we require 

(^)^i) = ^(ci)^(c 2 )L;(c 3 ) = 0, 

i.e., 

(l-2c 1 )(l-2c 2 )(l-2c 3 ) = 0, 
so we may assume c\ = 1/2. This means that x = 1/2 is a solution of the equation 
L*(x) = x, or Lj(l/2) - 1/2. But if fi satisfies the equation 2^(1/2) = 1/2, then 
it will also satisfy the equation involving the third iterate. Consequently we solve 
for fi: 0/4(1/2) = 1/2, where g^(x) = (L*(x) - x)/(L^(x) - x) as defined in the last 
section (this eliminates the root ji = 2 which gave rise to the superstable fixed point 
at x = 1/2). We obtain 

-?-(64 + 32/i + 16/x 2 - 24/i 3 - 4/i 4 + 6// - //) = 0. 
64 

Set 

p{a) =a 6 - 6a 5 + 4a 4 + 24a 3 - 16a 2 - 32a - 64, 
then Mathematica indicates that there is a single real root /x 0 larger than 1 + with 
exact value 



^ = 1 16 + 2^3 (ll + 2 • 2 2 / 3 (25 - 3^)^ + 2 . 2 2 / 3 (25 + 3^69)^ j . 

Mathematica is able to solve this equation exactly because it can be reduced to a 
cubic (and then the cubic formula may be used): 

Following Lee [26] , replace a by a + 1 and check (using Mathematica) that 

p(a + 1) = a 6 - 11a 4 + 35a 2 - 89 = b 3 - lib 2 + 356 - 89, 

a cubic in b = a 2 which can now be solved exactly for b, then for a, from which the 
original equation can be solved. It is seen that 

Ho = 3.8318740552 .... 

The other period 3-points may now be found since c\ = 1/2, c 2 = 1/^(1/2) = /xq/4 = 
0.95796. . ., and c 3 = £j 0 (l/2) = /xg/4(l - /x 0 /4) - 0.15248. . .. 

Example 2.6.4: The 3-Cycle when fi = 4 

When ji — 4, Z/4(x) = 4x(l — x). In this case we have 

g A (x) = 4096x 6 - 13312x 5 + 16640x 4 - 10048x 3 + 3024x 2 - 420x + 21, 
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and we can check using Mathematica (see also Lee [27]) 

g 4 (x/4) = x 6 - 13x 5 + 65x 4 - 157a: 3 + 189:r 2 - 105a; + 21 
= (x 3 - 7x 2 + Ux - 7)(x 3 - 6x 2 + 9x- 3), 

the product of two cubics. The solutions give a pair of 3-cycles, which we will show 
are given by the following theorem: 

Theorem 2.6.5 For the logistic map L 4 (x) = 4x(l - x) we have: 
(i) The 2-cycle is {sin 2 (7r/5), sin 2 (27r/5)} 
(it) The 3-cycles are 
{sin 2 (7r/7), sin 2 (27r/7), sin 2 (37r/7)} and {sin 2 (7r/9), sin 2 (27r/9), sin 2 (47r/9) }. 

Proof. Recall that the difference equation x n+1 = 4x n (l - x n ), n = 0, 1, 2, has 

. \ /; ' " 7***5 ^ 

the solution 

x n = sin 2 (2 n arcsin y^) . 

This was obtained by setting x n = sin 2 (6> n ) for some 9 n e (0, tt/2], (n = 1,2. . .) so 
that 

sin 2 (0 n+1 ) = 4sin 2 (0„)(l - sin 2 (0 n )) = 4sin 2 (0 n ) cos 2 (0 n ) = sin 2 (20„). 
We can use this to show that sin 2 (0 n+1 ) = sin 2 (4^ n _ 1 ), so in general 

x n = sm 2 (6 n )=sm 2 (2 n 8 0 ), where 9 0 = arcsin (y^). 
In particular we have 

6 1 = arcsin ( = I 29 » for 0 < °° < */4 
V V \ n-2e 0 for tt/4<0 o <tt/2 

Applying this to the situation where we have a 2-cycle {c 0 , Ci}, if a = sin 2 (^) we 
get 0 O is equal to 49 0 , n - 26 0 , 2n - 46 0 or tt - 46 0 . This gives 6 = 0 or 6 = tt/3 (giving 
rise to the two fixed points) or 0 O = tt/5, or 2tt/5 from which the result follows. A 
similar, but more complicated analysis gives the 3-cycles. □ 

Remark 2.6.6 Using the above ideas, if we want to find the period n-points of L 4 
we solve the equation L n A {x) = x. When x = sm 2 (6), this becomes 

sin 2 (<9) =sin 2 (2 n 6>). 
This gives rise to the two equations 

±9 = 2 n 6 + 2kn, or ± 9 = (2k + 1)tt - 2 n 0, for some k e Z. 
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This can be summarized as a single equation: 

/."TT 

±6 = 2 n 6 + kn^ 9 = — — , n = l,2,3...,fceZ 

2 n ± 1 

so that 

Per n (L 4 ) = {sin^^-) : 0 < k < 2"" 1 } U {sin 2 !^) : 0 < k < 2"" 1 }. 

It follows that L 4 has points of all possible periods. We shall see that the set of all 
periodic points constitutes a "dense" set in [0, 1] and each of these points is unstable. 



Exercises 2.6 

1. Recall the family of maps defined by S^x) = /isin(x) for x G [0, n] and /i G [0, tt]. 
Use the Mathematica commands below together with the ManipulatePlot command 
to estimate the values of /jl where periods two and three are born. Use 

h[a_,x_] := a*Sin[x] 
g[a_,x_] := h[a,h[a,x]] 
k[a_,x_] := h[a,h[a,h[a,x]]] 

Manipulate [Plot [{x,g[a,x]}, {x, 0, Pi}, PlotRange->{0,Pi}, 
AspectRatio->Automatic] , {a,0,Pi}] 

(Click on the + at the end of the Manipulate line to obtain the Animation 
Controls). 

2. Modify the bifurcation diagram of the logistic maps to give a bifurcation diagram 
for the family 5^, 0 < fi < tt. 

3. Do the same as in question 1 for the family C^{x) = ficos(x), x G [— 7r,7r] and 

fi€ [0,4 



4. Can you use the method of question 1 to estimate a value of fx for which 5 M has a 
superattracting 2-cycle, or 3-cycle? 



60 

(1 - x) 

5. Let g^x) = fix ) ' n > 0. 

{1+x) 

(a) Show that 9ll has a maximum at x = y/2- 1 and the maximum value is //(3-2 v / 2). 

(b) Deduce that g is a dynamical system on [0, 1] for 0 < n < 3+2\/2 (i e q (TO ll) C 
[0,1]). " u ' u ~ 

(c) Find the fixed points of g fl for n > 1. 

(d) Give conditions on // for the fixed points of ^ to be attracting. 

(e) Use Mathematica to graph gl and gj, and estimate when a period 2-point is 
created. 

(f) Use Mathematica to give a bifurcation diagram for g^ for 0 < // < 3 + 2^2. 



6] 



2.7 The Tent Family T M . 

We define another parameterized family in a piecewise linear manner: T M : [0, 1] — > 
[0, l],0</i<2, 

T M _ / 1* if 0 < x < 1/2 

\ /i(l-x) if 1/2<x<1. 

When ji = 2 we get the familiar tent map T(x) seen earlier. We now look at the 
parameter values 0 < ji < 2 (and later we shall see what happens when \± > 2, and 
where : M — » R). If 0 < /i < 1 we see that the only fixed point is c = 0. If // = 1, 
then all c G [0, 1/2] are fixed points and if 1 < fi < 2; then there are 2 fixed points. 
We look at each of these cases separately. 

1.0 

0.8 

0.6 
0.4 
0.2 



0.0 0.2 0.4 0.6 0.8 1.0 

Case 2.7.1 0 < /z < 1. We see that 0 is the only fixed point of and if 0 < x < 1/2, 
then 0 < T^(x) = fix < x, and if 1/2 < x < 1, then 

0 < T M (x) = //(l -x)<l-x<i<x, 

so in a similar manner to what we have seen previously, the sequence T™{x) is decreas- 
ing, bounded below by 0, so must converge to the fixed point 0. Thus Bt^(0) = [0, 1]. 

Case 2.7.2 /i — 1 Clearly if 0 < x < 1/2, then T\(x) = x, so is a fixed point, and if 
1/2 < x < 1, then x is an eventual fixed point since 

T?(x) = T l {l-x)=x. 

The fixed points are stable but not attracting. 
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Case 2.7.3 1 < /x < 2. There is a bifurcation at n = 1 with a second fixed point c 
created, 1/2 < c < 1. We solve T^x) = x to find this point: 



= - x) = so that c = 



1+fJL 



Since |7%(a;)| = > 1 in this case, the fixed point is repelling. 




Case 2.7.4 p = 2. In this case T 2 = T, the familiar tent map. The fixed points 0 and 
2/3 are repelling. This time the range of T is all of [0, 1]. T has the effect of mapping 
the interval [0, 1/2] onto all of [0, 1] and then folding the interval [1/2, 1] back over 
the interval [0, 1]. It is this stretching and folding that gives rise to the chaotic nature 
of T that we will see later, and is typical of many transformations of this type We 
will examine the periodic points of T in the next section. 



2.8 The 2-Cycles and 3-Cycles of the Tent Family 

At some stage for » > 1 a 2-cycle is created. It is interesting to use Mathematica 
to do a dynamic iteration of T, to see how this happens (use the ManipulatePlot 
command described in Exercises 2.6, but with 

h[x_,a_]:= a*Piecewise[{{x, x < 1/2}, {(1-x), 1/2 < x}}] 
in place of a*Sin [x] ). It can be checked that for » > 1, is given by the formula 

V 2 x if 0 < x < i 

Tl{x)=\= M-l*) if h<*<\ 

H{\-H + Iix) if 1 <X <1-J- 
A/ 2 (l~x) if \-±. <x <'{ 
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The 2-cycle is created when 7j(l/2) = 1/2, i.e., when 

nil - ti/2) = 1/2, or (//-1) 2 = 0, /z = l, 

so the 2-cycle appears when // > 1. For fi < 1, there are no period 2 points. The 
2-cycle {ci,c 2 } say, is unstable since 

\{Tl)\c 1 )\ = \Tl{c 1 )Tl{ci)\ = i?>\. 

For example, if we solve - fix) = x, we get c x = ^ and solving (j 2 (l-x) = x 
gives c 2 = as the period 2-points. Solving - fi + fix) = x gives c = ^ as 
the (non-zero) fixed point. 




00 02 04 °- 6 m ' o 7J for /i > 1 (but close to 1). 
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We now look for the smallest value of // for which there exists a 3-cycle. In a similar 
way to the situation for the 2-cycle, we look for where 2j(l/2) = 1/2 and we compute 
the value of the largest root of this equation. In general (see Heidel [19]) the smallest 
value of fi for which there exists a periodic orbit of period k is precisely the value of 
H for which 1/2 has period k. Using the formula for 2j(l/2) = //(l - ///2) above, 
then since - y,/2) < 1/2 for all fi, we get 

2j(l/2) = /x 2 (l - /x/2) = 1/2 when ^-2^ + 1 = 0, 

or (a* - 1)(^ 2 - /x - 1) = 0. Disregarding \l = 1 and solving the quadratic gives 
H = (1 + v5)/2 as the value of // where period 3 first occurs. 

In general it can be shown that for k > 3 odd, period k first occurs when is equal 
to the largest real root of the equation 

// - 2/i*" 1 + 2/- 3 - 2//~ 4 + . . . - 2/i + 1 

= - lX/- 1 - /- 2 - + /- 4 - A 5 + n k ~ e . • • + n - 1) = o, 

so has /x - 1 as a factor (see [19]). We will return to look at the tent family in Section 
6.4. 



Exercises 2.8 

1. Find /x e [0, 2] such that c = 1/2 is periodic of period 3 under the tent map T M (so 
that {1/2, /i/2, Ml - fi/2)} is a 3-cycle). 



2. For (j, > 1, show that T M has no attracting periodic points. 
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3. If L^x) = nx(l - x) is such that c = 1/2 is n-periodic for some n e Z + , prove 
that 1/2 is an attracting periodic point. Is it necessarily a super-attracting periodic 
point? 

4. Modify the bifurcation diagram of the logistic maps to give a bifurcation diagram 
for the tent family x € [0, 1] and fi € [0, 2]. 
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Chapter 3. Sharkovsky's Theorem 

3.1 Period 3 Implies Chaos. In 1975 in a paper entitled "Period three implies 
chaos" , Li and Yorke proved a remarkable theorem: 

Theorem 3.1.1 Let f : X — ► X be a continuous function defined on an interval 
X C R. If f(x) has a point of period three, then for any k = 1, 2, 3, . . ., there is a 
point having period k. 

This paper stirred considerable interest and shortly after it was pointed out that a 
Ukrainian mathematician by the name of Sharkovsky had in 1964 published a much 
more general theorem (in Russian) in a Ukrainian journal. His theorem was unknown 
in the west until the appearance of the Li- Yorke Theorem. To state his theorem 
we need to define a new ordering of the positive integers Z + . In the "Sharkovsky 
ordering", 3 is the largest number, followed by 5 then 7 (all of the odd integers), then 
2 • 3, 2 • 5, (2 times the odd integers), then 2 2 times the odd integers etc., finishing off 
with powers of 2 in descending order: 

3>5>7>--->2-3[>2-5t>--->2 2 -3>2 2 -5>---c>2"-3>2"-5c>---c>2 n >2"- 1 >---2 3 >2 2 >2>l. 

Sharkovsky's Theorem says that if a continuous map has a point of period k, then 
it has points of all periods less than k in the Sharkovsky ordering. The converse is 
also true in the sense that for each k e Z+ there is a map having points of period k, 
but no points of period larger than k in the Sharkovsky ordering. 

Theorem 3.1.2 (Sharkovsky's Theorem, 1964.) Let f : X X be a continuous 
map on an interval X (where X may be any bounded or unbounded subinterval ofR). 
If f has a point of period k, then it has points of period r for all r e Z+ with k>r. 

For example, this theorem tells us that if / has a 4-cycle, then it also has a 2-cycle 
and a 1-cycle (fixed point). If / has a 3-cycle, then / has all other possible cycles. If 
/ has a 6-cycle, then since 6 = 2-3,/ will have 2 • 5, 2 • 7, . . . 2 2 • 3, 2 2 • 5 . . .-cycles, etc. 

We omit the proof of this theorem, but shall prove the case where k = 3 (often 
known as the Li- Yorke Theorem because of their independent proof of this result in 
1975). James Yorke is a professor at UMD, College Park and Li was his graduate 
student. A few years ago, Yorke and Benoit Mandelbrot were awarded the Japan Prize 
(the Japanese equivalent of the Nobel Prize) for their work in Dynamical Systems and 
Fractals (Mandelbrot is regarded as the "father" of fractals, and we will review some 
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of his work later). In order to prove this result we shall need some preliminary lemmas. 
The first was proved in Chapter 1 (Theorem 1.2.9). 

Lemma 3.1.3 Let f : / R be a continuous map, I an interval with J = f(I) D I, 
then f(x) has a fixed point in I. 

Lemma 3.1.4 Let f : I -* R be a continuous map. If J C /(/) is a closed bounded 
interval, then there exists a closed bounded interval A C / with /(A) = J. 

Proof. Write J = [a, b] for some a, b € R, a < b. 
There are r, s € / with f(r) = a and f(s) = b. Set 

0 = y, where y G I, \y - r\ is a minimum, and f(y) = b, 

(P exists from continuity and /(/?) = 6: /3 is the point of / that is closest to r with 
the property that f{(5) = b). Similarly, set 

a = x, where x lies between /3 and r, |x - 0\ is a minimum, and /(x) = a 

(possibly equal to r), so as before /(a) = a. In this way, there is no c G (a,/?) 
where /(c) = a, or /(c) = 6. Then the closed interval K (either K = [a,0] or 
A' = [P,a]), has the property f(K) = [a, 6], since firstly, /(a) = a, /(/?) = 6, so by 
the Intermediate Value Theorem [a, b] C /(A"). On the other hand, if ic € /(A), then 
u' = /(*) for some z between a and /?. This must give f(z) e [a, b] since otherwise, 
another application of the Intermediate Value Theorem would contradict the choices 
of a and (3. q 

3.1.5 Proof of Sharkovsky's Theorem for k ; = 3 

We are assuming that / has a point of period 3, so there is a 3-cycle {a, b, c} with 
f( a ) = b, f(b) = c and /(c) = a. We assume that a < b < c (the other case with 
a < c < b may be treated similarly). 

We give the idea of the proof by showing why there must be points of period one, 
two and four. Let 

[a,b] = L 0 and [b,c] = L\. 

Observe that 

f(L 0 ) D U and f{L x ) DL 0 UL,. 



Case 1. / has a fixed point. 
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Since 

fiL^DLoULi D Li, 
Lemma 3.1.3 implies that / has a fixed point in L\. 

Case 2. / has a point of period 2. 

This time we use 

f{L\) 2 £o U Li D Lo, 
so by Lemma 3.1.4 there is a set B C Li such that f(B) = L 0 . We then have 

f 2 (B) = f(L 0 )DL 1 DB, 

so by Lemma 3.1.3, B contains a fixed point c of f 2 . c is a period 2-point of / (and 
not a fixed point of /) because 

/(c) € L 0 and c € L u so /(c) ^ c. 

Case 3. / has a point of period 4. 

The above two constructions do not illustrate the general method, but the following 
construction is easily generalized to any number of fixed points greater than 3. Our 
aim is to show that there is a subset B of L x which is mapped first by / into L u then 
into L x again, then onto L 0 and then onto L l5 so that f 4 (B) D B. Thus / 4 has a fixed 
point c in B, which cannot be a point of lesser period because /(c) 6 L\, / 2 (c) e L u 
f(c) € L 0 and / 4 (c) € B (so cannot have /(c) = c, / 2 (c) = c or / 3 (c) = c). 

It is useful to think of 5 copies of L 0 U L x with / mapping the first to the second 
etc. as shown: 



We find sets Bi, B 2 and B 3 as follows: 

f(L 0 )DL u f(Li) D L 0 ULi, 
so there exists B\ C L\ such that f(Bi) = L 0 . 
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There exists B 2 C L\ such that f(B 2 ) = B 1 and there exists B 3 C Li such that 
/(5 3 ) = 5 2 . Set 5 = B 3) then 

f(B 3 ) = f(B 2 ) = B 1 , and so / 3 (£) = L 0 , f 4 (B)DL 1 2B 3 . 

In other words / 4 (5) 5 £?, so there exists c G £?, a fixed point of / 4 , which is not a 
point of period 3 or less, so must be a point of period 4. □ 

In general if a function has points of period 4 the most we can deduce is that there 
are points of period 2 and fixed points. However, the following is true: 

Proposition 3.1.6 If f : I —> I is continuous on an interval I with 

f(a) = 6, f(b) = c, f{c) = d, f(d) - a, a < b < c < d, 
i/ien f(x) has a point of period 3, so also has points of all other periods. 
Proof. We may assume that 

f[a,b] = [b,c), f[b,c] = [c,d\, f[c,d} = [a,d\. 

In particular, there exists B\ C [c,d\ with f{B\) = [c,d\. There exists B 2 C [c,d] 
with /(ft) = M- 

Take AT X C ft with f{K x ) = ft, then 

/ 3 (^ a ) = / 2 (ft) = /[6,c] = [c,rf]D^ 1 , 

so / 3 has a fixed point in K\ which is not a fixed point of f(x). 



□ 

The above proofs can be summarized with the following type of result: 

Proposition 3.1.7 Let f : I —> I be a continuous map, and let h and I 2 be two closed 
subintervals of I with at most one point in common. If f(I\) D I 2 and f(I 2 ) D hUl 2 , 
then f has a 3-cycle. 

Proof. Exercise: Use the ideas of the proof of 3.1.5. □ 



3.2 Converse of Sharkovsky's Theorem 
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As we mentioned, for each m £ Z + in the Sharkovsky ordering of Z+ Sharkovsky 
showed that there is a continuous map / : / -+ / (/ an interval), such that f(x) has 
a point of period m, but no point of period k for k > m. The following were shown 
(as usual, / is either the real line or an interval): 

Theorem 3.2.1 For every k G Z + , there exists a continuous map f : I —> I that has 
a k-cycle, but has no cycles of period n for any n appearing before k in the Sharkovsky 
ordering. 

Theorem 3.2.2 There exists a continuous map f : I —> I that has a 2 n -cycle, for 
every n G Z + and has no other cycles of any other period. 

Strictly speaking, Sharkovsky's Theorem is really the combination of Theorem 
3.1.2, Theorem 3.2.1 and Theorem 3.2.2 (see [32]), and sometimes the latter two 
theorems are referred to as the converse of Sharkovsky's Theorem (see [15]). We look 
at some particular cases of this: 

Example 3.2.3 Define a function / : [1, 5] - [1, 5] as shown by the graph below (so 

/(I) = 3, /(2)=5, /(3)=4, /(4)=2, /(5) = 1, 

with f(x) piecewise linear between these points). Then / has a point of period 5, but 
no point of period 3. 




Proof. Clearly no point from the set {1,2,3,4,5} is of period 3, but this set is a 
5-cycle. Also, Theorem 1.2.7 tells us that f(x) has a fixed point c, so c is a fixed point 
for f. We shall show that f has no other fixed points. Suppose to the contrary that 
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P has another fixed point a. Now we can check that: 

/ 3 [1,2] = [2,5], / 3 [2,3] = [3,5], / 3 [4,5] = [1,4], 

so / 3 cannot have a fixed point in the intervals [1,2], [2, 3] or [4, 5], so a must lie in 
the interval [3,4]. In fact / 3 [3,4] = [1,5] D [3,4] and we show that / 3 cannot have 
another fixed point in [3, 4] . 

If a € [3, 4], then f(a) 6 [2, 4], so either f(a) e [2, 3] or f(a) € [3, 4]. If the former 
holds, then / 2 (q) e [4,5] and f(a) e [1,2] which is impossible as we have to have 
f 3 (a) = ae [3,4]. 

Thus we must have /(a) € [3,4], so that / 2 (q) g [2,4]. Again there are two 
possibilities: if / 2 (q) € [2,3] then / 3 (q) e [4,5], another contradiction, so that 
/ 2 (a)e[3,4]. 

We have shown that the orbit of a: {a, /(a), / 2 (a)} is contained in the interval 
[3, 4]. On the interval [3, 4] we can check that f(x) is given by the straight line formula 

f(x) = 10 - 2x, and /(10/3) = 10/3, 

so c = 10/3 is the unique fixed point of /. Also 

f(x) = -10 + 4*, f(x) = 30-8x, 

also with the unique fixed point x = 10/3. It follows that / cannot have any points 
of period 3. □ 

Remark 3.2.2 The above can be directly generalized to give a map having points of 
period 2n + 1, but no points of period 2n - 1, n = 2, 3, . . . . 



Exercises 3.2 

1. Use the ideas of Chapter 3 to show that if / : R — ► R is continuous and has a 
2-cycle {a, b}, then / has fixed point. 



2. Show that the map f(x) = (x - l/x)/2, x ^ 0, has no fixed points but it has 
period 2-points. Find the 2-cycle and by looking at the graph of f 3 (x), check to see 
whether it has a 3-cycle. Why does this not contradict Sharkovsky's Theorem? 
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3. A map / : [1,7] -> [1,7] is defined so that /(l) = 4, /(2) = 7, /(3) = 6, /(4) = 
5,/(5) = 3,/(6) = 2,/(7) = 1, and the corresponding points are joined so the map 
is piecewise linear. Show that / has a 7-cycle but no 5-cycle. 

4. If F A (x) = 1 - Xx 2 for x € R, show 

(i) F A has fixed points for A > -1/4. 

(ii) F\ has a 2-cycle for A > 3/4. 

(iii) The 2-cycle is attracting for 3/4 < A < 5/4. 

5. Show that an increasing function / : R -> R cannot have a 3-cycle. Can it have a 
2-cycle? Answer the same questions when / is decreasing. 
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Chapter 5. Devaney's Definition of Chaos 

In this chapter we will give Devaney's definition of chaos for one-dimensional maps 
and also for more general maps defined on metric spaces. One-dimensional maps are 
functions / : I — > R for some interval ICR. It turns out that maps whose periodic 
points form a dense set often have highly chaotic properties. 

5.1 The Doubling Map and the Angle Doubling Map 

Example 5.1.1 The Doubling Map. 

The doubling map D : [0, 1] — » [0, 1] is defined by 



D(x) = 2x(mod 1) 



2x ; 0<x<l/2, 
2x- 1 ; l/2<x$l. 




The doubling map D. 



It is useful to describe D in terms of the binary expansion of a real number in [0,1]. 
Let x G [0, 1] with binary expansion 



x = -a^as . . . where a* = 0 or 1. 



In other words, 



oo 




1=1 



Suppose that a\ = 0 in the binary expansion of x, then 



D(x) =2X = — + — + •.. = -€12(13 .... 
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On the other hand, if a\ = 1, then 

D(a?) = 2x-l = (a! + y + — + ...)-i = _ + — + ... = .a 2 a 3 .... 
We see that in general 

D(-aia 2 a 3 . . .) = -a 2 a 3 . . . and J D n (-aia 2 a 3 . . .) = -a n +ia n+2 a n + 3 

Consequently, if x = -aia 2 . . . a n aia 2 . . . a n a\ . . . has an expansion which repeats every 
n places, then D n {x) = x, so that x is periodic of period n. 

For example D 2 (-010101 . . .) = -010101 . . . , so is a point of period 2. We use this 
to show that the set of periodic points of D are dense in [0, 1]. This can also be 
used to count the number of periodic points of period n. Notice that since D'(x) > 1 
everywhere it is defined, all of the periodic orbits of D are unstable. 

Proposition 5.1.2 The periodic points of the doubling map are dense in [0, 1]. 

Proof. Let e > 0 and choose N so large that 1/2 N < e. If x G [0, 1] it suffices to 
show that there is a periodic point y for D that is within e oi x. Suppose that the 



binary expansion of x is 



then we set 



i=i Z 



y = -aia 2 . . . a^aia 2 . . . a^ai 



a point of period N. Then 



y\-o J— "Z- J 

5o 



1 



\x 



-y\ = \ E E *=w <e ' 



j=N+l 



j=N+l 



(where hi = 0, 1 or — 1). 

Example 5.1.2 The Angle Doubling Map. 



-J 



As before, we denote by C = {z = a + ib : a, b G M}, the set of all complex numbers. 
If z = a + ib G C, then its absolute value (or modulus ) is given by |z| = \/a 2 + b 2 . 
The conjugate of z is z = a — ib and we can check that zz = \z\ 2 . We can represent 
C using points in the (complex) plane {(a, b) : a, 6 G R}. The unit circle S 1 in the 
complex plane is the set 



S 1 = {z € C : |*| 
Points in S 11 may be represented as: 

z = e ld = cos 0 + i sin 6>, for some 6 € 
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Here 9 is the argument of z (written Arg(z)), and it is the angle subtended by the 
ray from (0,0) to (a, b) (when z = a + ib, a, 6 £ R), and the rea/ azzs. 

5 1 is a metric space if the distance between two points z,w £ S 1 is defined to be 
the shortest distance between th e two points, going around the circle. We define a 
map / : S l — » S l by [J(z) = z 2 ?j This map is called the angle doubling map because 
of the effect it has on 9 =Arg(z): f(e l °) = e 2ld . We see that the angle 9 is doubled. 
It is clear that there are a lot of similarities between the doubling map and the angle 
doubling map, and we shall examine this in the next section when we study the notion 
of conjugacy for dynamical systems. First we show that the periodic points of / are 
dense in S l . 

Consider the periodic points of f(z) = z 2 \ Solving z 2 = z gives z = 1 (we can 
disregard z = 0), f 2 {z) = z gives z 4 = z or z 3 = 1 and continuing in this way we see 
that the periodic points are certain nth roots of unity. 

Proposition 5.1.3 The periodic points of the angle doubling map f : S 1 — > S 1 are 

dense in S 1 . ^ ~ 




Proof. If f n (z) = z for some n G Z + , then z = z or z : 
V we want to find the (2 n — l)th roots of unity. This gives: 



y 2 n -l _ 



1. Write z = e ld , then 



■•■2- e 



e (2 n -l)ifl _ e 2kni 



for some fc G Z + , 



, giving the[2" - 1 jdistinct roots: z fe = e^/P"" 1 ), fc = 0, 1, 2, . . . , 2 n 



— n-roo tj" 

2, showing that 





Per n (/) = ( e 2 W(2^-i) : 0 < /c < 2 ri - 1}, n G Z + . 

These points are equally spaced around the circle, a distance 27r/(2 n — 1) apart, 
which can be made arbitrarily small by taking n large enough. It follows that the 
periodic points are dense in S l . □ 

5.2 Transitivity 

Sometimes, given / : X — > X (X a metric space), when we iterate x 0 G X, the 
orbit 0(xq) = {x 0 , /(xo), . . .}, spreads itself evenly over X, so that 0(x Q ) is a dense 
set in X. This leads to: 

Definition 5.2.1 / : X — > X is said to be f topological^ ) transitiy e_ if there exists 
x 0 G X such that O(x 0 ) is a dense subset of X. A trans itive point ioi f is a point 
x 0 which has a dense orbit under /. If / is transitive, then there is a dense set of 
transitive points, since each point in 0(xq) will be a transitive point. 



r % 



±4> 



J 




•2S 

in 



H.I 



4 




'I 
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Exa mple 5.2.2 The doubling map D : [0, 1] — > [0, 1] is transitive: To show this 
we explicitly construct a point xo G [0, 1] which has a dense orbit under D. xo is 
defined using its binary expansion in the following way: first write down all possible 
"1-blocks", i.e., 0 followed by 1. Then write down all possible "2-blocks", i.e., 00, 01, 
10, 11, then all possible '3-blocks", i.e., 000, 001, 010, 011, 100, 101, 110, 111, and 
then continue in this way with all possible "4-blocks" etc. (we could write them down 
in the order in which they appear as in the binary expansion of the integers). This 
gives: 

x 0 = -01 00011011 000001010 011 100101110111 00000001. .. , 
a point of [0, 1]. 

To show that O(x 0 ) is dense in [0, 1], let y G [0, 1] with binary expansion 

oo 

y = - y }y^; • = Yl §' Vi = 0 or *» 

and let 5 > 0. 

Choose N so large that < S. All possible finite strings of 0's and l's appear in 
the binary expansion of x 0 , so the string 2/12/22/3 • • ■ Vn must also appear in the binary 
expansion of Xq. 

It follows that for some r 6 Z + we have 

D r (x 0 ) = -2/12/22/3 • • • 2/yvfyv+i&iv+2 • • • , for some b N+ i, b N+2 , 

so that 



\D r (x 0 ) - 2/| = I • 2/12/2 • • • yNb N+ ib N+2 ... - -2/12/2 • • • VnVn+iVn+2 ■ ■ 



i =N+l 

This shows that any point y G [0, 1] is arbitrarily close to the orbit of x 0 under Z), so 
this orbit is dense in [0, 1] □ 

Remark 5.2.3 It is quite easy to see that for a transitive map / : X — » X on a 
metric space, given any non-empty open sets U and V in X there exists m G Z + such 
that 

Less easy to see is the converse of this statement, which holds for complete separable 
metric spaces (X is separable if it has a countable dense subset - see Chapter 7 for 
the notion of completeness). This is the Birkhoff Transitivity Theorem: 
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Theorem 5.2.4 A continuous map f of a complete separable metric space X is 
transitive if and only if for every pair U and V of non-empty open subsets of X there 
exists m 6 Z+ with U n f m (V) ^ 0. 

Using this result it is easy to show directly that the angle doubling map is transitive, 
however, the proof is beyond the scope of this text. 

5.3 The Definition of Chaos 

Devaney was the first to define the notion of chaos, saying that a function is chaotic 
if it has a dense set of periodic points, it is transitive and also has what is called 
sensitive dependence on initial conditions (known as the ''butterfly effect" in the 
popular literature). Subsequently it was shown that the first two requirements imply 
the third, so we define chaos as follows: 

Definition 5.3.1 Let / : X — > X be a map of the metric space X, then / is said to 

be chaotic if: r'S+rov*^ cW^ c ) 

(i) The set of periodic points of / is dense in X. *Q ( ^ — ' ^^'Lv/iV 

(ii) / is transitive. v >Vt ^j^** ' ' 

Examples 5.3.2 1. Homeomorphisms or diffeomorphisms of an interval ICR. 
cannot be chaotic as they ar e never t ransitive. The same type of considerations apply 
to functions such as sinx, cosx, arctanx and the logistic map for 0 < fi < 3. 

2. We have shown that the doubling map D : [0, 1] — ► [0, 1] has a dense set of periodic 
points and is transitive, so it is a chaotic map. 

3. The tent map T : [0, 1] — ♦ [0, 1] is chaotic. Recall that if x € [0, 1] has a binary 

,i rr f 'a 2 a 3 a 4 . . . if a x = 0 , , ' . r 
expansion x — • a\ a<i as — , then I x = < , , , . f _ where a- = 1 if 

a{ — 0 and = 0 if a* = 1. More generally we can see using an induction argument 
that ([47]): 



jm x _ | ' a n+l 0, n +2 • • • if Cl n = 0 

\ • a' n+1 a' n+2 «n+3 • • . if Q>n = 1- 



We can use this to write down the periodic points of T. For example, the fixed points 
are x = 0 and x = -1010... = 2/3, and the period 2-points are 



= -OllOOllO. 



= 2/5 and x 2 = 



•11001100. 



= 4/5. 
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The period 3-points are 

•010010010010... = 2/7, .100100100100... = 4/7, .110110110110... = 6/7 
and 

•001110 001 110... = 2/9, -011100 011 100... = 4/9, -111000111000... = 8/9. 

Notice that points of the form x = kjT e (0, 1), k G Z + , are almost fixed since they 
have a binary expansion of the form 

x = • a\ d2 as . . . , a n 0 0 . . . , 

where a n = 0 if k is even and a n = 1 if A: is odd. It follows that T n x = 0 if /c is even 
and T n x = 1 if k is odd. In particular 

'k - 1 fc 

2 n ' 2™ 



= [0,1]. 



The Intermediate Value Theorem now implies that there is a fixed point of T n in the 
interval [^-, Since such intervals can be made arbitrarily small and cover all of 
[0, 1], the set of periodic points must be dense in [0, 1]. 




0.0 0.2 0.4 0.6 0.8 1.0 



The fixed and period 2-points of the tent map. 

We use these ideas to show that the periodic points of T are numbers in [0, 1] of 
the form x = r/s where r is an even integer and s is an odd integer. 

Proof. If x 6 (0, 1) is a periodic point for T, then T n x = x for some n G Z + . There 
are two cases to consider: Suppose that x has a binary expansion 



x = • a\ a 2 as . . . a n a nH _i . . . , where a n = 0, 



Ill 

then 

T n x = • a n+ i a n + 2 . . . a 2n . . . , 
so we must have a x = a n+u a 2 = a n+2 , • . . a n = a 2n = 0, and 

x = • ai a 2 . . . a n ai a 2 . . . a n a\ . . . , where a n = 0, 
Rewriting this gives 

/ai a 2 a n _! \ 1 /a x a 2 a n _i \ 

*- ly + ^ + --- + ^J + FAy + ^ + --- + 2^j + --- 

= (2! + 2* + + ( 1 N 

V2 T ? T H l l-l/2» J 

_ aig^ +a 2 2 n - 2 + ... + a n -i2 _ r 
2" -I ~ s ' 

where r is even and s is odd. 

In the second case x = • a x a 2 a 3 . . . , a n a n+ i . . ., where a n = 1, so that 

T n x = • a n+1 a^ +2 . . . a! 2n . . . , 
and this gives a 2 = a 2 = < +2 , . . . , a n = a' 2n , . . ., and 

x = • ai a 2 . . . a n _i 1 a 2 . . . aj^ 0 a x — 
We can now argue as before, but using the first 2n terms of x. 

Conversely, suppose that x = r/s <E (0, 1) where r is an even integer and s is an odd 
integer. Since s and 2 are coprime, we can apply Euler's generalization of Fermat's 
Theorem to give 

2^ s) = l(mod s), or 2 P - 1 = ks for some p,fcgZ + , 
(1 < At < 2 P - 1), where 0 is Euler's function. Write the binary expansion of kr as 

kr = a x 2 p ~ l + a 2 2 p " 2 + . . . a p _ 2 2 2 + a p _ 2 2, a x G {0, 1}, 
which is even, so 

^ = («t2T» + a 2 2^ + . . . a p _ 2 2 2 + a p _ l2 ) (^-^ 

= (^I + ^ + I Q P-! \ ( 1 \ 

^ 2 T 2 2 2p-^' \1 - 1/2P / 

= • a! a 2 ... a p _i 0 ax a 2 . . . a p _! 0 . . . , 
which is a point of period p (or less). □ 



oo 



r 
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It is now clear that T is transitive since if we define xo G (0, 1) having a bi- 
nary expansion consisting of all 1-blocks, all 2-blocks, all 3-blocks etc., as before, 
except that we insert a single zero between every block then it follows that if x = 
• x\ X2 . . . x n x n +\ . . ., then T p x 0 = • X\ Xi . . . x n 6 n +i ... for some p > 0, i.e., every 
block will appear in the iterates of Xq. As before we see that T is transitive, so is 
chaotic. 

4. It is possible for a map to be transitive without being chaotic (although for 
continuous functions / on intervals in R, this is not possible: see [50]). For example, 
consider the irrational rotation R a : S 1 — > S 1 defined by R a (z) = a • z for some 
(fixed) a G S 1 . To say that R a is an irrational rotation means that a n ^ 1 for any 
n G Z + , i.e., a is not an nth root of unity for any n G Z + . It can be shown that every 
zo G S 1 has a dense orbit (a transformation with this property is said to be minimal). 
However, suppose that K£(z) = z, then a n z = z or a n = 1, a contradiction, so that 
i? a has no periodic points. R a is an example of an isometry: points always stay the 
same distance apart: 



\R a (z) — Ra(w)\ = \az — aw\ = \a\\z — w\ = \z 



w\ 



Note that if instead we have a n = 1 for some n G Z + , then R%{z) = a n z = z for all 
z G 5 1 , so that i?™ is just the identity map (every point of S l is of period n). 

5.4 Some Symbolic Dynamics and the shift map 

Recall that the set of all infinite sequences of O's and l's: 

E = {uj = (si,s 2j S3, • • •) • Si = 0 or 1}, 
is a metric space with metric defined by 

V! = (Si,$ 2 ,...)> ^2 = (*1,*2,...) S S. 

This metric has the following properties: 

Property 5.4.1 If u\ = (si,s 2 ,...)j ^2 = (*i,*2, • • •) G E, with ^ = £ i? i = 
1, 2, . . . , n, then d(u;i, u; 2 ) < l/2 n . C v>f^) 



Proof. 



OO | | OO | | OO 1 1 

k=l k=n+l k=n+l ▲ 



2 k 2 r 

k=n+l ▲ 



□ 



s 




2. -2.0-^) = z «v' 
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Property 5.4.2 If d(u u u 2 ) < l/2 n , then s t = U for i = l,2,...,n. 

Proof. We give a proof by contradiction: Suppose that Sj = t 3 for some 1 < j < n, 
then 



The shift map cr (sometimes called the Bernoulli sh ift) is an important function 
defined on E. 

Definition 5.4.3 The shift map a : E — > E is defined by 

so for example <j(l, 0, 1,0,...) = (0, 1, 0, 1, . . .) and a 2 (l, 0, 1,0,...) = (1, 0, 1,0,.. .), 
so that if = (1,0, 1,0,.. .) and u 2 = (0, 1,0, 1, . . .), then {ui } v 2 } is a 2-cycle for a. 
In this way, it is easy to write down all of the points of period n. Any sequence which is 
eventually constant is clearly an eventually fixed point of E, and any sequence which 
is eventually periodic (such as (1, 1, 1, 0, 1, 0, 1, 0, 1, . . .)), is an eventually periodic 
point. 

Proposition 5.4.4 The shift map a : E — ► E is continuous, onto, but not one-to-one. 



Proof. Clearly a is onto but not one-to-one. To show that a is continuous, let e > 0, 
then we want to find S > 0 such that if 



We shall see that it suffices to take 8 = l/2 n+l if n is chosen so large that l/2 n < e. 

In this case, if d(u>i,u) 2 ) < 8 = l/2 n+1 , then from Property 5.4.2, s t = U for 
i = 1, 2, . . . , n + L Clearly the first n terms of the sequences a(ui) and g(uj 2 ) are 
equal, so by Property 5.4.1, d(a(ui), a(v 2 )) < l/2 n < e, so that a is continuous. □ 

We can now prove: 

Theorem 5.4.5 The shift map a : E — E is chaotic. 

Q) ~~ — 1 

Proof. We first show that the periodic points are dense in E. Let/a; = s 2 , • • -)i e ^ 5 
then it suffices to show that there is a sequence of periodic points uo n G E with u n — > a? 
as n — > oo. Set 
<^i = . . .), a period 1-point for a, 




a contradiction. 



□ 



rf(^i ? ^2) < then d(a(cji), ct(cj 2 )) < e. 
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^2 = (si, 5 2 , 5i, 5 2 , . • .)j a period 2-point for <r, 
<^3 = «2 3 s 3 , si, $2, s 3 , . . .), a period 3-point for a, 
and continuing in this way so that 

= ( si,S2, . ■ . ,s n , si, . . .), a period n-point for a. 
Since a; and u; n agree in the first n coordinates, d(u,u n ) < l/2 n , so d(uj,uj n ) -> 0 
as n — > oo, or cj n — > u as n — » oo. r u\Zs* y — 

To^how that a is transitive, we explicitly construct a point u; 0 € E having a dense 
orbit under a. Set 

w 0 = ( 01 00 01 10 1 1 000 0 01 010 101 on.. ) 

""V""' V s v " f ' 

1-blocks all possible all possible 3 blocks 

2-blocks 

and co ntinuing in this way so that all possible n-blocks appear in u Q . To see that 
O(u 0 ) = E, letjo; = s 2 , s 3 , . . .) j e E be arbitrary. Let e > 0 and choose n so 
large that 1/2" < e, then since u> 0 consists of all possible n-blocks, the sequence 
(si, S2, . . . , s„) must appear somewhere in u> 0 , i.e., there exists A; > 0 with 

cT fc (u>o) = (si,s 2 ,...,s n ,...), 
so that a; and a (u; 0 ) 5®[51onJheJrst_n coordinates. It follows that 



This shows that the orbit of u Q comes arbitrarily close to any member of E, so it is 
dense in E. It follows that a is transitive and so it is chaotic. □ 

Remark 5.4.6 It can be shown that (topological) properties such as being totally 
disconnected, perfect etc. are preserved by homeomorphisms. 

The shift space E and the Cantor set C are homeomorphic metric spaces. In 
addition, the half open interval [0, 1) and the unit circle S 1 in the complex plane are 
homeomorphic. In particular, E and C will have identical topological properties, and 
so will S 1 and [0, 1). It follows that E and [0, 1] cannot be homeomorphic as C and 
[0, 1] are not homeomorphic (C is totally disconnected, but [0, 1] is not). 

Proof. E is given its usual metric, and C has the metric induced from being a subset 
of R, so that d(x, y) = \x - y\ for x,y eC. 

We define a map h : C -> E by h{- a x a 2 a 3 ...) = ( Sl , s 2) s 3 , ...), where a t = 0 or 2 
and Si = di/2. Clearly h is both one-to-one and onto. We show that it is continuous 
at each Xq 6 C. 



115 



Let e > 0 and choose n so large that l/2 n < e. Set S = l/3 n , then if \xq — x\ < 5, 
both xo and x must lie in the same component (sub-interval) of S n of length l/3 n . It 
follows that x 0 and x must have an identical ternary expansions in the first n places. 
Correspondingly, h(x 0 ) and h(x) must have the same first n coordinates. It follows 
that d(h(x 0 ), h(x)) < l/2 n < e, so h is continuous at x 0 . In a similar way we see that 
h~ l is continuous. 

To see that [0, 1) and S l are homeomorphic, define h : [0, 1) — ► S 1 by h(x) = e 27Tlx , 
then h is clearly one-to-one and onto. The map h wraps the interval [0, 1) around the 
circle, so that strictly speaking we are looking at [0, 1] with the end points identified, 
so in this way h becomes continuous. □ 



Exercises 5.4 



1. Let / : X — ► X be a transitive map of the metric space Show that if U 
and V are non-empty open sets in X, there exists m 6 Z + with U D f m {V) ^ 0. 

2. Let F : [0, 1) -> [0, 1) be the tripling map F(x) = 3x mod 1. Follow the proof for 
the doubling map (but use ternary expansions) to show that F is transitive and the 
period points are dense (find the periodic points), and hence show that F is chaotic. 

3. Let D : [0, 1) -> [0, 1) be the doubling map. Show that |Per n (D)| = 2 n . 

4. Use Proposition 4.4.7 to show that for a continuous increasing function / : [a, b] — > 
[a, 6], the periodic points cannot be dense in [a, b] (so a homeomorphism of [a, 6] 
cannot be chaotic). 

5.5 Sensitive Dependence on Initial Conditions 

We now show that the original definition of chaos due to Devaney [12] follows from 
the definition we have given. This result is due to Banks et al. [4]. 
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Definition 5.5.1 Let / : X -> X be defined on a metric space (X,d). Then / has 
sensitive dependence on initial conditioris^ti^re&dsts 5j> 0 such that for any x G X 
and any open interval U containing x and points other than x, there is a point y G U 
and n G Z + with 

d(rtx),r(y))>s. 



'aJ^ s ^J^ Tllis is the Precise definition of the idea that iterates of points close to each other 
may eventual ly be widely apart, so that a map has sensitive dependence on initial 



^T^Vy conditions if there exist points arbitrarily close to x which are eventually at least 
v/'V-/ distance 8 away from x. It is important to know whether we have sensitive depen- 
^/yi dence when doing computations as round-off errors may be magnified after numerous 
yjfy^* ji£ Orations. For example, suppose we iterate the doubling map, starting with x 0 = 1/3 
\^ Oj^ f and x 1 = .333. After 10 iterations we have D 10 (x 0 ) = 1/3 and £> 10 (xi) = .92, more 
than distance 1/2 apart. 

y/f Examples 5.5.2 1. The linear map / : R -> R, /(x) = ax, |a| > 1, has sensitive 
^ ' dependence since if x ^ y, 

l/ n W-r(2/)| = |a n x-a n y| = ^1^-2/1-^00 as n -> oo. 



However, clearly the dynamics of / is not complicated (f is not chaotic)." ^ ^^'^ 

2. The shift map a : E — > E has sensitive dependence on initial conditions since if 
o;i, <J2 G £ with cji 7^ u; 2 , then they must differ at some coordinates, say Si^U. Then 

= (s h Si+ u ...)> and a 2 " 1 ^) = (t h t i+u ...), 

so that 

OO I i 

d(a^M^^M) = J2 lSl+k - 1 ~ U+k ~ l1 = \ + other terms > 1. 

k=i 2 2 2 



3. The angle doubling map / : S 1 -> 5 1 , /(z) = z 2 has sensitive dependence since if 
we iterate z = e l6 ,w = e i(j) G 5 1 , their distance apart doubles after each iteration. 

4. The doubling map D : [0, 1] — > [0, 1] can be seen directly to have sensitive depen- 
dence on initial conditions. This also follows from the following theorem: 

Theorem 5.5.3 (Banks et al. [4]) Let f : X — > X be a chaotic transformation, then 
\ 'js f ^ as sens ^ ve dependence on initial conditions. 

We first prove a preliminary result: 
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Lemma 5.5.4 Let f : X — > X be a transformation which has at least two different 
periodic orbits. Then there exists e > 0 such that for any x G X there is a periodic 
point p satisfying 

d(xj k (p))>e, for all k G Z+ . 

Proof. Let a and b be two periodic points with different orbits. Then d(f k (a), f l (b)) > 
0 for all k and / (since we are dealing with finite sets). 

Choose 6 > 0 small enough that d(f h (a) ) f l (b)) > 2e for all k and L Then 

d{f\a\ x) + d(x, f l (b)) > d(f k (a), f l (b)) > 26 Vfc, I G Z + , 

by the triangle inequality. 

If x is within e of any of the points f l (b), then it must be at a greater distance than 
e from all of the points f k (a) and the result follows. □ 

Proof of Theorem 5.5.3 Let x G X and U be an open set in X containing x. 
The periodic points of / are dense in X, so there is a periodic point q of period r 

say: 

q eV = UnB 6 (x), 

Let p be a periodic point of period n for /, whose orbit is a distance greater than 
45 form x, and write 

Wi = £,(/*(?)). 

Now 

/*(p)€Wi,Vt^p6/-*(Wi).Vt, 

so the open set 

W = rWi) n r 2 (w 2 ) n • • • n r n (W n ) ^ 0. 

Since / is transitive, there is a point z G V with / fc (z) G W for some k G Z + . 
Let j be the smallest integer with k < nj, or 

1 <nj — k <n. 

Then 

/^(z)=/^-*(/*(«))e/^-*(W). 

But 

r*- k {w) = r j - k (r 1 (w 1 )nr 2 (w 2 )n- ■ -nr n (w n )) c r- k {r {nj - k) w n] . k = w nj . k ) 
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so that d(f n i(z)J n i- k (p)) < 5. Now f n i(q) = q, so 

d(F j (q),r i (*)) = d(q,r i (z)) 
> d{x, f nj - k (p)) - d(f nj - k (p), f nj (z)) - d(q, x), 
by the triangle inequality since 

d( q , x) < d(x, r>- fc (p)) - d(r- k ( P ), r{ Z )) + 

So 

«!(r j (9),/ nJ W)> 45-5-5 = 25. 
This inequality implies that either 

d(f nj (x),r>(z))>5, 

or 

for if f nj (x) were within distance < 5 from both of these points, they would have to 
be within < 25 from each other, contradicting the top inequality above. So one of 
the two, z or q will serve as the y in the theorem with m = nj. □ 



f * i'f »a tU*. • 
- /a oj^Vv 




Chapter 6. Conjugacy of Dynamical Systems 
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Two metric spaces X and Y are the "same" (homeomorphic) if there is a homeo- 
morphism from one space to the other. In this chapter we study the question of when 
two dynamical systems are the same. Given maps / : X — ► X and g : Y — > Y, we 
require them to have the same type of dynamical behavior, e.g., if one is chaotic, then 
so is the other, there is a one-to-one correspondence between their periodic points etc. 
One obvious requirement is that the underlying metric spaces should be homeomor- 
phic. We have seen a lot of similarities between the logistic map L 4 (x) = 4x(l - x) 
and the tent map T(x) and this will be examined in this chapter together with other 
examples such as the shift map and circle maps. 

6.1 Conjugate Maps^L^V^ 1 

This "sameness" is given by the idea of conjugacy, a notion borrowed from group 
theory, where two members a and b of a group G are conjugate if there exists g G 
G with ag = gb. One of the central problems of one-dimensional dynamics and 
dynamical systems in general is to be able to tell whether or not two dynamical 
systems are conjugate. We will see that if one map has a 3-cycle and another map 
has no 3-cycle (for example), then the maps cannot be conjugate, or if one map has 2 
fixed points and the other has 3 fixed points, then they are not conjugate. These are 
examples of conjugacy invariants, which give criteria for maps to be non-conjugate. 
A generally harder question is establishing conjugacy between maps. 

Definition 6.1.1 ^1 Let / : X — > X and g : Y —*Y be maps of metric spaces. Then 
/ and g are said to be conjugate if there is a homeomorphism h : X — > Y such that 

The map h is a co njugacy between / and g. Obviously conjugacy is an equivalence 
relation. { ^JS) Vepfo</Vt 

^ If in the above definition we only require the map h : X — > Y to be continuous, 
then we say that g is a factor of /. If in addition, h is an onto map, then we say that 
g is a quasi-factor of /. . * -v") » ^QJ^Ji/' 

Exercises 6.1 (J)y \p 4— trw*W&*. rek***- 

C 

1. Prove that if / : X — > X and g : Y — » Y are conjugate maps of metric spaces, then 
/ is one-to-one if and only if g is one-to-one, and / is onto if and only if g is onto. 

h 
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2. Prove that if / and g are conjugate via h and / has a local maximum at x 0 , then 
g has a local maximum or minimum at h(xo). 

3. Suppose that h : [0, 1] -» [0, 1] is a conjugacy between /, g : [0, 1] -> [0, 1] where 
/(0) = f(l) = 0 and #(0) = g(l) = 0. Show that h is increasing on [0, 1]. Deduce 
that h maps the zero's of / to the zero's of g. 

4. ([6]) The function T n {x) = cos(n arccos(x)) is the nth Chebyshev polynomial. Show 
that T n is conjuagte to the map A n : [0, 1] -> [0, 1], the piecewise linear continuous 
map defined by joining the points (0, 0), (1/n, 1), (2/n, 0), (3/n, 1), . . . , ending with 
(1, 1) if n is odd, or (1,0) if n is even. Use the conjugacy map h : [0, 1] -> [0, 1], 
h(x) = cos(pix). (In [6], there is a generalization of this to maps T\, where A > 1 is 
a real number) . 

6.2 Properties of Conjugate Maps 

It is often easier to show that certain maps are chaotic indirectly by showing that 
they are conjugate to chaotic maps and using the following result: 

Proposition 6.2.1 // / : X X and g : Y -> Y are maps conjugate via a conjugacy 
h : X — > Y: h o / = g o h, then 

1. h o f n = g n o h for all n G Z+ , (so f n and g n are also conjugate). 

2. // c is a point of period m for f, then h(c) is a point of period m for g. c is 
attracting if and only if h(c) is attracting. 

3. / is transitive if and only if g is transitive. 

4. / has a dense set of periodic points if and only if g has a dense set of periodic 
points. 

5. f is chaotic if and only if g is chaotic. 

Proof. 1. hof 2 = hofof = gohof = gogoh = g 2 oh, and in the same way 
h o p = g 3 o h, and continuing inductively the result follows. 
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2. Suppose that /'(c) 7^ c for 0 < z < m and / m (c) = c, then /i o / 2 (c) ^ h(c) for 
0 < i < m since /i is one-to-one, and so o h(c) ^ h(c) for 0 < i < m. In addition, 
h o / m (c) = g m o /i(c), or /i(c) = g m (h(c)), so /i(c) is a period m point for g. 

We shall only show that if p is an attracting fixed point of / (so that there is an 
open ball B e (p) such that if x 6 B e (p) then f n (x) — > p as n — > cx>), then is an 
attracting fixed point of g. 

Let V — h(B e (p)), then since h is a homeomorphism, 1/ is open in Y and contains 
Let y G V, then h~ l (y) 6 B e (p), so that f n {h~ 1 (y)) — > p as n — > oc. 

Since /i is continuous, h(f n (h~ l (y))) — > h(p) as n — > oc, i.e., 

g n (y) = hof n oh~ l (y)->h(p) i as n -> oc, 
so /i(p) is attracting. 

3. Suppose that 0(2) = {z, /(z), / 2 (z), . . .} is dense in X and let V C F be a non- 
empty open set. Then since h is a homeomorphism, h~ l (V) is open in X, so there 
exists k e Z + with f k (z) e h~l{V). 

It follows that h(f k (z)) = g k (h(z)) G V, so that 

0(/z(z)) = {/i(z),^(z)),^(^)),--.} 
is dense in Y, i.e., g is transitive. Similarly, if g is transitive, then / is transitive. 

4. Suppose that / has a dense set of periodic points and let V C Y be non-empty 
and open. Then /i -1 (V0 is open in X, so contains periodic points of /. As in (3), 
we see that V contains periodic points of g. Similarly if g has a dense set of periodic 

points, so does /. 
r> v * ) * ■ 

s Sy^ ■ 5. This now follows from (3) and (4). 1 r □ 

Example 6.2.2 We remark that sensitive dependence on initial conditions is not a 



A one to have sensitive dependence, but the other not: Consider T : (0, oc) — > (0, oo), 
7 zJlp T(x) = 2x and S : R -> R defined by 5(x) = x + In 2. If if : (0, oc) -> R is defined 



y/ conjugacy invariant. It is possible for two maps on metric spaces to be conjugate, 



\ jj^ V by ! H{x) — In xj then if is a homeomorphism and we can check that H o T = S o if 

so T and 5 are conjugate, T has sensitive dependence, but S does not. 
( / It can be shown however, that if T : X — > X is a map on a compact metric space 
^ (for example X = [0, 1]) having sensitive dependence, then any map conjugate to 
a? ' T also has sensitive dependence. — r~ x 
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It can also be shown that the property of having negative Schwarzian derivative is 
not a conjugacy invariant. 

Example 6.2.3 1. The logistic map L 4 : [0. 1] -> [0, 1]. L 4 (x) = -l.r( 1 -,r) is conjugate 
to the tent map T : [0, 1] - [0, 1], T(x) = j %) j* J J ^ VJ . 

Proof. Define /i : [0, 1] -> [0, 1] by = sin 2 (7nr/2). We can see that h is one-to- 
one, onto and both h and h~ l are continuous, so it is a homeomorphism (it is not a 
diffeomorphism as h'(l) = 0). Also 



and 



U o h(x) = U (sin 2 ( ^)) = 4sin 2 (^) (l - sin 2 (^)) = sin 2 ^*), 

h o T(x) - h(Tx) - ( h ^ if 0 < x < 1/2 2 

1 j \ h(2-2x) if 1/2<*.<1 = sm M> 



so L 4 o ft = h o T and L 4 and T are conjugate. ©tiw L9z2 5f*\©CbS^ □ 
2. The doubling map D : [0, 1] -> [0, 1] is a quksf-factor of the shift map a : E -> E 

c r 

Proof. D is a factor of a since ■ c «*^ * ^ ©*Ao 

hoa = Doh, 

where A : E — > [0, 1] is defined by 

OO 

h(x 1 ,x 2l x 3 ,...) = . Xl x 2 x s . .. = Y^£i 

^ 2* 

is a continuous function (note that /i is not a homeomorphism since it is not one-to- 
one: for example h(l, 0, 0, . . .) = 1/2 = h(0, 1, 1, 1, . . .), but (1, 0, 0, . . .) ? (0, 1, 1, 1, . . .) 
in E). In addition /i(E) = [0, 1], so that h is onto and D is a quasi-factor of a □ 

y^. The logistic map L 4 is a quasi-factor of the angle doubling map f • S 1 -» S 1 

fU) = z 2 . 

Proof. Define h : S 1 - [0;l] by h(e ix ) = sin 2 x, then 

L 4 o /i(e-) = L 4 (sin 2 x) = 4sin 2 (l - sin 2 ;c) = sin 2 (2x), 

and 

hof(e ix ) = h(e 2ix ) =sm 2 (2x). 

h is clearly onto and continuous, but it is not one-to-one: h(e ix ) = h(e~ ix ), so L 4 is a 
quasi-factor of /, but h is not a conjugacy. □ 
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We can now show that many of the above maps are chaotic. In order to do this 
we need to weaken the conditions of Proposition 6.2.1. If we drop the requirement 
that h is a homeomorphism, but just require it to be continuous and onto, then we 
can show that if / is chaotic then so is g. In other words, if g is a quasi-f actor of / 
where / is chaotic, then g is also chaotic. This result will be useful in showing that 
a number of well known examples are chaotic. \A*i\-i<*jcA+r- 

Proposition 6.2.4 Leth : X -> Y be surjective (onto) and continuous. If f : X — ► X 
and g : Y -> Y satisfy h o f = g o h and f is chaotic, then g is chaotic. 

Before proving this proposition we need a lemma concerning continuous functions 
on metric spaces: 

Lemma 6.2.5 Let h : X — > Y be a continuous function of metric spaces and A a 
subset of X, then h(A) C h(A). 

h(x). We can find a sequence 



yeh(A). 

□ 

Proof of Proposition 6.2.4 Use Per(/) and Per(g) to denote the periodic points of 
/ and g respectively, then we saw earlier that h(Per(/)) C Per(p). Since / is chaotic, 
Per(/) = X, and since h is onto, h(X) = Y. Then using the lemma we have 

Y = h(X) = h(¥^jj) C h(Pev(f)) C P^), 

so that Per(g) = Y. In other words, the periodic points of g are dense in Y. 

/ is transitive so there exists x 0 G X with O f (x 0 ) = X (where we use the subscript 
to distinguish the orbits with respect to / and g) . Now 

h (O/(x 0 )) = h{f n (x 0 ) :neZ+} = {h o f n M : n 6 Z + } 
= {/o^ 0 ):n6Z + } = O,(/i(x 0 )) 5 

so that 

Y = h(X) = h(Ojfa)) c MoTM) = O p (M*o)), 
so that /i(x 0 ) is a transitive point for g. □ 



Proof. Let y G /i(i4), then there exists x e A with y 
x n e A with linin—oo x n = x. 

Then /i(x n ) G and since /i is continuous 

lim h(x n ) = h(x) = y so that 
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It is easily seen that Proposition 6^4 remains true if we replace the requirement 
that h be onto by requiring that h(X) be dense in Y. 

Theorem 6.2.6 The tent map T : [0, 1] -» [0, 1], the logistic map L 4 (x) = 4x(l - x), 
the angle doubling map f : S 1 -> 5 1 , /(z) = z 2 , and tfie Doubling map D : [0, 1] -* 
[0, 1] are all chaotic. 

Proof. Suppose that we can show that the angle-doubling map / is a quasi-factor of 
the doubling map D, then we have: 

The tent map T is conjugate to the logistic map L 4 , which is a quasi-factor of f(z), 
which is quasi-factor of D, which is a quasi-factor of a, the shift map. It has been 
shown that the shift map is chaotic. The result now follows. 

It therefore suffices to show that for / : S 1 -> S 1 , f(z) = z 2 and D : [0, 1] -> [0, 1], 
D(x) = 2x (mod 1), / is a quasi-factor of D. Define h : [0, 1] -> S 1 by h{x) = e 2nix , 
then h is continuous, onto and one-to-one everywhere except that h(0) = 1 = h(l). 
Now 

r 

/ o h(x) = f(e 2nix ) = e 4 ™, 

and 

h o D(x) = h(2x (mod 1)) = e^ 2x ( modl » = e 4 ™, 
or h o D = f o h, and / is a quasi-factor of D. □ 

Exercises 6.2 

1. If p(z) = z 3 on S 1 , show that ^ is the angle tripling map. Find the period 
points and show that they are dense in S 1 . Show that g is conjugate to the map 
F : [0, 1) -> [0, 1), F(x) = Sx mod 1 of Exercises 5.4. 



1. If T 3 : C -»• C is the tent map with n = 3, but restricted to the Cantor set. Show 
that T 3 is conjugate to the shift map a : E -> E. 



3. Prove that the doubling map D : [0, 1) -> [0, 1), D(x) = 2x (mod 1), and the angle 
doubling map / : S 1 -> S\ f(z) = z 2 , are conjugate. 
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4. Is the shift map a : E — > £ conjugate to the doubling map D? (It can be shown 
that the shift map a is conjugate to T 3 restricted to the Cantor set). 



5. Let U : [-1. 1] -> [-1, 1] be defined by £/(x) - 1 - 2x 2 and T 2 : [0, 1] -> [0, 1] be 



Prove that # : [0, 1] — > [-1, 1], H(x) = - cos(nx) defines a conjugacy between these 
maps. 

6. Prove that the map F : [-1,1] -> [—1,1], F(x) = 4x 3 - 3x is conjugate to 
T : [0, 1] -> [0, 1], T(x) = 3x mod(l), via the conjugacy /i : [0, 1] -> [-1, 1], - 
cos(7rx). (Hint: cos(3x) = 4cos 3 (x) - 3cos(x)). Deduce that F is chaotic. This 
question is related to Exercise 6.1, # 4 concerning the Chebyshev polynomials. 



6.3 Linear Conjugacy 

It is sometimes the case that the conjugacy between two real (or complex) functions 
is given by a map with a straight line graph (an affine map). This is called a linear 
conjugacy, and is stronger than the usual notion of conjugacy. 

Definition 6.3.1 For functions /:/—>/ and g : J — > J defined on subintervals of 
R, we say that / and g are linearly conjugate and that h is a linear conjugacy if h 
maps / onto J where h(x) = ax + b for some a, b G R, a ^ 0 and h o f = g o h. 

The following example gives a criterion for two quadratic functions to be linearly 
conjugate: 

Example 6.3.2 Let 



the tent map 




F(x) = ax 2 + bx + c and G(x) = rx 2 + sx + 1, 



where a^0 and r ^0. If 
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then F and G are linearly conjugate via the linear conjugacy 

(x) = -x + 

Proof. 

7 771 / \ 7 / 2 1 l \ a(ax 2 + bx + c) b-s 
h o F\x) = h(ax z + bx + c) = — - ^ 

r 2r 

0 a& 2ac 
= —x 2 + —x + 
r r 

G o h(x) = G(-x +~^r)= r(-x + ~^r) 2 + s(-x + ^r) + t 



and 




fa 2 2 a 
= r {-2 x2 + 2" 



,( 6 _ 5 ) (6-«)= 



sa bs — s 2 



+ —x + 

r 



and we see that these are equal if 



a 2 2 , (6 - s) 2 + 2fe - 2s 2 + Art 

= — x H x + 



4r 



c = 



o 2 - s 2 + 2s - 2b + Art 
4a 



+ < 



□ 



For example, if F is defined on the interval [0, 1], then 

ft(0)~ and Ml) = 2a + 6 " S 



2r 



so if a/r > 0, then F is conjugate to G on the interval 



2r ' 
6 — s 2a + b - s 



2r ' 



2r 



Example 6.3.3 1. If L,{x) = fix(l - x), and Q c (x) = x 2 + c and c = then 
Lf, on the interval [0,1] is linearly conjugate to Q c on the interval [-/i/2, fi/2}. If 
= 4, this shows that L 4 (x) = 4x(l - x) on [0, 1] is conjugate to Q c (x) = x 2 + c on 
the interval [-2, 2] when c = -2. In particular, Q_ 2 on [-2, 2] is chaotic. 

Proof. We apply Example 6.3.2 with 

a = -fx, b = fi, c = 0, r = 1, s = 0, 7j = c. 

In this case /i(0) = fx/2 and = -fi/2 and we can check that the conditions of 

2u — u 2 

the example hold when c = ^ . □ 
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2. On the other hand, if fi = 2, we see that L 2 (x) = 2x(l — x) on [0, 1] is conjugate 
to Qq(x) = x 2 on [—1,1]. Recall in Exercises 1.1 the difference equation x n+ i = 
2x n (l - x n ) transforms to y n +i = y 2 n on setting x n = (1 — y n )/2. This is just the fact 
that L 2 and Qo are conjugate via h(x) = —2x + 1. 

3. We can check that the logistic map L 4 is conjugate to F : [—1,1] — > [—1,1], 
F(x) = 2x 2 - 1. 

Exercises 6.3 

1. Show that every quadratic polynomial p(x) = ax 2 + bx + d is conjugate to a unique 
polynomial of the form f c (x) = x 2 + c. 

2. Prove that the logistic map L 4 is conjugate to F : [—1, 1] — ► [—1, 1], F(x) = 2x 2 — l. 

3. Show that the logistic map L^(x) = /ix(l — x), x G [0, 1] is conjugate to the logistic 
type map F^(x) = (2 — fj)x(l — x) (/i ^ 2), via the linear conjugacy (which is defined 
on the interval with end points |^ and 2Z~) : 

7 / x M 1 — M 

/iW = — ^— x + - — 
v } 2-a 2-a 



4. Let f c (x) — x 2 -\- c. Show that if c^>^JL/4, there is a unique \i > 0 such that f c is 
conjugate to L fl (x) — /ix(l — x) via a map of the form h(x) = ax + b. 

5. (a) If f a (x) = ax, /^(x) = bx, a,b e R, defined on R, when are they linearly 
conjugate? 

(b) Shew that fi/ 2 and /x/ 4 are conjugate via the map h(x) = 




if x > 0 
if x<0 
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6.4 Conjugacy and the Tent Family 

We saw in Section 2.7 that for ft > (1 + y/E)/2, the tent map T M has a point of 
period three, so by Sharkovsky's Theorem, it will have points of all possible periods. 
In this section we use a certain conjugacy to show that for ft>l,T ll will have points 
of period 2" for each n > 1. Our argument is based on that in [19]. We first show 
that the interval [1/(1 + fi),n/(l + ft)] is invariant under T 2 when 1 < ft < y/2. 

The formula for T 2 in Section 2.7 gives 



fix if 0 < x < i 

— 2a 



^ ' ) fix + fi-fl 2 if I < x < l-± 



fi -fix if 1 - i < x < 1 



Proposition 6.4.1 For 1 < // < y/2, The restriction 



"' l l + i u'l + /i J l l + ft'l+fi> 
is well defined. 

Proof. Note that for this range of values of ft, 1/(1 + ft) < 1/2 and ft/ {I + ft) > 1/2, 
so that 7^(1/(1 + ft)) = ft/ (I + ft) is an eventual fixed points since T^(ft/(l + ft)) = 
fi(l + ft). 

Let x e [1/(1 + ft), ft/ {I + ft)}, then from the formula for T 2 and the fact that 

1_ _J_ _ft_ 1 
2ft < l + ft < l+ft <1 ~2jl' 

we see that on this interval the minimum value of T 2 occurs at x = 1/2 and this „• r 

^ (*) > t; 2 (i/2) = /i(i - ft/2) > ~, 

i+ ft 

since this is equivalent to 

fi - fi - 2ft + 2 < 0, or (ft- l)(fi - 2) < 0, 
where 1 < ft < \[2. 

Again assuming that x G [1/(1 + //), /x/(l + ft)], we see that on the other hand if 
x < 1/2, then T» = ftx > ft /(I + ft) > 1/2, so Tftx) = ft(l - ftx) < Mj - + 
ft)) = ft/(l + ft), so T%(x) 6 [1/(1 + ft), ft/(l + ft)}. 
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If instead x > 1/2, then 



W = „(i-,)>Mi- T ^) = T ^>|, 



SO 



1 + 1 + //' 



T£(x) = Ml -/x(l-x))=/i(l-/i + /ias) < - /i + 

so again T^x) € [1/(1 + ji), /x/(l + /i)]. □ 

We use the above proposition to show that T p and T 2 ^ are conjugate when T 2 ^ is 
restricted to a suitable invariant subinterval. 

Proposition 6.4.2 For 1 < \i < y/2, T 2 restricted to the interval 

1 /x 



i + M i + /x_ 

is conjugate to on [0, 1]. 

Proof. From Proposition 6.4.1, we see that the given interval is invariant under T 2 . 
Now we show that we actually have a linear conjugacy h(x) = ax + b: 

/ioT 2 = 7>o/i,, 

where 



1 



and 



where we can check that 



a = 



1 + // 

l + fi 



[0,1], 



1-/X 



, 5 = 



W-^— ) = 0, and h(— *— ) = 1. 

If 0 < x < 1/2, then since h~ l (x) = ar/a - 6/a, we can check that 1/2 < /i _1 (;r) < 
fj,/(l + n) < 1 - 1/2/x, so that 

/ioT 2 o/i- 1 (x) = /ioT 2 f--- > i 
M ^ \a a J 

= h (/J?(- ) +/i - /x 2 J = M x - n 2 b + a(fi - /j 2 ) + b = /x 2 x = ^(x). 
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Similarly we can check that if 1/2 < x < 1, then 1/2/i < 1/(1 + < < 1/2, 

so that 

hoT l° h~ x {x) = ho T 2 (- --) = h(fj,- fi 2 (- - -)\ 

= a// - // 2 x +^6 + 6 = /i 2 (l - x) = 7>(a;), 
i.e., in both cases we have /i o T 2 o h~ l {x) = 7>(x), giving the desired conjugacy. □ 
Corollary 6.4.3 For 1 < fi < 2, restricted to the interval 

_ fi - 1 ' /i - 1 J ' 

is conjugate to on [0, 1] . 

We apply these results to give us information about the period points of T^. 
Proof. Replace /i by ^JJi in the previous result. □ 
Theorem 6.4.4 For 1 < fi < 2, T„ has a 2 n -cycle for each n e Z+. 

Proof. We have seen that for each p > 1, has a period 2-point distinct from the 
fixed point of T„. In particular, as fj 2 > 1, 7> has a period 2-point distinct from 
the fixed point of 7>. But by the last result, 7> and T 2 are conjugate, so 7j has a 
period 2-point distinct from the fixed point of T 2 . This must be a period 4-point for 
Tp, for if not it would be a period 2-point, giving a fixed point for T* 

Continuing this argument, starting with a period 2-point for 7> and the conjugacy 
between T 2 2 and we deduce that has a period 8-point. In this way, for each 
n e Z+ we see that has a period 2 n -point. □ 

Example 6.4.5 Consider the case where fi = 2, then we see that T 2 , the standard 
tent map, is conjugate to T 2 ^ when it is restricted to the interval [\/2~ - 1, 2 - v^l . 
This implies that has the same dynamics as T 2 on this subinterval. For example, 
it must have a 3-cycle, say {ci,c 2 ,c 3 }, where the a's are distinct and T^(ci) = Ci. 
It follows that c : is a point of period 6 for T^, and in this way we deduce that 
has 2£>cycles for each k e Z + . We saw earlier that /i = (1 + y/E)/2 is where period 
3 is born for the tent family. In particular has no 3-cycle, but if a = (1 + >/5)/2 
and using the fact that Tfa (suitably restricted) is conjugate to T a , it follows that 
must have points of period 6. 
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a 2 1 a 3 1 

Remark 6.4.6 1. Suppose that a > 1 and Q < then — - — = 1 > 

l + /z 3 - 2' 1 + fi 3 1 + fi 3 ~ 

-, so that 



We see that we have a 3-cycle: 



2 S 

A* ^ 



1 + /x 3 ' 1 + // 3 ' 1 + /i 3 



! 



/i 2 1 

This 3-cycle appears when /i > 1 and j-^ — - < -, or equivalently 

/i 3 -2/i 2 + l >0, or (/z- 1)(// 2 -/i- 1) > 0. 

We see that this will happen when // > (1 + v / 5)/2. A similar analysis can be done 

a 3 1 

for other periodic orbits. For example, if a > 1 and < - we get a 4-cycle, and 

l+/i 4_ 2 to J ' 

this happens when fj 3 — fi 2 — ji — 1 > 0. 

2. Suppose that 1 < // < 2, then if £ G [/x — /i 2 /2,/i/2], we can check that T^x) G 
[// - // 2 /2, ///2], so that this interval is an invariant set. For 1 < // < v^2, the smallest 
set invariant under is a collection of subintervals of [/i - /i 2 /2,/i/2]. If /i > \/2 
this becomes all of the interval [fi - /i 2 /2, ///2], called the Ju/ia se£ of T M (named after 
one of the early pioneers of chaotic dynamics, Gaston Julia, who worked on complex 
dynamics in particular in the early 1900's). For fi = 2, the Julia set is all of [0, 1]. 
The bifurcation diagram for T^, fi > 1 gives us some insight into what is happening 
here. 

3. The conjugacy between T 2 and L 4 can be constructed by consideration of the 
periodic points of these maps. Since the period points are dense for each of these 
maps, by carefully ordering them according to their ordering in [0, 1], we can define a 
map h by defining it on the periodic points, h is then defined on a dense subset of [0, 1] , 
into a dense subset. This map can be continuously extended to a homeomorphism 
of [0, 1] with h(0) = 0, h(l) = 1. In this way it can be shown that the conjugation 
between T 2 and L 4 is unique. See the exercises for a proof that the conjugation 
between T 2 and L 4 is unique. 



6.5 Renormalization 



132 



The previous example shows that T* restricted to the interval 



fi — 1 p? — p 




[ i n 1 


y-iy-i. 







is conjugate to T^ 2 on [0, 1] (where we are replacing p, by p? in 6.4.2). How do we 
arrive at this conjugacy? Notice that for n>\,T (1 has a fixed point p^ = p,/(/j, + 1) 
and another point Vil = 1) with T^p^) = T^fa), so it is eventually fixed. Let 

us look at the graph of T 2 restricted to the interval [p^p^]. Inside the square shown 
we see that the graph obtained looks like an "upside-down" version of and we 
consider the possibility that T 2 restricted to the interval [p^pj is actually conjugate 
to T M (or in fact 7». 

Define a linear map hp : \p»,p u ] -» [0, 1] of the form h„(x) = ax + b in such a way 
that h^) = 0 and h lx {p (JL ) = \. We can check that 

MaO = T-— (*-**,), and ^(x) = (p M -p^x + p 

hp expands the interval [p M ,p M ] onto the interval [0, 1] and changes the orientation. 
This is exactly the conjugacy defined in Proposition 6.4.2. 
We define a renormalization operator of by 

(RT,){x) = h,oTloh-\x). 

What we actually showed in the previous section is that (RT^)(x) = T^(x), giving us 
the conjugacy claimed. This procedure can be continued for T*, 7j etc, and similar 
considerations can be made with the logistic map L M (see [12] for more details). 

6.6 Conjugacy and Fundamental Domains 

We have seen that two dynamical systems / and g with different dynamical proper- 
ties cannot be conjugate. On the other hand, sometimes we have dynamical systems 
having seemingly very similar dynamical properties and which we would like to show 
are conjugate. This is sometimes possible using the notion of fundamental domain, a 
set on which we construct a map h in an arbitrary manner and show that it extends 
to a conjugacy on the whole space. We first illustrate this idea with homeomorphisms 
/, g : R — ► R. We look at a fairly straightforward case where both homeomorphisms 
are order preserving and have no fixed points (in fact lie strictly above the line y = x). 

Proposition 6.6.1 Let f,g:R-*Rbe homeomorphisms satisfying f(x) > x and 
g(x) > x for all xER. Then f and g are conjugate. 
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Proof. The idea for the proof is a follows: Select x 0 G R arbitrarily and consider the 
2-sided orbit 

Of(x 0 ) = {f (x 0 ) : n G Z} = {..., X - U *o, *i, *2, . . .}. 

Since /(x) > x for all x, this is an increasing sequence: . . . x_i < x 0 < Xi < x 2 < . . . , 
so that the sets 

••.,[^-i,^o),[^o,^i),[xi,x 2 ),..., 
are disjoint and their union is all of R. We must have lim^oo x n = oo since otherwise 
the limit would exist and would have to be a fixed point. There are no fixed points 
since f(x) > x always. 

The set / = [x 0 ,/(x 0 )) = [x^x\) is called a fundamental domain for /. Set 
J = [xo, g(x 0 )) and define a map h : J -> J arbitrarily as a continuous bijection (e.g., 
we can set h(x 0 ) = x 0 and h(f(x 0 )) = g(x 0 ) and then linearly from / to J). 

Now every other orbit of / intertwines with O/(x 0 ): if y G G (x 0 ,xi), then y n = 
f n (yo) e / n (/), so lies between ^ and x i+1 . It follows that every orbit has a unique 
member in the interval [x 0 , x x ) and we use this to extend the definition of h to all of 
R. 

If x G / n (7) we define /i(x) by mapping x back to / via /~ n , then using h(f- n (x)) 
which is well defined, and then mapping back to g n (J) using gf n . I.e., if x G / n (7), 
n G Z, define 

h(x)=g n ohof- n ( x ). 
In this way, /i is defined on all of R. We can check that h is one-to-one. It is onto 
because h(f n (I)) = g n (J)) for each n, and we can check that it is continuous. Finally, 
because of the definition of /i, if x G R, then x G f n (I) for some n G Z, so x = / n (?/) 
for some y € I. Then 

p o - o h o /-(x)) - o h o - h o /(x), 

so that / and g are conjugate. □ 

Examples 6.6.2 1. The above argument can be generalized to the situation where 
f(x) and g(x) are homeomorphisms with corresponding fixed points. To be more 
precise, consider a homeomorphism / : [0, 1] -> [0, 1] which is orientation preserving, 
so that /(0) = 0, /(l) = 1 and / is increasing. Suppose that / has fixed points 
(in addition to 0 and 1), at c u c 2 , . . . , c n , then f 2 has the same collection of fixed 
points (no additional fixed points as / cannot have points of period 2 or higher). If 
f(x) > x for c k < x < Cfc+i, then we can use the argument of the proposition to 
construct a homeomorphism between / and / 2 , and do the same for each interval 
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[q, q+i] (treating the case where f(x) < x in an analogous way). In this way we see 
that / and f 2 are conjugate maps. 

2. Consider the logistic maps L„(x) = fix(l - x) for various values of /i G (0,4] and 
x G [0, 1]. We first show that for 0 < /i < A < 1, L M and L x are conjugate. There 
is a slight complication here as these maps are not increasing, but they do have an 
attracting fixed point at 0, and we saw earlier that the basin of attraction is all of 
[0, 1]. We first deal with the interval on which the maps are increasing, [0, 1/2], and 
look at the restriction of the functions to this interval. 

Our aim is to construct a homeomorphism h : [0, 1] -> [0, 1] with the property 
L x o h = h o L M . Take (L„(l/2), 1/2] - (///4, 1/2] as a fundamental domain for 
L M and (L x (1/2), 1/2] = (A/4,1/2] as a fundamental domain for L x . Define h : 
(/x/4, 1/2] (A/4, 1/2] by h{\/2) = 1/2 and /i(/i/4) = A/4 and then linearly on the 
remainder of the interval. 

Set / = (/i/4, 1/2] and J = (A/4, 1/2], then since 0 is an attracting fixed point, the 
intervals L£(J) and L n x (J) are disjoint for n G Z+, and their union is all of (0, 1/2]. 
Extend the definition of h so that it is defined on (0, 1/2] by; 

h(x) = L n x ohoL; n (x), for x G ££(/). 

We can now check that h is continuous and increasing on [0, 1 /2] when we set h(0) = 0. 

Now define h on (1/2, 1] by setting h(l - x) = 1 - h(x) for x G [0, 1/2), clearly 
giving a homeomorphism on [0, 1]. Then 

L x (h(l - x)) = L x (l - h(x)) = L x (h(x)) = h(L,(x)) = h(L,(l - *)), 
so that h is the required conjugation. □ 

3. A similar proof shows that L M and L x are conjugate whenever 1 < ji < A < 2. 
Look at the intervals [0, 1 - l//i] and [1 - l//i, 1/2] separately and the fact that 1 - l//i 
is an attracting fixed point, then use the symmetry about the point x = 1/2. 

However, these maps cannot be conjugate to L 2 since any conjugating map h : 
[0,1] -> [0,1] must have the property that h(l/2) = 1/2 (see the exercises). This 
leads to a contradiction. 

4. The maps L 4 and L M , /i G (0, 4) cannot be conjugate since L 4 : [0, 1] -> [0, 1] is an 
onto map, but L M is not (see the exercises). 



Exercises 6.6 
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1. (a) Let a,fr G (0, 1) and f a (x) = ax,f b (x) = bx be dynjamical systems on [0, 1]. 
We saw in Exercises 6.3 that these maps need not be linearly conjugate. Prove that 
f a and f b are conjugate (Hint: Use the method of examples in this section). 

(b) Let g : [0, 1] — > [0, 1] be continuous, strictly increasing with #(0) = 0 and g(x) < x 
for all x G (0, 1]. Prove that g is conjugate to f a for any a G (0, 1). 

2. Consider the tent map T^. 

(i) Show that x — 1/2 is an eventual fixed point for 

(ii) Use Section 6.4 to show that there is a subinterval of [0, 1] on which is 
conjugate to T2. 

(iii) Deduce that has periodic points of period 2k for any k > 1, but no points of 
odd period greater than 1. 

(iv) Prove that if /i > \/2, then T M has a 3-cycle. 

(v) Prove that there is an interval on which is chaotic. 

3. Let 0 < A, fx < 1. If h : [0, 1] — ► [0, 1] is an orientation preserving homeomorphism 
with h o L fi (x) = L\ o h(x) for all x G [0, 1], show that h(l/2) = 1/2 (Hint: h is a 
conjugation between two different logistic maps with h o L^{x) = L\ o h(x). Note 
that this equation also holds if we replace x by 1 — x. Use this to deduce that 
h{x) + h(l -x) = l for all x G [0, 1]). 

4. Use exercise 2 above to deduce that /i(/i/4) = A/4. It follows that the orientation 
preserving homeomorphism of Example 6.6.2 can be extended in only one way from 
[/i/4,1/2] to [A/4,1/2]. 

5*. Prove that for the logistic map L M , if 0 < /z < 2, then L M is conjugate to L£, 
the composition of with itself. Show that this is not true for /i > 2. What is the 
corresponding result for the tent family? 



6. Let S^(x) = /isin(x). Prove that if 0 < // < A < 1, then and S x are conjugate. 
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7. Prove that the rotation Ra : S 1 -> 5 1 , # a (z) = az is conjugate to the map 
T a : [0, 1) -> [0, 1), T Q (x) = x + a (mod 1), when a = e 2nia . Can Z? a be conjugate to 
Rb for a ^ b? 



8. Prove that T Q : [0, 1) -* [0, 1), T a (x) = x + a (mod 1) is conjugate to its inverse 
map T~ x (x) = x-a (mod 1). Can T a be conjugate to T*l 

9. The aim of this exercise is to show the uniqueness of the conjugacy between the 
tent map T 2 and the logistic map L 4 . 

(i) Check that this conjugacy k : [0, 1] -> [0, 1] is given by, 

2 

k(x) = - arcsin(y / x); T 2 ok = koL 4 . 



(ii) Suppose that h : [0, 1] -+ [0, 1] is another conjugacy between T 2 and L A , then 
h(0) = 0, h(l) = 1 and h is a strictly increasing continuous function (why?), i.e., h is 
an orientation preserving homeomorphism of [0, 1] . 

(iii) Show that h maps the local maxima (respectively minima) of T£ to local maxima 
(respectively minima) of L\. 

(iv) Use the fact that any such conjugation is order preserving to show that h(x) = 
k(x) at all local maxima and local minima. 

(v) Use the continuity of h and k to deduce that h(x) = k(x) for all x e [0, 1]. 

(vi) Deduce that there is no C 1 conjugacy between T 2 and L A . 



10. Use the above exercise to show that if L 4 (x) = 4x(l - x), then L\ has turning 
points at sin 2 (kn/2 n+1 ), for k = 1, 2, ... , 2" +1 - 1. 

11. Use the fact that the conjugation between T 2 and L 4 is unique to show that if 
4> : [0, 1] -» [0, 1] is a homeomorphism satisfying L 4 o 0 = 0 o L 4 , then <j>[x) = x for all 
x E [0, 1], i.e., 0 is the identity map (hint: first show that k o 0 is also a conjugation 
between T 2 and L 4 ). 
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12*. Let be the tent map. Show that if \i > y/2, then for each open interval 
U C [0, 1], there exists n > 0 such that 

[22(1/2), T„(l/2)] C T2(U). 

(Hint: Use the fact that |T M (J7)| > /jl\U\ if U does not contain 1/2, so that the length 
keeps increasing. We claim there exists m > 0 such that T m (U) and T m+1 (U) both 
contain 1/2, for if not, \T m+2 (U)\ > n 2 \U\/2 for all m G Z + , a contradiction, since 
this eventually exceeds 1). 
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Chapter 7. Singer's Theorem 

We shall now show that the logistic map L M : [0, 1] -» [0, 1], L^x) = (xx(l - x) 
with 0 < fi < 4 has a t most one attracting cycle . We use a result due to Singer (1978) 
which is applicable to maps having a negative Schwarzian derivative. Recall that a 
map / : K -» R is C 3 if f"'(x) exists and is continuous. 

7.1 The Schwarzian Derivative Revisited 

Recall the Schwarzian derivative of f(x) is: 



Sf(x) = 



f'"(x) 



/"(*) 



1 2 



= n*) - \ [F(x)f , 



f'(x) 2 [f'(x))\ 
where F(x) = 

Our first goal is to show that many polynomials have negative Schwarzian deriva- 
tives. 

Lemma 7.1.1 Let f(x) be a polynomial of degree n for which all the roots of its 
derivative f'(x) are distinct and real. Then Sf(x) < 0 for all x. 

Proof. Suppose that the derivative of f(x) is given by 

f'(x) = a(x - n)(x - r 2 ) • • • (x - r n _i), 

where a G R, then 

f"(x) n ~ l i 

and so 

n-1 

F'(x) = -Y— - . 

tr(*-n) 2 

Now put this into the Schwarzian derivative: 

Sf(x) = F\x) - \[F{x)f < 0. 

Lemma 7.1.2 Assume that f is a C 3 map on R, then 

(i) S(f o g)(x) = Sf{g(x)) ■ {g'{x)f + Sg(x). 

(ii) IfSf<0 and Sg < 0, then S(f o g) < 0. 



□ 
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(iii) IfSf < 0, then Sf k < 0 for- all keZ + . 

Proof, (i) As above we have F(x) = Q^-, so set G(x) = and = ^ 

JK X ) g( x ) h'(x) 

where /i = fog. Then 

= /'(</(*)) • g'(x), h"{x) = (f"(g(x))g'(x) 2 + (f'(g(x))g"(x), 

so that 

H(x) = f"(9(x)W(x)* + f(g(x))g"(x)) 
1 ' ' f'(9(x))g'(x) 



This gives 
so that 



f"(g(x))g'(x) g"(x) 

/'(<?(*)) W) = F{9{x))g {x) + G{x) - 

H'(x) = F\g{x))g'{xf + F(g(x))g"(x) + G'(x), 



F'( 9 (x)) - -F\g(x)) 



S{fog){ x ) = H\x)- l -H\x) 

g'(x) 2 + F(g(x))g"(x)-F(g(x))g'(x)G(x)+G'(x)~G 2 (x) 



= Sf(g(x)).g'(x) 2 + Sg(x), 

since G(x) = g"{x)/g'{x). 

(ii) is now immediate, and (iii) follows using induction. □ 

Example 7.1.3 Let g(x) = aX + ^ a, 6, c, d G R, a linear fractional transformation. 

cx + a 

A direct calculation shows that Sg(x) = 0 everywhere in its domain. It now follows 
from Lemma 1 that if h(x) = g(f(x)) = then = 

We now prove a version of Singer's Theorem. Recall that if / : R -> R is a 
continuous function and c G R is an attracting fixed point or attracting periodic 
point, then the basin of attraction B f (c) is an open set. Denote by W the maximal 
open interval contained in B f (c) which contains c (called the immediate basin of 
attraction of c) . 

Theorem 7.1.4 Let f : R — > R be a C 3 map with negative Schwarzian derivative. If 
c is an attracting periodic point for f, then either: 

(i) the immediate basin of attraction of c extends to oo or —oc, or 
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(ii) there is a critical point of f (i.e., a root of f(x) = 0), whose orbit is attracted to 
the orbit of c under f. 

Proof. We first look at the case where c is a fixed point of /. Suppose that its 
i immediate basin of attraction is the open interval W and that (i) does not hold, then 
W is a bounded set. 

Thus W = (a, b) for some a, b G R. Because of the continuity of / and the fact 
that /(a), f(b) (a, 6), there are three possibilities for f(a) and f(b). 

Case 1: f(a) = f(b) - this will happen for example if a is a fixed point and b is an 
eventual fixed point. It follows from the Mean Value Theorem that (a, b) contains a 
critical point of /. 

Case 2: f(a) = a and f(b) = 6, then by the Mean Value Theorem, there are points 
xi e (a,c) and x 2 G (c,6) such that f'(xi) = f'(x 2 ) = 1 (see picture). But since c 
is an attracting fixed point, |/'(c)| < 1. It follows that either f'(x 0 ) = 0 for some 
x 0 G or has a minimum value /' Oo) > 0. In the latter case we have 

f (x Q ) > 0, f"(x 0 ) = 0 and f'"(x 0 ) > 0, so that Sf(x 0 ) > 0, contradicting the 
Schwarzian derivative being everywhere negative. 

Case 3: f(a) = b and f(b) = a. Here c is fixed by f 2 and so are a and 6, so that (f 2 )' 
has a zero x 0 in (a, 6) (as in Case 1). But 

(/ 2 )'(*o) = or tf w)/w = ° 5 

so either x 0 or /(x 0 ) is a root of but both lie in (a, 6). 

Now suppose that c is a point of period k, then f k (c) = c, an attracting fixed point 
for /*, then from our earlier arguments, the immediate basin of attraction of c (for 
f k ) contains a critical point of f k , say x 0 : 

(f k )'(x 0 ) = f({x 0 )f'(f(x 0 )) ■ ■ ■ fif- 1 (x Q )) = 0, 

so that f{f m (x Q )) = 0 for some 0 < m < k. In this case f m (x 0 ) € f m (W) C W, the 
basin of attraction of c. □ 

Example 7.1.5 Consider the map f(x) = x - x 5 . We see that f'(x) has two real 
roots: ±(l/5) 1 / 4 , and f"(x) and f'"{x) are both continuous. We have f"{x) = -20x 3 , 
and f"'(x) = -60x 2 . 
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Substituting these into the Schwarzian derivative gives: 

Sf(x) = (-60z 2 )/(l - 5x 4 ) - 3/2[(-20:r 3 )/(l - 5a: 4 )] 2 = + 5x% 

and this is always negative . We can check that the critical points are in the basin of 
attraction of the fixed point x — 0. 

Example 7.1.6 The map f(x) = 3x/4 + x 3 cannot have a negative Schwarzian 
derivative. It has fixed points x = 0 and x = ±1/2, 0 being attracting with bounded 
basin of attraction, but / has no critical points. 

7.2 Singer's Theorem 

Corollary 7.2.1 Let f : [0, 1] -> [0, 1] be a C 3 map with Sf(x) < 0 for all x. The 

basin of attraction of an attracting cycle con tains 0, 1 or a critical point of f(x). 

Proof. If J = (a, 6), 0 < a < b < 1, is the basin of attraction of an attracting 
cycle, then we have seen above that it must contain a critical point of /. Any other 
attracting cycles will be of the form [0, a) or (6, 1], so will contain 0 or 1. □ 

Example 7.2.2 We now see that the logistic map = fj,x(l - x) } 0 < fi < 4, 

x e [0, 1], has at most one attracting periodic cycle. If 0 < \i < 1, 0 is the only 
attracting fixed point, having basin of attractions [0,1]. For 1 < /x < 4, has 
exactly one critical point x 0 = 1/2. 

Since 1^(0) = fi > 1, the fixed point 0 is unstable; therefore [0, a) cannot be 
a basin of attraction. Furthermore, L^(l) = 0 and hence (6, 1] is not a basin of 
attraction either. Since SL^x) < 0 everywhere (at x = 1/2, L^x) = — oo), 

we conclude that there is at most one attracting periodic cycle in (0, 1) and the result 
follows. 

If we look at the bifurcation diagram of the logistic map, it follows, for example 
that where we see six horizontal lines, we have an attracting 6-cycle, and not two 
attracting 3-cycles. 

Remark 7.2.3 In a similar way we get Singers Theorem [45], proved by David Singer 
in 1978. This theorem is actually a real v ersion of a theorem about complex polyno- 
mials proved by the French mathematician Gaston Julia in 1918 [22]: 

Oj (2) (X I 

Suppose that f is a C 3 map on a ^closed ii^ejyal I such that Sf(x) < 0, for all ] 

x e I. Iff has n critical points in I, then f has at most n + 2 attracting cycles. 

LJh^: Ooi*t P is bo^> ^ pm*s* *W a OnboA foKb. 
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t^X^J^ Example 7.2.4 Let G{x) = Aarctan(^), A^O. Then G'(x) = A/(l + x 2 ). Clearly 
] Yv=o G(x) ha s no critica l points. Now, if |A| < 1, then c = 0 is an asymptotically stable 
S(^c^_ g xec l p 0 i nt with basin of attraction (-00,00). If A > 1, then G has two attracting 

fix ed points X\ and x 2 with basins of attraction (-00, 0) and (0, 00), respectively. 
beeaurO^e Finally, if A < -1, then G has an attracting 2-cycle {x u x 2 } with basin of attraction 



Exercises 7.2 



1. If G(x) = Aarctan(x), (A / 0), show that the Schwarzian derivative is 



SG(x) = 



-2 



(1 + x 2 ) 2 ' 

Use this to verify the conclusion from Example 7.2.4. 




2 L (i+xT-2i 



ir 
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Iff 



^ T^^wvrotl System AAJz^ 



Stability of Two- Dimensional Maps 

Is evolution a matter of survival of the fittest or survival of the 
most stable? 

A. M. Waldrop 

4.1 Linear Maps vs. Linear Systems 0u+«tz --tv(A)=^\+ / )vz 

Recall from linear algebra that a map L : R' 2 —> R 2 is called a linear transfor- 
mation if 

1. L(Ui + U 2 ) = L{U X ) + L(U 2 ) for UuU 2 € R 2 

2. L(aU) = aL(U) for U € R 2 and a € R. 

Moreover, it is always possible to representj^with a given basis for R 2 ) by 
a matrix A. A typical example is 



L 



x\ i ax + by 
y J \cx + dy 



which may be written in the form 



L 



x 



x \ _ fa b 



or 

W) = AU, (4.1) 

where£/=(*) aod ^ tA"*C\ *l V^jfcxW* 

By iterating L, we conclude that £ n (c/) = A n t/. Hence, the orbit of U 
under /'is given by — 

L {U,AU,A 2 U,....A n U,...} (4.2) 

Thus, to compute the orbit of U, it suffices to compute A n U for n <E Z + . 
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172 ^(q^Oi^d -jpC^O ^< J Discrete Chaos 

Another way of looking at the same problem is by considering the following 
two-dimensional system of difference equations 

x(n + l)= ax(n) + by(n) 
y(n + l) = cx{n) +dy(n), 

or 

U{n+l) = AU(n). (4.4) 
By iteration, one may show that the solution of Equation (4.4) is given by 

U(n) = A n U(Q). (4-5) 

So. if we let U 0 = U(0), then L n {U Q ) = t/(n). 

The form of Equation (4.3) is more convenient when we are considering 
applications in biology, engineering, economics, and so forth. For example, 
x(n) and y(n) may represent the population sizes at time period n of two 
competitive cooperative species, or preys and predators. cr-v^\ x^^' 

In the next section, we will develop the necessary machinery to compute A n 
for any matrix of order two. The general theory may be found in [32, 33, 60]. 



4.2 Computing A n 

Consider a matrix A = (a 0 ) of order 2x2. Then, p(A) = det(A - XI) is called 
the characteristic polynomial of A and its zeros are called the eigenvalues 

of A. Associated with each eigenvalue A of A a nonzero eigenvector V € R 2 
with AV = XV. 



Example 4.1 

Find the eigenvalues and the eigenvectors of the matrix 



A = 



23' 
1 4 



D 



SOLUTION First we find the eigenvalues of A by solving the character- 
istic equation det(A - XI) = 0 or 



2 - A 3 
1 4- A 



= 0 



which is 



A 2 - 6A + 5 = 0. 
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Hence, Ai = 1 and A2 = 5. To find the corresponding eigenvector Vi, we solve 
the vector equation AV\ = AVi or (A - X\I)V\ = 0. 
For Ai = 1, we have 



13\/vii\ (0 



Hence, v u + 8^21 = 0. Thus, v u = -3v2i. So, if we let ^21 = 1, then 
t>ii = -3. It follows that the eigenvector V\ corresponding to A] is given by 

;-.')■ 

For A2 = 5, the corresponding eigenvector may be found by solving the 
equation (A - A 2 /) T /2 - 0. This yields 



Vi = 





Thus, -3i>i2+3i>22 = 0 or v\2 = ^22- It is then appropriate to let v\2 = ^22 = 1 
and hence V2 = y j ^ . I 

To find the general form for A n for a general matrix A is a formidable task 
even for a 2 x 2 matrix such as in Example 4.1. Fortunately, however, we 
may be able to transform a matrix A to another simpler matrix B whose nth 
power B n can easily be computed. The essence of this process is captured in 
the following definition. 



DEFINITION 4.1 The matrices A and B are said to be similar if there 
exists a nonsingulai 1 matrix P such that 

/UpBPJ = p-*ap = b. ^ a?=PB 

We note here that the relation "similarity" between matrices is an equiva- 
lence relation, i.e., 

1. A is similar to A. (ref/c^WrLy) <\>-^&o\ 

2. If A is similar to B then B is similar to A. ^^y^eWj) ^Jf^ 

3. If A is similar to B and B is similar to C, then A is similar to C. tfctr&M 

# ft g ^ C ^ ^ C (£Y*ns rt *Wtj) 

The most important feature of similar matrices, however, is that they possess 
the same eigenvalues. 



( invert 



{hie) 



1 A matrix P is said to be nonsingular if its inverse P~ l exists. This is equivalent to saying 
that det P ^ 0, >vhere det denotes determinant. 
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T'J-fVOTt.El^I 4 1 

Let A and B be two «m«or matrices. Then A and B have the same e,gen- 
values. 

PROOF Suppose that P" ' AP = B or A = PBP^ Let A be an 
eigenvalue of A and V be the corresponding eigenvector. Then, Af - ay 
PBP 'V. Hence, B(P-'V) = A(P">V). Consequently, A is an e.genvalue of 
B with P' l V as the corresponding eigenvector. I 

The notion of similarity between matrices corresponds to Hnearcoj^y, 
which we have enured in Chapter 3. In other words, two hnear maps 
conjugate if their corresponding matrix representations are similar Thus he 
linear maps UM on R 2 are linearly conjugate if there ex.sts an mverfble 
map h such that 

L\ o h — h o L-2 

h ~ o L\° h = L 2 . 

The next theorem tells us that there arc three simple -canonical" forms for 
2x2 matrices. 



x'f real matrix. Ifen -4 to similar to one of the following 

matrices: 



2. 
3. 



PROOF Suppose that the eigenvah.es A, and Aj are eal. Jhen, we 
1 have two cases to consider. The first case is where A, + A 2 . In ttaa** 
3 may easily show that the corresponding eigenvectors V, ami V s ar Imea y 
^ independent (Problem 10). Hence, the matm P = (ViM), re., ^ matr,^ 
-1 P whose columns are these eigenvectors, is nonsingular. Let f 



1 (ef 



Then, 



W AP=PJ. ^ 

Comparing both sides of Equation (4.6), we obtain 

AVi = eVi + gVl 
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Hence, 

A1K1 = eV l +gV 2 . 

Thus, e — Ai and g = 0. 

Similarly, one may show that / = 0 and /? = A 2 . Consequently, J is a 
diagonal matrix of the form (a). 

The second case is where Ai = A2 = A. There are two subcases to consider 
here. The first subcase occurs if we are able to find two linearly independent 
eigenvectors V\ and V2 corresponding to the eigenvalue A. This subcase is 
then reduced to the preceding case. We note here that this scenario happens 

when (A - \I)V = 0 for all V G IR 2 . In particular, one may let V\ = 

and V 2 — ( j V which are clearly linearly independent. 

The second subcase occurs when there exists a nonzero vector V2 G M 2 such 
that (A - XI)V2 / 0. Equivalently, we are able to find only one eigenvector 
(not counting multiples) V\ with (A - \I)V\ = 0. In practice, we find V 2 by 
solving the equation 

(A-\I)V 2 - V v 

The vector V2 is called a gcenerali zed eigenvector of A. Note that AV\ = AVi 



and AV2 = AV2 + Vi. Now, we let 



and P- 1 AP = J. Then, 



Comparing both sides of Equation (4.7) yields 



(47) 



J =(11). (4.8) 

The matrix J is in a Jordan form. 

Next, we assume that A has a complex eigenvalue Ai = a + i(3. Since A is 
assumed to be real, it follows that the second eigenvalue A2 is a conjugate of 
Ai, that is, A 2 = a - ifi. Let V — V\ + iV 2 be the eigenvector corresponding 
to A]. Then, 

AV = AiV 
A(Vi + ivy = (<* + i/3)(Vi + Wt). 

Hence, 

AVj = aVi PV 2 

AV 2 = 0V l +aV 2 , 

letting P = (Vi, V 2 ) we get P~ 1 AP = J. Hence, 




4P = PJ. 



(4.9) 
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Comparison of both sides of Equation (4.9) yields 

J = 



a 0 
-0 a 



(4.10) 



Theorem 4. 2 gives us a simple method of computing the general form of A n 
for any 2 x 2 real matrix. In the first case, when P~ l AP = D = ( ^ ® V 
we have 



A n = (PDP- l ) n 
= PD n P' 1 



(4.11) 




In the second case, when P~ l AP = J = 

A n = PJ n P~ l 



A 1 

0 A 



then 



fi" = P 



A n nA"- 1 
0 A" 



P 



- i 



(4.12) 



Equation (4.12) may be easily proved by mathematical induction (Problem 11). 



In the third case, we have P~ l AP = J = 



Q 0 



. Let u> = arctan (0/oc). 



Then smw== /?/|Ai|. Now, we write the matrix J in the form 

J = \\i\ 



<*/|Ai| j9/|Ai| 
-/?/|Ai| o/|Ai| 



cosu; sinu 
- sinu; cos u 



By mathematical induction one may show that (Problem 11) 

cos nu sin nu \ 



J n = |Ai| n 



~ Sill nu> COS 7iU> 



and thus 



I 



cos nu> sin nu; 
- sin nu) cos nu; 



P 



3 



c 



Example 4.2 

Solve the system of difference equations 



X(n+l) = AX(n) 



(4.13) 
(4.14) 




(4.15) 



where 
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SOLUTION The eigenvalues of A are repeated: Ai = A 2 = 2. The only 
eigenvector that we are able to find is V\ = y 3 Y To construct P we need to 
find a generalized eigenvector V2. This is accomplished by solving the equation 

(A - 21) V2 = V\ . Then, V2 may be taken as any vector ( ^ V with 3y - 2x = 

1. We take V 2 = ( j J. Now if we put P = ( jj j J, then P~MP - J - 



2 1 

0 2 




. Thus, the solution of Equation (4.15) is given by 

[3 l\ (T n2 n ~ 1 \ ( 1 
- 1,2 1^0 2» J 1,-2 3 J ^0 



REMARK 4.1 If a map / : R 2 - K 2 is given by /(A" 0 ) = AXq, then 
/ n (X 0 ) - 4 u *o - PJ n P' l X 0 . In particular, if Xq = (J J, then / n (X 0 ) = 

2 n f 1 " 371 ) for all n e Z 4 . I 



(fry X< ^ <J 



Exercises - (4.1 and 4.2) 

In Problems 1-5, find the eigenvalues and eigenvectors of the matrix A and 
compute A n . 



1. A = 



-4.5 5 \ Rew^S, , v 



4.5 -1 



2. A= 7 tU- 



2.25 1.5 



3. ,4 = 



f 8/3 1/3 \ n-^-bcUc rtj 

l " 4/34/3j ® if ^ 

/ 2 3\ c *' 

4 • ^ = ( -3 2 ) T-tr« *rCfl) =: c*+d**n*« 

117 -a trfW-vii , 



5. a= 
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6. Let L : R 2 — > R 2 be defined by L(X) = AX, where A is as in Problem 1. 

'!)• 



Find L n 



7. Solve the difference equation X(n + 1) = AX(n). where A is as in 
Problem 3, and X(0) = ^ j J . 

8. Solve the difference equation X(n + 1) = /4A"(n). where A is as in 
Problem 4, and X(0) = X 0 . 



9. Let / : R 2 -» R 2 be defined by f{X) = with .4 as in Problem 5. 

'!)■ 



Find f n 



10. Let A be a 2 x 2 matrix with distinct real eigenvalues. Show that the 
corresponding eigenvectors of A are linearly independent. 

1 1 • (a) If J = ( q I ) , show that ./" = ( A 0 " nA A 'l" 1 ) . 



(b) If J = ( - J q ) ' show that J = l A ' 



sin iiuj cos no; 



where 



|A| - /a 2 +/i 2 ,o;-arctan(f). 
12. Let a matrix A be in the form 

0 1 )■ 

V ~P2 -Pi J 

(a) Show that if A has distinct eigenvalues Ai and A2, then 



— =( A 0A° 2 ). 



where P — 



a, a 2 )- 

(b) Show that if A has a repeated eigenvalue A, then 

Al\ 
OA/' 



P~ l AP = 



where P — 



1 0 
A 1 



(c) Show that if A has complex eigenvalues Ai = a+ifl and A 2 = a—ifJ, 
then 



/ 1 n \ 

where P = 



Stability of Two- Dimensional Maps 



179 



4.3 Fundamental Set of Solutions 

Consider the linear system 

X{n + l) = AX{n), (4-16) 

where A is a 2 x 2 matrix. Then, two solutions Xi(n) and X 2 (n) of Equation 
(4.16) are said to be linearly independent if X 2 {n) is not a scaler multiple of 
Xi{n) for all n G Z + . In other words, if ciXj (n) + c 2 X 2 (n) = 0 for all n € Z + , 
then ci = c 2 = 0. A set of two linearly independent solutions {Xi(n), X 2 (n)} 
is called a fundamental set of solutions of Equation (4.16). 



DEFINITION 4.2 Let {X\ (n), X 2 (n)} be a fundamental set of solutiom 
of Equation (4-16). Then 



is c 



X(n) = fciXi(n) + fc 2 X 2 (n), fci,fc 2 € R .. (4.17) 



caMed a general solution of Equation (4-16). 



Finding Xj (n) and X 2 (n) is generally an easy task. We now give an explicit 
derivation. 

In the sequel Ai , A 2 denote the eigenvalues of A: \\ , V 2 are the corresponding 
eigenvectors of A. 

We have three cases to consider. 7 
C ase (i) 

Suppose that P~ l AP = J = f* 1 £ J • Then a general solution may be given 
bv 



X(n) = A n X(Q) = P^P" 1 *^) 



k 2 



where 



- p-'xio) 



Then, 




X(n) = fciA?Vi+fc 2 A 2 % 



(4.18) 



Here, Xi(n) = A?Vi and X 2 {n) = AJ> l V 2 constitute a fundamental set of 
solutions since in this case V\ and V 2 are linearly independent eigenvectors. 
Note that one may check directly that A?Vi and X 2 l V 2 are indeed solutions of 
Equation (4.16) (Problem 13a). 
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'A 1 
0 A 



Suppose that P~ x AP = J = 
given by 

x{n) = p^p-^ro) 



Then, a general solution may be 





^A n nA n_1 ^ 
V 0 A " > 








\MnJ = A;iA n Vi+A:2(nA n - 1 Ki 


+ \ n V 2 ) 



(4.19) 

Hence. *i(n) - A^Vj and X 2 (n) = X n V 2 + n\ n ~ l Vi constitute a funda- 
mental set of solutions of Equation (4.16) (Problem 13b). 



Case (iii) 

Suppose that P~ l AP = J 
solution may be given by 



ft 0 
-6 ft 



If w= arctan(/?/a), then the general 



X(n) = PJ n P~ l X(Q) 
= (V r iK 2 )|A 1 | n [ 



cos no; sin nuj 

- sin no; cos nu> 



X^V= |Ai|"[&] cosnw -T^sinnu;)Vi~~J 
+ sin no; + & 2 cos no;) V^], |f 



fc 2 



(4.20) 



Hence, X^n) = |Ai|"[($cosnu;)Vi - (^sin(nu;))K 2 ] and X 2 (n) 
= (l^il n [($ sin(nw))Vi + (^cos(nw)]V2 constitute a fundamental set of solu- 
tions (Problem 13c). 

Example 4.3 

Solve the system of difference equations 




X(n + 1) = AX(n), X(0)= (lY 



w 



here 



A = 



-2 -3 
3 -2 



0 



SOLUTION The eigenvalues of A are Ai = -2+ 3i and A 2 = — 2— 3i. The 
corresponding eigenvectors are V = ^ ~ 1 j and V = f ~* V respectively. 

This time, we take a short cut and use Equation (4.20). The vectors 
V\ and V 2 referred to in this formula are the real part of V, V\ = ( ~ 1 J , and 
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123.69°. Thus. 



X{n) = (13) n / 2 



(k\ cos nu> + k'2 sin nto) 



+ (—k] sin nu) + k2 cos nto) 



X(0) = 



] 



= fa 



0 



+ A: s 



Hence, fci = l,fc 2 = 2. Thus, 



-I 



A» = (13) n/2 (cos + 2 sin nu>) ^ + (- sinnw + 2 cos no;) ^ 



/inui/2 (—cosruv - 2 sin no; \ 
I - sin + 2 cos nu> / ' 



4.4 Second-Order Difference Equations 

A second-order difference equation with constant coefficients is a scalar equa- 
tion of the form 

u(n + 2) + pyu{n + 1) + p 2 u(n) = 0 (4.21) 

Although one may solve this equation directly, it is sometimes beneficial to 
convert it to a two-dimensional system. The trick is to let u(n) = X\{n) and 
u(n + 1) = x 2 (n). 

Then we have 

xi(n + 1) = x 2 (n) 

x 2 (n + l) = -p 2 x\{n) -pix 2 (n) 

which is of the form 

X(n+l) = AX(n) (4.22) 

where 

X(n)=(*f\), and ° 1 V 

The characteristic equation of A is given by 

A 2 +P1A+P2 - 0. (4.23) 

Observe that we may obtain the characteristic Equation (4.23) by letting 
u(n) = A" in Equation (4.21). Thus, if X 1 and A 2 are the roots of Equation 
(4.23), then m(n) = A? and u 2 {n) = A£ are solutions of Equation (4.21). 

Using Eqs. (4.18), (4.19), and (4.20), we can make the following conclusions: 
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1. If Ai ^ A 2 and both are real, then the general solution of Equation 
(4.21) is given by 

(4.24) 



u(n) = ciA" + C2A2, 



1 



2. If A 1 = A-2 = A. then the general solution of Equation (4.21) is given by 

(4.25) 



u(n) = c\X n + c 2 nA n , 



3. If A] = a + i A2 = o -i/3, then the general solution of Equation (4.21) 

3 



is given by 



where 



arc t an {/3/a 



(4.26) 



Example 4.4 

Solve the second-order difference equation 



*(// + 2) + 6x(n + 1) + 9a:(/z) = 0, x(0) = i, x(l) = 0. D 

SOLUTION The characteristic equation associated with the equation is 
given by A 2 + (iA + 9 = 0. 

Hence, the characteristic roots are Ai = A 2 — -3. The general solution is 
given by 

x(n) = #-3) n + c 2 n(-3) n 

x(0) = 1 = ci 

x(l) = 0 = -3ci - 3c 2 . 

Thus, C2 = — 1 and, consequently, 

x(n) = (-3) n -n(-3) a 
= (-3) n (l-n) I 



Exercises - (4.3 and 4.4) 



1 Solve the system 



xi(n+ 1) 
X2{n + 1) 
with x\ (0) 



-xi(n) + ar 2 (n) 

2rr 2 (n) 

1, * 2 (0) = 2: 



Stability of Two- Dimensional Maps 183 

2. Find the general solution of the system 

X(n-U) = AX{n), where A = 

3. Solve the problem X(n + 1) = AX(n), where 4 = 

4. Solve the system 



X(ri + 1) = AX{n). with /l = 




5. Solve the system x(n + 1) = Ar(n), with A = 

6. Solve the difference equation 

x(n + 2) - bx{n + l) 6 x{n) = 0, x{0) = 2 

(a) By converting it to a system, 

(b) Directly as it is. 

7. Solve the difference equation 

F(n t- 2) = F(n + 1) + F(n), F(l) - 1, F(2) - 1. 
(This is called the Fibonacci sequence.) 

8. The Chebyshev polynomials of the first and second kind are defined as 
follows: 

T Tl (x) = cosfncos-^x)), 

U n {x) = / — _ sin[(n4- ljcos" 1 ^)], for \x\ < 1. 
Vl - x l 

(a) Show that T n {x) satisfies the difference equation 

T n+2 (x) - 2xT n+l {x) + T n (x) = 0. T 0 (x) = 1, T x {x) = x. 

(b) Solve for T n {x). 

(c) Show that 

U n+2 (x) - 2xU n+1 (x) + Un{x) = 0, U 0 (x) = 1, Ui(x) = 2x. 

(d) Write down the first four terms of r n (a;) and U n {x). 

9. Solve the equation x(n + 2) + 16.x(n) = 0. 
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10. Let A be a 2 x 2 real matrix with distinct eigenvalues Ai and A 2 . Prove 
that the corresponding eigenvectors V\ and V 2 are linearly independent. 

11. Let A be a 2 x 2 real matrix with a repeated eigenvalue A. Let V\ be an 
eigenvector corresponding to A and let V 2 be a generalized eigenvector. 
Show that V\ and V 2 are linearly independent. 

12. Let A be a 2 x 2 real matrix with complex eigenvalues \\ = a + i(5 and 
A 2 =a-i/3. Suppose that V = Vi f / V' 2 is the eigenvector corresponding 
to Ai. Prove that the matrix P = (Iq.VA) is nonsingular. {Hint: It 
suffices to show that Vi and V 2 arc li nearly independent.) 

13. (a) Show that Xi(n) and X 2 (n), obtained from Equation (4.18), are 

solutions of Equation (4.16). 

(b) Show that X\(n) and X 2 {n), obtained from Equation (4.19), are 
solutions of Equation (4.16). 

(c) Show that X\(n) and X 2 (n), obtained from Equation (4.20), are 
solutions of Equation (4.16). 

In Problems 14 and 15. consider the iionliomogeneoiis equation 

Y{n + l) = AY(n) + g(n) (4.27) 

where A is a 2 x 2 matrix and g is a function defined on Z + . 

14. Show that 

H - i 

Y{n) = A n Y(0) + ^A n ~ k ~ l g{k). (4.28) 

(This is called the variation of constants formula.) 

15. Use Formula (4.28) to find the solution of Equation (4.27) with 

16. Solve the equation y(n + 2) - 5y(n 4 I) + 4y(n) = 4 n . 



, 4.5 Phase Space Diagrams fi^ i s^>^^ >u ^ 

, c ^ C One of the best graphical methods to illustrate the various notions of stability 
QO^ X is the phase portrait or the phase space diagram. Let / : K 2 -» K 2 be a 

given map. Then, starting from an initial point X Q = f^o))' W6 P ^ 0t ^ 
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sequence of point Xq, f{Xo), f 2 {Xo)< f 3 (Xo), . . . and then connect the points 
by straight lines. An arrow is placed on these connecting lines to indicate 
the direction of the motion on the orbit. In many instances, we need to be 
prudent in choosing our initial points in order to get a better phase portrait. 

In this section, we consider linear systems for which f(X) = AX, where A 
is a 2x2 matrix. Observe that if A - I is nonsingular, i.e., det(A-I) ^ 0, then 

the origin ^jjj is the only fixed point of the map /. Equivently, X* = 

is the only fixed point of the system 

X(n+1) - AX(n) (4.29) 

As stipulated in Theorem 4.2, there exists a nonsingular matrix P such that 
P~ l AP = J where J is one of the forms (1), (2), or (3) in Theorem 4.2. If 
we let , 

X(n) = PY(n) 0 | (4.30) 

in Equation (4.29), we obtain ><^+0= P^-r p~'><( 

K(n + 1) = jy(n). 3r (4.31) 

Our plan here is to draw the phase portrait of Equation (4.31), then use 
the transformation (4.30) to obtain the phase portrait of the original sys- 
tem (4.29). 

(I). We begin our discussion by assuming that J is in the diagonal form 

J = ( ^ f ], where A] and A2 are not necessarily distinct. Here we have 
V U A 2 / 

two linearly independent solutions: 

Y 1 (n) = X{Vu and Y 2 (n) = A£K 2 , where VJ = ^Jj and K 2 = ^Jj 

are the eigenvectors of j4 corresponding to the eigenvalues Ai and A 2 , respec- 
tively. 

Observe that Y\ (n) is a multiple of V\ , and thus must stay on the line ema- 
nating from the origin in the direction of V\\ in this case, the x axis. Similarly, 
y 2 (n) must stay on the line passing through the origin and in the direction 
of V2; in this case, the y axis. These tw r o solutions are called straight line 
solutions. The general solution is given by 

W-Mi (J), ™>)=(£)- ( 4 - 32 ) 

We have the following cases to consider: 

1. If I Ai| < 1 and |A 2 | < 1, then all solutions tend to the origin as n — > 00. 
Observe that if |Ai| < |A 2 | < 1, then, |A"| goes to zero faster than |A 2 |. 
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FIGURE 4.1a 

(a) A sink: 0 < Aj < A 2 . 
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FIGURE 4.2a 

(a) A source: Ai > A2 > 1. 



And consequently, any solution Y(n) in the form (4.31) is asymptotic 
to the straight line solution Y 2 {n) = A 2 (see Fig. 4.1a). 

On the other hand, if |A]| > |A 2 |, then Y{n) is asymptotic to Y\{n) = 
A?(j) (see Fig. 4.1b). 

Phase portraits 4.1a and 4.1b are called sinks. 

2. If |Ai| > 1, and |A 2 | > 1, then we obtain an source as illustrated in 
Figs. 4.2a and 4.2b. 

Note that if |A X | > |A 2 | > 1, then Y(n) is asymptotic to Y 2 {n) = A^Jj 

(the y axis) when n -* -00 and is dominated by Yi(n) = A" f * J when 
n — > oc. 7 

3. If fAii < 1 and |A 2 | > 1, then we obtain a saddle (Fig. 4.3). In this case, 
when n — 00 Y(n), is asymptotic to Y 2 (n) = A 2 as n -> 00 and is 

asymptotic to Yi(n) = A" ^ * J as n -» -oc. Similar analysis is readily 
available for the case |Ai| > 1 and |A 2 | < 1. 

4. If Ai = A 2 , then 
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FIGURE 4.3 

A saddle: 0 < Ai < 1,A 2 > 1. 




FIGURE 4.4 

A sink: 0 < Ai < A 2 < 1. 



Hence, every solution Y(n) lies on a line passing through the origin with 
a slope k 2 /ki (see Fig. 4.4). 

Observe that in each of the four subcases, the presence of a negative eigen- 
value will cause the solution Y(n) to oscillate around the origin and the phase 
portrait will not look as nice as in Figs. 4.1a 4.4. 

(II). Suppose that J is in the form 

'-(»)■ 

Then, we only have one straight line solution, Y\{n) = A n V'i = A n 
The general solution is given by 

r(n) = fc 1 A n (j)+A; 2 (nA^ 1 (j) + A" ( J)) 

Now, if |A| < 1, then,F(n) -> 0 as n -> oo, since lim nA n-1 = 0 (by 
L'Hopital Rule). Since the term h\ n (®J tends to the origin, as n oo, 

faster than the term (ki\+k 2 n)X n ( J J, our solution Y{n) tends to the origin 
asymptotic to the x axis (see Fig. 4.5a). In this case, the origin is called a 




FIGURE 4.5a 

(a) A degenerate sink: Ai = X- 2 = A. 0 < A < 1. 



y 2 




*1 




FIGURE 4.5b 

(b) A degenerate source: Ai = A 2 = A, A > J . 
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degenerate sink. Figure 4.5b depicts the case when |A| > 1 and in this case, 
the origin is called a degenerate source. 
(III). Suppose that J is in the form 

'-(■«)■ 

In this case, we have no straight line solutions due to the presence of cosn/5 
and sinn/3 in the solutions Y\{n) = |Ai | n (fci cosn/5 + /^sinn/^) f^J and 

y 2 (n) = |Ai| n (-A;i sinn/3-f focosn/?) f j J- The general solution is given by 

, / fci cos n0 + /J2 sin n/3 \ 

r[n)-\Ai\ y_ klSinn p + k2COSnj3 j 

with Y(0) = ^ y- Define an angle 7 by setting COS7 = fci/ro, and 

sin 7 = ^ 2 /r 0 . where r 0 = y/fc? + k%. Then 

j/i(n) = |Ai| n r 0 cos(nu; - 7) 
y . 2 (n) = -|Ai| n r 0 sin(7iu> - 7). 

Thus, the solution in polar coordinates is given by 



r(«) = \]y 2 M) +yj{n) 

= r 0 \\i\ n (4.33) 

\yi(n)/ 

= -(nu--7) (4.34) 
It follows from Eqs. (4.33) and (4.34) that 

1. If |Ai| < 1, then we have a stable focus where each orbit spirals to- 
ward the origin [Fig. 4.6(a)]. On the other hand, if |Ai| > 1, then we 
have an unstable focus, where each orbit spirals away from the origin 
[Fig. 4.6(b)]. 

2. If I Ai| = 1, then we have a center, where the orbits follow a circular 
path [Fig. 4.6(c)]. This is due to the fact that yf(n) + y$(n) = rj. 

To this end. we have obtained the phase portraits of Equation (4.31), which 
may be called "canonical"' phase portraits. To obtain the phase portraits of 
the original system (4.29), we apply (4.30), i.e., we apply P to the orbits 

of Equation (4.31). Since P (^j = (V u V 2 ) (J) = V u and P ( J J = 
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(Vi, V ^[°ij = aPP^K P t0 the orbits Y ^ amounts t0 rotatin S the 
coordinates; the x axis to V\ and the y axis to V 2 . In other words, the straight 
line solutions are now along the eigenvectors \\ and V- L Using this observa- 
tion, one may opt to sketch the phase portrait of Equation (4.29) directly 
and without going through the canonical forms. 

The set of points on the line emanating from the origin along V\ is called 
the stable subspace W s ; the set of points on the line passing through the 
origin in the direction of V 2 is called the unstable subspace W n . Hence, 

W" = {X € R 2 : A n X 0 as n - oo}, (4.35) 
W u = [X e R 2 : A~ n X - 0 as n - oo} (4.36) 

The following example illustrates the above-described direct method to 
sketch the phase portrait. 



Example 4.5 

Sketch the phase portrait of the system X(n + 1) = AX{n). where 

H.25!) D 

SOLUTION The eigenvalues of A are Ai = §, and A 2 = 5; the corre- 
sponding eigenvectors are V\ = ^ ^ and V 2 = f _ 2 j J , respectively. Hence, 

we have two straight line solutions, X\(n) = (1.5) n ^ J and X 2 (n) 

= (0.5) n f J^)- The general solution is given by X(n) = /ex (1.5)" + 

fc 2 (0.5) n Note that x ( n ) is asymptotic to the line through V x = 

(see Fig. 4.7). I 



4.6 Stability Notions 

Consider the map / : R 2 -» R 2 and let X* = J be a fixed point of /; i.e., 
f(X*) = X*. 

Our main objective in this section is to introduce the main stability notions 
pertaining to the fixed point x*. Observe that these notions were previously 




FIGURE 4.6 

(a) Stable focus: |Ai| < 1. (b) Source: |Ai| > 1. (c) Center: Ai, 2 = Q ± i/3, 
|Ai l2 | = l. 
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FIGURE 4.7 

Saddle: Ai > 1, 0 < A 2 < 1. Stable and unstable subspaces: W\ W u . 
introduced in Chapter 1. The only difference in R 2 is that we replace the 



absolute value by a convenient "norm on 



2 . Roughly speaking, a norm 



of a vector (point) in K' 2 is a measure of its inagnit 
follows. 



ude. A formal definition 



DEFINITION 4.3 A real valued function on a vector space V is said to 
be a, norm onV , denoted by | |, if the following properties hold: 

1. \X\ > 0 and \X\ = 0 if and only if X = 0, for X € V. 

2. \aX\ = |q||X| for X GV and any scalar a. 

3. \X + Y\ < \X\ + \Y\ for X, Y 6 V (the triangle inequality). 

In the sequel, we choose the t x norm on R 2 defined for X = ^ 

\X\ = \xi\ + \xi\ 



as 



(4.37) 



For each vector norm on R 2 there corresponds a norm || || on all 2 x 2 
matrices A = {ay) defined as follows 

(4.38) 



||>l||=sup{|i4X|:|X|< I}- 

It may be easily shown that for A' 6 E 2 , 

\AX\ < \\A\\ \X\ 



(4.39) 
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Let p{A) be the spectral radius of A defined as p(A) = max{|Ai|, |A 2 | : 
Ai, A 2 are the eigenvalues of A}. Then it may be shown that for our selected 
vector norm 

= max{(|o u | + |a 2 i|), (|oi 2 | + |a 22 |)}. (4-40) 

(For a proof see [49].) 

For example, if X = Q), |X| = 3. And for the matrix A = (^_ 2 _ 4 j, 

= max{3,7} = 7. The eigenvalues of A are A] = -2, A 2 = -1. Note 
that p(A) = max{| - 2|,| - 1|} = 2, and thus < PHl 

It is left to the reader to prove, in general, that p{A) < \\A\U for any matrix 
A (Problem 14). 

Without any further delay, we now give the required stability definitions. 



DEFINITION 4.4 A fixed point X* of a map f : K 2 -> R 2 is said to be 

1. stable if given e > 0 there exists 6 > 0 such \X - X m \ < 6 implies 
\f n {X) -X*\ <e for all n G Z + (see Fig. 4.8a). 

2. attracting (sink) if there exists v > 0 such that \X - X*\ < v implies 

lim f n (X) = X*. It is globally attracting if v = oc (see Fig. 19). 

n—>oo 

3 asymptotically stable if it is both stable and attracting. It is globally 
asymptotically stable if it is both stable and globally attracting, (see 
Fig. 112(a)) 

4. unstable if it is not stable (see Fig. 4.8b). 



REMARK 4.2 In [91], Sedaghat showed that a globally attracting fixed 
point of a continuous one-dimensional map must be stable. Kenneth Palmer 
pointed out to me that this result may be found in the book of Block and 
Coppel [12] Moreover, the proof in Block and Coppel requires only local 
attraction (see Appendix for a proof). The situation changes dramatically in 
two- or higher dimensional continuous maps, for there are continuous maps 
that possess a globally attracting unstable fixed point . We are going to present 
one of these maps and put several others as problems for you to investigate. 
I 



Example 4.6 

Consider the two-dimensional map in polar coordinates 
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FIGURE 4.8b 

(b) X* = 0 is an unstable fixed point. 
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FIGURE 4.9 

x * = ( J J is an unstable globally attracting fixed point. 



Then, 



9 



9) " 1 (27t)( 1 " r ">& 2 ' " 



Clearly lim g n (j) = LM = (j)- Thus, each orbit is attracted to the 



J Y However, if 0 = Six, 0 < (5 < 1, then the orbit of ( JJ 1 will 



fixed point 

spiral clockwise around the fixed point (jj before converging to it. Hence, 
f is globally attracting but not asymptotically stable (see Fig. 4.9). D 



4.7 Stability of Linear Systems 

In this section, we focus our attention on linear maps where /(X) = AX, and 
A is a 2 x 2 matrix. Equivalently, we are interested in the difference equation 

X(n + 1) = AX(n). (4.41) 
For such linear maps, we can provide complete information about the sta- 
bility of the fixed point X* = (z). The main result now follows. 
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THEOREM 4.3 

The following statements hold for Equation (4-4 V : 

(a) If p(A) < 1, then the origin is asymptotically stable. 

(b) If p(A) > 1, then the origin is unstable. 

(c) If p{A) == 1, then the origin is unstable if the Jordan foivn is of the form 
, and stable otherwise. 

PROOF Suppose that p(A) < L Then it follows from Eqs. (4.18), (4.19), 

(4.20). that Inn X(n) = 0. Thus, the origin is (globally) attracting. To prove 
stability, we consider three cases. 

(i) Suppose that the solution J<(n) is given by Equation (4.18). This is 
the case when the eigenvalues of the matrix A are real and there are two 
linearly independent eigenvectors. 

X <" ) = P (o f A°) P "' X(0) ' 

Hence, 

\X(n)\ < \\P\\ \\P- 1 \\ P (A)\X(0)\ 
< M\X(0)\ 

where M = \\P\\ \\p- l \\ p(A). 

Now, given e > 0, let S - e/M. Then |X(0)| < 5 implies that \X(n)\ < 
Mb = e. This shows that the origin is stable. 

(ii) Suppose that the solution X(n) is given by Equation (4.19). This 
case occurs if the matrix A has a repeated eigenvalue A and only one 
eigenvector. 

^(r,)=p( A 0 nnA A ';r l )p-'x(o) 

l^(n)l<li/ 5 lll|P^ 1 ||(«|A|"- 1 + |Ar)|A(0)|. 

Since n|A| n ~> 0 as n oo, there exists iV e Z + , such that the term 
(n|A| n_1 + |A| n ) is bounded by a positive number L. Hence, 

\X(n)\ < M\X{0)\ 

where M = L\\P\\ \\P~ l \\. 

The proof of the stability of the origin may be completed by seting 
S = e/M as in part (a). 
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(iii) If A is not in the form * J . then it is either diagonalizable to J = 

V where |Ai| < 1 and |A 2 | - 1 or |A 2 | < 1 and |Ai| = 1. In either 
case, J n is bounded and hence the origin is stable. I 
The proofs of parts (b) and (c) are left to you as Problem 11. 



Exercises - (4.5 - 4.7) p 

In Problems 1-9, sketch the phase portrait of the system X(n+ 1) = AX{n), 
when A is the given matrix. Determine the stability of the origin. 



1. 



2. 



3. 
4. 
a. 
6. 
7. 
8. 
9. 
10. 
11. 



1/2 0 \ 
0 1/2) 

1/5 (A 
0 2) 

02) 

-1/2 1 
, 0 -1/2. 

' 0.5 0.25 N 
-0.25 0.5 , 

8 ; 

^ 1 3\ 

:i9 



2 ^ * Discrete Chaos 

12. Show that if X* is a fixed point of a linear map / on R 2 and is asymp- 
totically stable, then it must be globally asymptotically stable. 

13. Complete the proof of Theorem 4.3. 

14. Prove that for any 2 x 2 matrix A,p(A) < \\A\\. 



4.8 The Trace-Determinant Plane 
4.8.1 Stability Analysis 

Table (4.1) provides a partial summary of everything we have done so far 
In this section we provide another important way of presenting these results,' 

, Til T determma i nt P lan ?' wher « we employ pictures, rather than a 
able, rhis turns out to be a better scheme when one is interested in studying 
bifurcation in two-dimensional systems 

The following two results provide the framework for using the trace-deter- 
minant plane. Recall that for matrix A = fa «») , tr A = a „ + ^ and 
det A = a u a 22 + a V2 a 2x . 

^ THEOREM 4.4 

Let A = (a tj ) bea2x'2 matrix. Then p(A) < 1 if an d only if 



\tr A\ - 1 < det A < 1. 



(4.42) 



PROOF 



(i) Assume that p{A) < 1. If Al> A 2 are the eigenvalues of A, then JAJ < 1 
Ttf A \ K n x2 characteristic equation of the matrix A is given by 

7 i V ~^ 11+a 22M-(a„a 22 -a 12 a 21 ) = 0,orA 2 -(^4)A + 
det A = Q. Hence the eigenvalues are 



Ai - - [rr- A + /{tr A) 2 - 4 det J] , 
A 2 = ~ [tr A - /(tr~ A) 2 - 4 det ij . 



Case (a) Al and A 2 are real roots, i.e., (tr A) 2 -AdetA>Q. Now -1< Ai, 

A 2 < 1 implies that 



"2 < tr A ± y/(tr A) 2 - A det A < 2 
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TABLE 4.1 

Partial summary of the stability of linear systems. 
Type Eigenvalue 



201 



Saddle 



Sink 



Source 



Center 



Phase Portrait 



0 < Ai < 1 < A 2 




0 < A 2 < Ai < 1 




A 2 > Ai > 1 






Spiral Sink \ = a±i0, \X\ < 1, 0 f 0 




Spiral Source \ = a±iP, |A| >_L£jjO 



X = a±i(3, |A| = l t 0?Q 




"Oscillatory" Saddle -1< A| < 0, A 2 < -1 





"Oscillatory" Source Ai > 1, A 2 > -1 
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or 



-2 - tr A < y/(tr A) 2 - 4 det A < 2 - tr A (4.43) 
-Z- tr A < -^(tr A) 2 - 4 det A < 2 - tr A. (4.44) 
Squaring the second inequality (4.43) yields 

1 - tr A 4 det A > 0. (4.45) 

Similarly, if we square the first inequality in (4.44) we obtain 

1 + tr A + det A > 0> (4.46) 

Now from the second inequality (4.43) and the first inequality in 
( 1.44) we obtain 2 4 tr A > 0 and 2 - tr A > 0 or \tr A\ < 2. Since 
(tr A) 2 - 4 det A > 0, it follows that 

det A < (tr A) 2 / 4 < 1. (4.47) 

Combining (4.45), (4.46) and (4.47) yields (4.42). 
Case (b) Ai and A 2 are complex conjugates, i.e., 

[tr A} 2 - 4 det A < 0. (4.48) 

In this case we have A li2 = i j/r ,1 ± i^/4detA-(tr A) 2 } and 

l A i2 ,2 (tr A) 2 4 deZ A (tr A) 2 

|Ai| =|A 2 | =— — 4 — i—J-^rfeiA 

4 4 4 

Hence 0 < det A < 1. To show that inequalities (4.45) and (4.46) 
hold, note that since det A > 0 either (4.45) (if tr A > 0) or 
(4.46) (if A < 0) holds. Without loss of generality, as sume t hat 
tr A > 0. Then (4.45) holds. From (4.48), tr A < 2^deTA. If 
we let det A = x, then x € (0, 1) and f(x) = 1 4 x - 2^ < 
1 + detA-tr A. Note that /(0) = 1 and /'(0) = 1 - indicate 
that / is decreasing on (0, 1) with range (0, 1). This implies that 
1 + det A - tr A > f(x) > 0 and this completes the proof. 

(ii) Conversely, assume that (4.42) holds. Then we have two cases to con- 
sider. 

Case (a) (tr A) 2 - 4 det A > 0. Then 

|Ai, 2 | = - tr A ± y/(tr A) 2 - 4 det A 



1 



< g \ tr A ± V(det A + l) 2 -4det l| 

< i (de< A 4 1 4 x/(rfei /I - l) 2 ) 
= ^ (dei A+l - (det, A - 1)) = 1. 
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Case (b) {tr A) 2 - 4 det A < 0. Then 



|Ai, 2 | - ^ tr A ± iy/4 det A - {tr A) 2 
= - 1 y/{tr A) 2 + 4 det A- {tr A) 2 
= Vdet A < 1. 



I 

As a by-product of the preceding result, we obtain the following criterion 
for asymptotic stability. 



COROLLARY 4.1 

The origin in Equation (4.41) is asymptotically stable if and only if condition 
(4.42) holds true. 

Note that condition (4.42) may be spelled out in the following three in- 
equalities: 

I + tr A + det A > 0, or 



det A > -tr A - 1 (4.45)' 



1 - tr A + det A > 0, or 



det A > tr A - 1 (4.46)' 



det A < 1. (4.47)' 



(4.49) 



(4.50) 
(4.51) 



Viewing det A as a function of tr A, the above three inequalities give us the 
stability region as the interior of the triangle bounded by the lines det A = 
-tr A - 1, det A - tr A - 1, and det A = 1, as shown in Fig. 4.10. 

Next we delve a little deeper into finding the exact values of the eigenvalues 
of the matrix A along the boundaries of the triangle enclosing the stability 
region. The following result provides us with the needed answers. Let Ai = 
\ (tr A + yjltr A) 2 - 4 det A^j, A 2 = \ (tr A - yf(pr A) 2 - 4 det l) be the 
eigenvalues of A. 



THEOREM 4.5 

The following statements hold for ami 2 x 2 matrix A. 
(i) If \tr A\ - 1 = det A, then we have 

(a) the eigenvalues of A are Ai = 1 and A 2 = det AiftrA>0, 

(b) the eigenvalues of A are A 2 = -1 and A] = -det A if tr A < 0,. 
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(ii) If\tr A\-l< det A, and det A = 1, then the eigenvalues of A are e 
where 0 = cos~ ] {tr A/2). 



±i6 



PROOF 



(i) Let \tr A\-l = det A. Then (tr A) 2 - 4 det A = (det A + l) 2 > 0. 
This implies that the eigenvalues are real numbers. Moreover, A w - 
\(tr A± /(tr Af - 4 det A) = \ (tr A ± (det A - 1))- 

(a) If tr A > 0, then tr A~l + det A, and consequently. 



Al,2 = 



det A 



(b) If tr A < 0, then tr A = -1 - det A, and consequently, 

Al,2 = 



1 

det A 



(ii) Let \tr A\-l< det A, and det A = 1. Then (tr Af - 4 deM < 
(dei A + l) 2 —4 = 0. Hence, the eigenvalues are complex conjugates. 
Moreover, 

Xi,2 = \ ((tr Af/A ± JitetT^WW) 
= \tr A±sJ\-(tr Af. 

Thus |A,, 2 | = s/(tr A)V4 + 1 - (tr .4)^/4 = !• Furthermore, fl = arctan(A u ) 

= tan -r WIX^ 11 ) m C0S ' i{tr m < WWCh giVS AW = C± ' e = C0S " * 
isin#. 



4.8.2 Navigating the Trace-Determinant Plane 

The trace-determinant plane is effective in the study of linear systems with 
parameters. It provides us a chart of those locations where we can expert 
dramatic changes in the phase portrait. There are three critical Jo . L *T 
denotes the trace and D denote the determinant. Then there are thre cr*^ 
lines: D = tr A _ i, £> = _ tr A - 1, and D = 1; they enclose the stability 
region in the trace-determinant planes. 

We now illustrate our analysis by the following example. 
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detA = 1/4(tr A) 
I detA =J^T) 



X : =-detA, 

K =-1 




< det A = tr A - 1 

A., =1, \ 2 =detA 



complex eigenvalue 



tr A »-t -T 

real eigenvalue 



FIGURE 4.10a 

(a) The stability region for Equation (4.41) is the shaded triangle. 

|Tl-l <0<1 

ItI-i < o oU 

lTl<D+l 
-D-l <T <D*l 
ftA) <l^^j^ ^i^^? |p>^T-\j °U [ P>T-ll «Ki [b<lj 

o< A. < 1, K > 1 



X 2 <-1, -1 <x 1 < 0 




complex eigenvalue 

X,> 1, 0<X 2 <1 
tr A 



X,> 1, X 2 <-1 



FIGURE 4.10b 

(b) The determination of eigenvalues in different regions. 
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det A 




oscillatory 
source 



FIGURE 4.10c 

(c) Description of the dynamics of Equation (4.41) in all the regions in the 
del-trace plane. 



Example 4.7 

Consider the one-parameter family of linear systems X(n + 1) = AX{n) t 
where 




which depends on the parameter a. As a varies, the determinant of the matrix, 
det A, is always 2a - 1, while the trace of the matrix, tr A, is always 0. As 
we vary the parameter a from negative to positive values, the corresponding 
point (T,D) moves vertically along the line T = 0. Now if D < -1, which 
occurs if 2a - 1 < -1 or a < 0, we have a degenerate case, Ai = 1 and 

A 2 = -1 with corresponding eigenvectors and ^jj. Thus every point 

on the y-axis is a fixed point and every other point in the plane is periodic of 
period 2. For 0 < a < 5, we have a sink, and for \ < a < 1 we have a spiral 
sink. At exactly a = 1 we have a center, and if a > 1 we have a spiral source 
(see Fig 4.10b). 

The values of a where critical dynamical changes occur are called bifurcation 
values. In this example, the bifurcation values of a are 0, 5, 1. D 



or 



Exercises - 4.8 

In problems 1-6, consider the one-parameter families of linear systems 
X(n + 1) = AX(n), where A depends on a parameter a. In a brief essay, 
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discuss different types of dynamical behavior exhibited by the systems as a 
increases along the real line, modeled after Example 4.7. 



1. A = 

2. A = 

3. A = 
A. A - 
5. A = 



' 2 3\ 
v l + o4j 

'3 + a 2\ 
s -2 3 J 

'2 + a A 
v 0 2) 

/a+1 1 \ 
a a 4- 1 ) 

'a + ltt^a' 
1 a + 1 



g A= /" + lv/T^,_l<a<l 



In problems 7-9, we consider two-parameter families of the linear system 
X(n + 1) = AX{n), where 4 depends on two parameters a. 3. 

In the afc-plane, identify all regions where the system possesses a saddle, a 
sink, a spiral sink, and so on. 



7.* A - 



f Q + 1 l\ 

K 0 2) 

'a + 1 0 
k 5 a + 1 



8.* 4 : 

In the trace-determinant plane (Fig. 4.10d). 

1. Show that we have a saddle in regions in (3) and (5). 

2. Show that we have an oscillatory saddle in regions (7) and (8). 

3. (a) Show that we have an oscillatory source in region (6). 
(b) Show that we have a spiral source in region (3). 



X4.9 Liapunov Functions for Nonlinear Maps 

In 1892, the Russian mathematician A. M. Liapunov (sometimes transliter- 
ated as Lyapunov) introduced a new method to investigate the stability of 
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Exercises - (4.1 and 4.2) 

1. eigenvector (\) is associated with Ai = 3, (J is associated with 



.2, 



/_9 . Qi j. A 0 . 1 

\ _ 1 ATI _ * ° + on * O 2 "- 1 
A 2 - 2 , A — I _on-t-l , _3_ on+1 _ __L 



3. Aj = A 2 = 2 
/ 1 \ 



V 2 



-2) 

'2"(n+3) 



4" — 3 3 

n2 n + 1 2n(3n-l) 

3 3 



5. A, = A 2 = =i^Si 

2 2 
] 



7. X(n) 



-v^sin^ +cos^ -2>/3sm ^ 
^sin 2 ^ coB2p + vfi8hi2gs 

/ 2 n ~ 1 (n + 2)\ 
, 2 n (l-n) J 



9 r»/^-W°W "^sin^ \ 



Exercises - (4.3 and 4.4) 



3. A(n)— ^ i2 n + 2A>23 nJ 

5 *7nl-~ ^i(0) + 3- 1 x 2 (0)\ 
5> X(n) ~ V 3^ 2 (0) ) 

7. X(n) = ki cos n6 + k 2 smnd where 0 = tan -1 (4) « 76 c 
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f i n 

15. K(n) == I 2 „ _ j 



Exercises - (4.5 - 4.7) 

1. The origin is asymptotically stable since p(A) = \ < 1 (Theorem 4.13). 




Ai = A2 — \ 
Phase portrait: origin is asymptotically stable. 

3. Since p(A) = 2 > 1, it follows by Theorem 4.13 that the origin is asymp- 
totically stable. 










-^^> 



Aj = A2 = 2 
Phase portrait: origin is unstable. 
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5. Eigenvalues of A are Ai = ^ + Jt, A 2 = 5 - j*- 



ia.i-m + ^5-t <1 

Hence by Theorem 4.12, the origin is asymptotically stable. 




A stable focus. 



7. The eigenvalues of A are Ai = 2, A 2 - 3. The corresponding eigenvectors 



are V\ = f _ 5 j, V 2 



-3 



, straight-line solutions are 



The origin is unstable since p(A) = 3 > 1. 




Ai = 2, A 2 - 3 
Unstable node. 



Answers 
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9. The eigenvalues of A are Ai = i, A 2 = -i. The origin is stable but not 
asymptotically stable. 




Exercises - 4.8 

1. (i) a > 4 : oscillatory source 

(ii) a = 4 : unstable origin (bifurcation value) 

(iii) a < 4 : a saddle 

4 

(iv) a = — : unstable origin 

3 

(v) a < -- : spiral sources 

3. Bifurcation value: a = -3 
a < -3 : a saddle 
a = -3 : unstable oscillatory 

7 _ 2 

5. (i) a >= — where the eigenvalues are A 2 = - 1 and Ai = - det A = - 

6 7 
and we have an ^dilatory" stable origin, for a < --, we have an 

oscillatory saddle. 

(ii) a = -1, the eigenvalues are Ai = A 2 = 0 and every point is a fixed 

point, for -- < a < -1, we have a sink, and for -1 < a < 0, we 

have a spiral sink. 



Answers 



(iii) a = 0, where the eigenvalues are Ai = 1, A 2 = 1, where we have 
eigher a stable or unstable origin, for a > 0, we have a saddle. 

7. Region 5: a < b and 3a + 6 > b 
Region 6: a < b and 3a + 6 < 6 
Region 7: a > b and 3a + 6 < b 

Region 1: b > ---(a - I) 2 
4 

Region 8: 6 < 3a + 6 and 6 < 2a + 1 

Region 4: b < a and b > --(a - I) 2 and b < 2a -r 1 



X Exercises - (4.9) 

L Let V(xi,x 2 ) = #1 f 

3. Let V == X\Xo< 
5. Let - x 2 + 4y 2 . 

7. Let V(x,y) = x 2 + y 2 . 
11. Let V(x, y) = ary. 



Exercises - (4.10 and 4.11) 

3. |a| < 1 



5. a/3 < 1 





7. |6j < 1 — a < 2 



